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MELBOURNE 


THE CONSTANT FLUX PROBLEM IN NON-UNIFORM 
EXPONENTIAL MEDIA 


By P. R. Witson* 
[Manuscript received April 4, 1960] 


Summary 


A description is given of radiative transfer, under conditions of constant flux, in a 
semi-infinite medium in which the attenuation coefficient increases exponentially with 
depth, and has small sinusoidal variations in a direction parallel to the surface. 


The directional intensity at the surface is evaluated numerically in several cases, and 
the extrapolation to other cases is discussed. 


I. INTRODUCTION 

In astrophysical problems the equation of transfer has usually been applied 
to plane parallel semi-infinite media, uniform except for possible variation with 
depth. Recently, however, investigations of the granulation in the solar 
photosphere and of the structure of the umbra and penumbra of sunspots have 
focused attention on the need for solutions for non-uniform media. Giovanelli 
(1959) investigated the problem in which the attenuation coefficient, scattering 
parameter, and source function are independent of depth, but exhibit small 
sinusoidal variations as functions of one coordinate parallel with the surface. 

In any discussion of stellar atmospheres, however, it is necessary to consider 
the extent to which the increase of attenuation coefficient with depth affects 
the solution of the equation of transfer. In the present paper the solution of the 
constant flux problem is undertaken for a plane parallel semi-infinite medium in 
which the attenuation coefficient, in addition to small sinusoidal variation 
parallel to the surface, increases exponentially with depth. 

The transfer equation is solved here in terms of the total intensity, which is 
then integrated numerically to obtain the directional intensity at the surface. 


II. THE EQUATION OF TRANSFER 
Consider an isotropic medium in which x is the attenuation coefficient, xA 
the absorption coefficient, x@, the coefficient of single scattering, and S, the 
source function, is the ratio of emission per unit volume and solid angle to the 
attenuation coefficient. For a non-uniform medium, x, A, @, and S are scalar 
functions of position. The intensity of radiation, J, is a function of position and 
direction, while the total intensity J, defined by 


Ts | 7aQ, 
4m 


* Department of Applied Mathematics, University of Sydney. 
+ GIoVANELLI, R. G. (1959).—-Aust. J. Phys. 12: 164. 


is a function of position only. 


A 
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From the transfer equation, Giovanelli (1959) obtained an equation for J 
embodying an approximation similar to that of Eddington, 


VA =o VI. + 1(aT 408). eee (1) 
If the medium is in radiative equilibrium, then 
AnrS=dJ, 
and (1) simplifies to 
VI =_Vd.Vu AO re (2) 


The following discussion is confined to this conservative case. 


III. CHARACTERISTICS OF THE MEDIUM 
If the attenuation coefficient of a semi-infinite medium increases 
exponentially with depth below the surface and varies sinusoidally in one 
direction parallel to the surface, x must have the form 


ne (La-o Con te), Gea eect cee (3) 


the z-axis being normal to the surface z=0 and positive outwards. The simple 
model of a stellar atmosphere, in which the attenuation coefficient tapers off 
exponentially to a negligibly small value, corresponds to the limiting case where 
% tends to zero. A non-zero value of x9, with x—0 for z>0, provides a simple 
model for an atmosphere whose attenuation coefficient initially increases rapidly 
with depth and then more slowly according to (3). 

Since the horizontal structure size is determined by 27/1, and the scale 
height by v—1, the appearance of the medium will be found to depend largely on 
whether v/l is large or small as compared with unity. 


IV. SOLUTION OF THE EQUATION OF TRANSFER IN THE 
CONSERVATIVE CASE 


To solve (2) for a medium in radiative equilibrium carrying constant mean 
flux, J is expressed in the form 


J (7,2) =2j,(2) cos nla. 


Here jo(z) is the mean value of the total intensity at depth 2, while j,(z) is the 
amplitude of the variations in J(«,z) having the same periodicity as the structure 
of the medium. For n>2 the functions j,(z) represent distortions in the intensity 
distribution, and if « is sufficiently small they may be disregarded. Then 


oJ (w,2)o29,(2) -1-9,(2) CON Ut amc ee eee (4) 
where j,/jo is assumed small. Substituting (4) with (3) into (2) yields 


jo(1 +-« cos lx) +-(j1 -2j,){cos la +4a(1 +cos 2lx)} =j,.4012(1 —cos 21x) 
—v{jo+4ofi +(j1 + ajo) cos la+daj; cos Qla}, 6... cece eee (5) 
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where jo=dj,/de and j”=d?j,/de?. Equating the coefficients of cos nla when 
n=0 and n=1 yields 


Most WO) tees eye 204,) 0,0 isa (6) 
and 
a(jo-+vjo) +(jitvyii—l2j,)=0. oo. eee ee (7) 
The coefficients of cos 2l” correspond to the distortion terms and may be ignored 
for sufficiently small «. 


Elimination of j) and use of the approximation «?=0 yield a differential 
equation for 7,, having the solution 


4,—C,e"* +D,e%, CO ic Ooo: Ciabeo 0 Onrney 00 (8) 


where p =3{ —v+ /(v? +4l*)}, 9 = —3{v+ /(v?+41?)}, and C, and D, are constants 
of integration. Substitution in (6) or (7) yields the solution for j,, 


jo =Cp t+ Doe” +ha(CyeP?+Dye%), occ ce dec eecees (9) 
where O, and D, are also constants of integration. Thus 
J (a@,2) =C,) +Doe—” + (Ce? +D,e”)(4a+cos lx). .... (10) 


The constants of integration Oj, 0;, D), and D, are determined by the 
constant flux requirement and the boundary conditions. Deep within the 
medium the following conditions should hold. 

(i) The mean intensity j, may increase linearly with optical depth but at 
no greater rate, aS’ may be inferred by comparison with known solutions for 
uniform media. 

(ii) The form of J(a,z) should be self preserving, i.e. a jiljg Should 


exist and, from (4), should be small. As z+—oo, the mean ‘optical depth 
behaves as e~”. However, since lim e”/e—” is infinite, condition (i) requires. 


D, —()) F . 
Alternatively, from (6) and (7) 
a Pa pao Dye” 
hal a Eo ae 4¢D,e?+D e-” 


Since « is small, condition (ii) also requires D,=0. 
At the surface z=0, the inward radiation flux is assumed zero. As usual 
in the treatment of the plane parallel case, this leads to the condition 
xJ (7,0) = —3dJ (x,0)/dz. 


Substituting for J at 20, cos la=-+1, we obtain Cy and 0, in terms of Do, 
whence finally 


TOE to ea ae a pcoslut. .... (11 
Flare) =Dy) 37 +6 2—] ae rae COs 4 iss 8 


Here also terms of order «? have been ignored. 


464 P. R. WILSON 
\ The flux per unit area across any layer parallel to the surface is 


1 dJ(x,z2) 


ers hg 
Deep within the medium this has the value 


vDy 
3% (1+ cos lax) 


and, if the net flux is independent of x) and v, (11) takes the form 


J (x,2) =A} $+ (e-*--1) eP? cos la ; (12) 


a 
3/2-+p/[% ¥ 
where the constant A replaces D,v/x) and is determined solely by the mean value 
of the net flux. In particular, the coutrast at the surface, 7,(0)/j)(0), given by 
(12), is ne 
dr lig= —&Xo| (Xp +#P)- Bile, BO ww wt ae wy aw Lie Se im oa (13) 


Thus an intensity maximum is associated with an attenuation minimum. 


It is of interest to note that in the limit as v0 (i.e. if the medium is uniform 
with depth), then p—l and 


j = .— — Yeo. lz 
os J (&,2) AN Xp 3D Ux cos wt. 


This is similar in form to Giovanelli’s solution for a medium which is uniform with 
depth. 


V. THE DIRECTIONAL INTENSITY AT THE SURFACE 
If 0 is the angle between the emergent radiation and the z-axis, and ¢ is 
the azimuthal angle referred to the w-axis, the directional intensity I(#,z,6,0) 
at the surface (2=0) is found by integrating J(#,z) along a path cutting the surface 
at #, in a direction given by 8 and @. The position coordinate s of any point on 
the path measures the distance from the suriace to that point along the path, 
and is taken as positive below the surface. This yields the relations 


z—=—s cos 8, ) 
Aw = —s sin 6 cos 9. J 


In a volume of the path ds at s, the radiation scattered or emitted into unit 
solid angle in the direction of the path is (1/47)x(s)J(s)ds, where from (3) and (14) 


%(S) =x ge © 971 + cos I(a—es)], .......... (15) 


ce representing sin§cos@. Let @,(t) be the fraction of the radiation which 
has suffered no further scattering or absorption at a distance t from the surface. 
Then in a further distance —dt, 


d®,(t) = —O,(t)x(t)(—dt).. 
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After integrating this equation the fraction remaining at the surface, ®,(0), 
is found to be 


® (0) =exp (; ae plit +« cos ) cos(la +))} —evs 08 91 + & cos Y cos(la —les +h) 


where tan )=(l/v) tan 0. Thus the directional intensity at the surface is 


al co 
I(x,0,0,¢) =F DAOC IIS 08st Peeeett eee (17) 
0 
where 
Tis)=A} 3 +o" cos 0 __].) — 3)2 he e—Ps cos 8 egg Hoes) mea 1S) 
= 0 


from (14) and (11). 

It will be noted from (17) and (11) that the profiles of Z(#,0,0,¢) and J(a,2) 
depend only on the values chosen for the ratios x,/v, v/l and on the value of a. 
Although the structure size of each profile depends on the actual values of J and v, 
every position coordinate in (17) and (11) is multiplied by / or v (or p, where 
p=t{—v++1/(v?+4l*)}). Thus by specifying only the ratios v/l and x,/l, and 
evaluating I at various values of lw, (17) may be applied to media of any structure 
size merely by using the appropriate value of 1. 


VI. CALCULATIONS AND DISCUSSION 
In the special case of the normal intensity at the surface, 9=o=0 and 
(17) becomes 


Ay iN A a cos. le ~& v —ply 
= = e ed hae alt d 
B2),2,8) rat 1+. cos Z| At fees : (Fee iy 
» pee aR (19) 


where x’=x,(1+« cos lx) and 7(s)==(x’/v)(e%—1). The integral in (19) may 
be related to the incomplete gamma function. In three cases approximate 
forms are obtained. 


(i) If v/ils1, —p/v~ —(l/v)?~0, and 


5 A 
SS — — 1 Ch moet Deca Decapoda dc 20 
I~s, = a cos 1x) (20) 


(ii) If v/l<1, —p/v~—I/v~—oo, and 


b A 3 
P= Ap Nn ee ene € My | Rares Cech ers corsa 21 
Ix~5 =(! 52 008 Lt) (21) 


(iii) If x,/l<1, «/(3/2+p/%)<1, and 


Se Te te iol eee 22 
I~ Af 54 008 L). Deas enigh (22) 
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It is useful to denote the ratio of the fractional variation in the surface intensity, 
AI/I,, to the fractional #-variation in attenuation coefficient, «, by the symbol £. 
Thus @ describes the relative contrast at the surface. In case (1) above, B=1 d 0, 
while in (ii) and (iii), 80-6. The case (iii) (%)/1<1) is important In astrophysics 
and it may seem surprising at first that @ has the value 0 °6 for all non-zero values 
of v/I (the case v/l<1 and x/l<1 is trivial). However, it will be shown that, 
when the radiation emerges at an angle to the normal, the relative contrast 
decreases as v/l decreases. 


O14 (ii) 


@=15° 
(i) 
§=0c° 
Pz (iii) 
6 =30° 
(i) 
i @=45° 
+10 
(v) 
@ =60° 
0:08 
° m2 7 37 |2 27 


Ix 
Fig. 1—The directional intensity at the surface, I(x,0,6,0), 
against la for p=0 and (i) 0=0, (ii) 0=15°, (iii) 9=30°, 
(iv) 8=45°, (v) 6=60°. In all cases x,/J=10-5, v//=1, and 
a=0»1. 


For all cases involving non-zero values of 0 and 9, I(#,0,0,9) must be 
evaluated numerically for various values of the coordinates 7, 0, and 9, and of 
the parameters of the medium x,/I, v/l, and «. In Figures 1, 2, and 3, the results 
are shown for a typical case in which x,//=10-5, v/l=1, and «=0-1. The 
directional intensity I(x,0,0,¢) is plotted against lz for various values of 6 and 9. 
The curves obtained are obviously periodic, with wavelength 2z/l. The 
characteristic features, the average intensity I), the relative contrast. 8, and the 
difference in phase between J(0) and Z(0), &, are recorded in Table 1 for the curves 
of Figure 1. Also tabulated is the ‘‘ darkening ’’, Z,(0)/I,(0), which is compared 
with Eddington’s approximate relation I,(0)/I,(0)=(2+3 cos 0)/5. 
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The contrast decreases with 6, as expected. The phase difference between 
I (8) and I(0) arises from the location of the source below the surface ; it may 
easily be’shown to be consistent with a simple model in which the sources of radia- 


0-14 
(i) 
6 =15° 
(i) 
§8=0° 
(iii) 
0-12 6=30° 
I (iv) 
6=45° 
0-10 
(v) 
6 =60° 
0-08 
fo) m/2 7 37/2 an 


Ix 
Fig. 2.—As for Figure 1, but 9=45°. 


tion are distributed over a plane parallel to the surface at about unit optical depth, 
varying sinusoidally in one direction. As g—90°, —-0, as expected. The 
value of x,/l chosen in this case satisfies the condition for a continuously tapering 


TABLE 1 
FEATURES OF FIGURE | 


) 0° 15° 30° 45° 60° 
I,(0) 0-1333 0-1305 0-1222 0-1095  0-0930 
8 0-615 0-590 0-442 0-292 0-124 
E/2n 0 0-45 1-01 1-71 2-50 
,(0)/I(0) 1 0-979 0-921 0-821 0-697 
(2+3 cos 6)/5 1 0-980 0-920 0-824 0-700 


medium. As expected, almost identical curves are found for x)/J=10~* and 
1O=%; 

To facilitate the extrapolation of the curves of Figures 1, 2, and 3 to other 
values of v/l, many other cases have been evaluated. In Figure 4, £8 is plotted 
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against y, where y=log,,) v/l, for x/1=10-°, a=0-1, and various qalussh os 0. 
For all but large values of 0, the relative contrast 6 is approximately 0-6 when 


0-14 
(i) 
6=0° 
(ii) 
6 =15° 
0-12 
(iii) 
8 =30° 
| 
(iv) 
6=45° 
0+10 
(v) 
6=60° 
0-08 
° m2 7 anf2 A 


Fig. 3.—As for Figure 1, but 9=90°. 
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Fig. 4.—The graphs of 8 against ~, where Y=log,, v/l for p=0° and (i) 8=8°, 

(ii) 0=15°, (iii) 0=30°, (iv) 92=60°. In these cases, ~=0-1 and x,//=10-5. 

The points plotted thus, ©, are values of 6 calculated for «=0-2, 0=15°, 
and x,/J=10-5, 


v/l is large (i.e. small scale heights and coarse structures), but decreases to zero as 
v/l decreases; the greater the value of 0, the greater the rate of decrease. Points 
have also been plotted for 96=15° and «=0:2. The values of 6 undergo little 
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change provided $> 0-2. Thus the amplitude of the variations in intensity is 
directly proportional to the fractional variation in attenuation coefficient, «, 
provided « is small. 


VII. ConcLusions 

The emission of radiation from the surface (z=0) of a semi-infinite medium, 
below which the attenuation coefficient takes the form ( 3), has been examined 
in the conservative case (radiative equilibrium) assuming constant flux across 
planes parallel to the surface. 

The total intensity is given by (12) and the variation in total intensity across 
the surface by (13). These expressions reduce to the form given by Giovanelli 
(1959) in the special case of a medium uniform with depth. 

The behaviour of the directional intensity across the surface may be 
described by the relative contrast 8 (@=(1/a)(AI/I,)). The extent to which 
the increase in attenuation coefficient with depth affects 8 may be summarized 
thus : 

(i) For radiation normal to the surface, 0-6 <@ <1, depending on the ratios 
%/! and v/l. In the important case of a rmedium which tapers off continuously 
(%/1<1), 6 =0-6 for all non-zero values of v/l. 

(ii) The behaviour of radiation emerging obliquely is shown for a special 
case in Figures 1-3 and Table 1. Figure 4 shows that for a medium tapering off 
continuously, 0<%<0-6. For sufficiently coarse structures and small scale 
heights (v/l>1), 6 =0-6 for all but large angles of emergence, while in all cases, 
the relative contrast falls off as the angle of emergence increases. 
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MAGNETIC BAYS AT MACQUARIE ISLAND 
By C. 8S. ROBERTSON* 
[Manuscript received April 4, 1960] 


Summary 


Polar magnetic bays are the most striking feature of magnetic records obtained 
by the Australian National Antarctic Research Expeditions at Macquarie Island. 
Negative bays are more numerous than positive bays and have greater amplitudes 
and durations. As a rule, negative bays commence about magnetic midnight whereas 
positive bays commence about 5hr earlier. A daily reversal in direction of the bay- 
producing currents is indicated and the time at which this reversal takes place appears 
to vary with the seasons. During magnetically disturbed periods the ratio H bay 
amplitude to Z bay amplitude increases, indicating a northward movement of the 
bay-producing currents. At the same time auroras seen from Macquarie Island lie 
further north in the sky than usual. 


I. INTRODUCTION 

Since 1948 the Australian National Antarctic Research Expeditions have 
maintained a scientific station at Macquarie Island. A geomagnetic observatory, 
manned by officers of the Bureau of Mineral Resources, Geology and Geophysics, 
has operated there since 1950. 

Macquarie Island (54-5°S., 159° E.) lies about 1° north of the southern 
auroral zone (cf. Bond and Jacka 1960) and the geomagnetic records are 
characteristic of that region. ‘‘ Bay ’’-like indentations are prominent features 
on the magnetograms. In Figure 1 the record of October 24, 1954 is typical of 
‘“‘ disturbed ’? periods while that of October 30, 1954 shows a less disturbed 
record. In the early part of this record a bay (or pair of bays) is formed by 
increase in the horizontal (H) and vertical (Z) components. Later in the day 
another bay is formed by decrease in H and Z. In the following discussion the 


first type (increase) will be termed a positive bay and the second type (decrease) 
a negative bay. 


II. CHARACTERISTICS OF BAYS 

Most of the data used in this paper are from records of the sunspot-minimum 
year 1954. During that year 435 negative bays and 172 positive bays were 
recorded ; their distribution throughout the year is illustrated in the histograms 
of Figure 2. These show that most bays are recorded at or near the equinoxes— 
this is of course true of magnetic disturbances in general. Negative bays are 
more numerous and usually have amplitudes two or three times as great as 
those of positive bays; in 1954 the average amplitude of negative bays was 
about 285 gammas. Although the beginnings and endings of many bays are 
poorly defined, it is possible to say that the average duration of negative bays 
(approximately 100 min) is roughly twice that of positive bays. 


* Bureau of Mineral Resources, Geology and Geophysics, Melbourne. 
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Fig. 1—Tracings of Macquarie Island magnetograms. Scales of HD; 
and Z variations are given on the right. The vertical line indicates the 
time (13°4 hr G.M.T.) of local midnight. For comment refer to text. 
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The amplitudes of all 1954 bays were plotted against their durations, to 
find whether these two quantities were related. Apart from the fact that most 
positive bays have amplitudes less than 80 gammas and durations less than 
50 min, no trends could be discerned. 

Figure 3 shows the variation through the year of monthly mean and extreme 
range of commencement times of bays recorded in 1954 ; there is an obvious 
grouping of the commencement times of negative bays around ‘‘ magnetic 
midnight ” (as defined by McNish (1936)). Almost all positive bays occur before 
magnetic midnight and most of them occur several hours before magnetic 
midnight. Although these diagrams suggest that there is considerable over- 
lapping between the commencement times of negative and positive bays, it is. 


NUMBER OF BAYS RECORDED MONTHLY 


JAN, FEB. MAR. APR. MAY JUNE JULY AUG, SEPT. OCT. NOV. DEC. 


Fig. 2._Variation of frequency of occurrence of bays in the H 
magnetograms through 1954 at Macquarie Island. The dashed 
line refers to negative bays, the full line to positive bays. 


usually found that on individual days there is no overlapping ; all bays before a 
certain time are positive and all bays after that time are negative. For example, 
on the magnetogram of February 27, 1954 (Fig. 1) there is a sudden reversal of 
the direction of disturbances from positive to negative at the time marked 
“7”, The transition is not always sudden but it is possible on about a quarter 
of all days to pick such a transition time by inspection of the records. 


Several exceptions to this rule were recorded in the sunspot maximum 
period 1957-58 ; during large positive disturbances on stormy days a pronounced 
negative bay lasting about 1 hr sometimes occurred at about 7 or 8 hr G.M.T., 
some time before the general change-over to negative bays. Only one such 
case was observed during 1954. 

Figure 4 illustrates the pattern of variation of transition times throughout 
the sunspot minimum year 1954 and the sunspot maximum year 1957-58. There 
appears to be a marked seasonal effect ; the transition from positive to negative 
bays is a few hours earlier at the equinoxes than at the solstices. In the data 
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of 1954 the individual points cluster closely about the monthly mean ; in 1957-58 
the scatter is much greater but the seasonal variation of monthly means is 
comparable. 


24 
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Fig. 3.—Variation of commencement times of bays in the H magneto- 

grams through 1954 at Macquarie Island. Circles denote monthly 

means for negative bays, crosses for positive bays; the vertical lines 

indicate the extreme range of commencement times. The horizontal 
line is at “‘ magnetic midnight ”’. 
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Fig. 4.—Variation through the years 1954 and 1957-58 of 
monthly mean times of transition from positive to negative bay 
activity in the H component at Macquarie Island. 


It was thought that the transition times might be related to the degree of 
magnetic storminess and that the apparent seasonal effect might result from 
the larger number of disturbed days near the equinoxes. This could be so if the 
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transitions occurred earlier on disturbed days than on quiet days. To test this 
possibility transition times were plotted against the daily K-index sums; this 
was done separately for equinoctial and solstice months. The plot for solstice 
(or other than equinoctial) months in 1954 shows some tendency for the points 
to form an elongated group, suggesting a correspondence between K sum and 
transition times. However, this tendency is not repeated in the plot for the 
equinoctial months of 1954, nor does it appear in either of the plots for the period 
July 1957 to June 1958. It seems probable then that there is a real seasonal 
variation in transition times. 

So far this part of the paper has been concerned with magnetic bays which 
are evident on the H and Z component traces. Such bays are invariably 
accompanied by declination (D) component disturbances, but these are generally 
oscillatory in character and not of simple bay shape. The record of July 17, 
1954 (Fig. 1) shows a different type of bay. It is evident that there is little 
disturbance on either H or Z but there is a pronounced negative bay in the 
declination component D. In 1954 about 80 such bays were recorded. All of 
these were negative (i.e. the westerly declination increased) and all occurred 
in the intervening period between the times when positive H, Z bays and negative 
H,Z bays were most common. 

From the foregoing facts about magnetic bays at Macquarie Island it is 
possible to draw a number of fairly definite conclusions. It is clear from the 
reversal in bay sign that the bay-producing currents flowing near the southern 
auroral zone change direction from eastward to westward in the early evening of 
each day and that about this time currents often flow northward over the 
recording station producing bays in declination only. 

These facts may be explained qualitatively by several different theories 
(e.g. Maeda 1957) which postulate that currents are due to an excess of ions of 
one sign carried along in some way by ionospheric winds. It appears that the 
bay-producing current system rotates so that it presents the same general pattern 
as viewed from the Sun. This could result from a wind pattern fixed relative 
to the Sun while bay-producing ionization takes place in bursts, giving rise to 
considerable variation in current intensity. The normal diurnal variation may 
be produced by less intense but more regular ionization operating in conjunction 
with the same wind system. The large excess of negative bays over positive 
bays indicates a preponderance of bay-producing ionization on the night side of 
the Earth. The apparent seasonal variations in the change-over times from 
positive to negative bays suggests that magnetic bay studies may make useful 
contributions to the understanding of upper atmospheric wind systems. 


III. CORRELATION BETWEEN BAYS AND AURORAS 
It has long been recognized that the polar bay type of magnetic disturbance 


is often accompanied by visible aurora. Some aspects of the relation between 
the two phenomena have been investigated. 


The magnetic records and records of auroral observations for the year 1954 
were studied closely to determine whether simultaneous onset of magnetic bays 
and the auroras generally occurred. Of 77 nights on which the sky was clear 
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and it was possible to decide whether the first bay activity coincided with the 
first auroral activity, good agreement was noted on 35 nights. On 16 more 
nights there was agreement in the absence of both magnetic and auroral activity ; 
as absence of activity is comparatively rare in each phenomenon the probability 
of chance coincidence is low. Of the remaining 26 nights 14 would have shown 
agreement but for the observation of auroral glow at a time when there was no 
magnetic bay activity. Glows are frequently the upper limits of very distant 
auroras so that the magnetic effects of currents associated with them may be 
very weak; these cases cannot be regarded as conclusive. 


Summarizing these findings : 


Bays and auroras commenced together 35 nights 


Neither was recorded all night 16 
Results inconclusive 14 
Only one recorded, or times disagree 12 

Total 77 


It is clear that in most cases magnetic bays and auroras do commence at about 
the same time. 

Once bays have begun, magnetic fluctuations are usually so frequent that 
further correlation with the auroras is more difficult. It seems likely that a 
burst of ionization impinging on the ionosphere in the auroral zone produces the 
visible aurora and at the same time creates the condition in the ionosphere 
necessary for a bay, which commences simultaneously with or soon after the 
visible aurora. Arcs and bands are the auroral forms most frequently associated 
with the commencement of bays. 

Observations at Macquarie Island also demonstrate that the bay-producing 
currents and the aurora are more or less coincident in space. On most days the 
Z component trace closely resembles the H component trace, as for example 
on October 31, 1954 (Fig. 1). Easterly or westerly currents several degrees 
polewards of the recording station would account for this. However, on many 
magnetically disturbed days the H and Z traces are observed to be quite 
dissimilar: a general increase in the value of the H component may be 
accompanied by a general decrease in the value of the Z component 
or vice versa. These effects would be caused by easterly or westerly 
currents situated equatorwards of the recording station. For example, on the 
magnetogram for October 24, 1954 (Fig. 1) a positive bay or series of bays between 
~ about 5 and 10-5 hr is accompanied by depressed values in Z. After 10-5 hr 
a series of pronounced negative H bays occurs and is accompanied by generally 
increased values of Z lasting until about 17-5 hr, after which H disturbances 
again tend to be positive and Z disturbances negative. 

It may be noted that, although the Z component trace shows an average 
increase between 10-5 and 17-5 hr, it exhibits small variations which correspond 
in form and sign to baysin H. This suggests the continued presence of a varying 
current located polewards from the recording station in addition to a current 


situated equatorwards from the station. 
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In 1954 there were 15 days on which the Z trace remained neutral or became 
generally positive during negative H bays. All except one of the 15 days were 
very disturbed magnetically and, during all on which the night sky was clear, 
auroras were observed overhead or to the north of the island, rather than only 
to the south as was normal. Moreover, these 15 days include the times when all 
the most brilliant auroral displays of the year were seen. 

These observations at Macquarie Island strongly support the view that 
magnetic bay-producing currents and the aurora are closely related in both 
time and space. 
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MOTION OF THE AURORA AND MAGNETIC BAYS 
By F. R. Bonp* 
[Manuscript received April 4, 1960] 


Summary 


Observations from the ANARE station at Macquarie Island show that the develop- 
ment of an auroral display is characterized by a slow northward drift, east to west 
longitudinal motion, and an associated positive bay in the horizontal component 
magnetogram in the evening hours. This is followed by a sudden change in structure 
of the aurora and appearance of a negative bay. The negative bay persists for some 
hours, accompanied by west to east longitudinal motion and slow southward drift of the 
aurora. 

The co-latitude of the northern limits of the aurora is shown to be strongly dependent 
on level of geomagnetic disturbance. 

It is suggested that the observed motions are the mass motions of electrons which 
constitute the magnetic bay producing current within the boundaries of the auroral 
form. 


I. INTRODUCTION 

The data used in this paper are from the records of visual observations of 
the aurora at the Australian National Antarctic Research Expeditions’ (ANARE) 
station at Macquarie Island and from volunteer observers’ reports in eastern 
Australia collected by the Auroral Data Centre. 

The auroras considered were located on or to the north of the southern 
auroral zone which lies approximately 1° south of Macquarie Island (Bond and 
Jacka 1960). In this region auroras are usually aligned or extended roughly 
along parallels of geomagnetic latitude. The observations were recorded mainly 
during the period of the International Geophysical Year 1957-58. 

Attention is confined to the large-scale motions of the aurora as a whole 
across the sky and associated geomagnetic disturbance. Changes of internal 
structure such as movements of rays along an “arc” or ‘‘ band” are not 
considered. 


II. MAGNETIC BAYS 

A study of magnetic bays recorded at Macquarie Island during 1954 has 
been reported by Robertson (1960). Examination of the magnetic records for 
1958 showed trends similar to the records for 1954. 

Positive bays in the horizontal component H occurred less frequently 
than negative bays, their duration was shorter and their amplitude less. Many 
positive bays commenced during daylight, mainly between 14” 30™ and 18" 305 
local time ; while negative bays occurred during the night hours, commencing 
mainly between 182 30™ and 04" 30™. The time of change-over from positive 
bay to negative bay most frequently occurred between 21 00™ and 225 00m, 


* Antarctic Division, Department of External Affairs, Melbourne. 
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III. Motion OF THE AURORA 
(a) Typical Development of a Display 
Climatic conditions at Macquarie Island are extremely unfavourable for 
auroral observing ; the sky is very rarely clear for more than a few hours at a 
time. However, an examination of the Macquarie Island observations over the 
period 1951-1959 indicates that the main pattern of auroral behaviour may be 
described in terms of two phases. 


The first phase commences as a succession of one or more ares or bands 
appearing above the southern horizon. Aligned roughly along parallels of 
geomagnetic latitude, these arcs or bands may show rayed structure of intensity 
1 or 2 and very rarely intensity 3 (on the usual 0-4 scale). There is, then, a 
general northwards progression of the display. This may be a continuous drift, 
sometimes with a succession of advances or retreats ; or it may be a discontinuous 
progression involving fading of the aurora at one position and a new appearance 
at another. A new arc or band may form in advance of an existing display. 
An arc or band may disappear only to be replaced by another further to the south 
and the advance be repeated. When this phase has a duration exceeding 30 min 
it is almost invariably associated with a positive bay in the H magnetogram. 


In general, the northward drift continues until the beginning of the second 
phase which is characterized by a quite sudden increase in intensity and develop- 
ment of bright rayed structure. This may be followed by flaming. The discrete 
structure is replaced by moving diffuse rayed forms and pulsating and diffuse 
surfaces which latter slowly drift or progress to the south. This second phase is 
always associated with a negative bay in the H magnetogram. 


This main sequence may be halted at any stage and one or more stages may 
be absent. Occasionally after the second phase has begun ares reappear for a 
short time interval and the whole or part of a sequence may be repeated one or 
more times. 

The lateral drifts of individual forms and lateral progression of the display 
as a whole referred to above, take place typically at a rate of a few degrees of 
latitude per hour (cf. Jacka 1953) corresponding to speeds of the order of 
100 m/sec. Canadian observations by Kim and Currie (1958) on individual 
auroral forms show a modal value of north-south drift speed of 50-100 m/sec. 
They also find some evidence for an increase of drift speed with increase in 
geomagnetic disturbance. 


During times of very marked geomagnetic disturbance the complicated 
auroral situation over Macquarie Island is difficult to interpret in detail, but in 
these cases the pattern of the aurora as described above seems to apply to the 
observations from southern Tasmania and, for the major magnetic storms, to 
the observations from further north in Australia. 


The patterns of the auroral display as observed at Macquarie Island appear 
to be very similar to those described by Heppner (1954) and Elvey (1957 ) for 
the approximately magnetically conjugate location, College, Alaska Heppner 
notes that in general the first appearing arcs and bands approximately coincide 
in t'me with the recording of a positive bay in H, while the change-over to a 
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negative bay is associated with the appearance of bright rayed structure and 
Subsequent break up in the are nearest the equator. 

The Macquarie Island 1958 aurora records were examined in detail, in 
conjunction with the magnetograms to see to what extent a detailed corres- 
pondence, similar to that recorded by Heppner, existed. 

Auroral reports for 62 nights were examined. It was found that on 11 nights 
auroral behaviour was quiet and positioned well to the south of Macquarie Island ; 
these were all nights of little magnetic activity. On 46 nights, auroral behaviour 
and magnetic behaviour did not contradict the Heppner patterns, although in 
some cases the auroral record was somewhat meagre. Of the remaining 5 nights, 
on 3 of them the auroral behaviour is difficult to classify, while on 2 nights bright 
rayed displays occurred at times when the magnetograms showed positive bays. 
On both occasions these bright rayed displays were at least 1° of latitude to the 
south, rather than overhead or to the north, which is more usual. 


(b) Co-latitude of Auroras and K, 

Using data obtained during the I.G.Y., synoptic maps have been prepared 
showing, at 15-min intervals, the distribution of auroras in the region of Macquarie 
Island and eastern Australia. Auroral positions were computed from recorded 
elevations and azimuths of points on the lower border, assumed at a constant 
height of 105 km. 

From the maps the co-latitude of the most northerly lower border crossing 
geographic longitude 152° E. was extracted for each quarter hour and tabulated 
against the prevailing K,-index. These data are represented in Figure 1, which 
shows clearly the equatorwards appearance of the northernmost auroral form 
with increasing magnetic disturbance. Over the range of K, values 1-7, the 
mean co-latitudes in Figure 1 are comparable with those reported for homo- 
geneous arcs at Macquarie Island by Jacka (1953). A similar relation with the K- 
index recorded at Macquarie Island is apparent in the radio-echo observations 
tabulated by Gadsden (1959a). 


(ce) Longitudinal Motion 

Observations of longitudinal motion of auroras were made by G. Cowling, 
a member of the ANARE at Macquarie Island in 1958. Mr. Cowling also carried 
out a preliminary analysis. 

Measurements of elevation and azimuth to points on the lower border were 
made using an open sights theodolite and successive positions of the aurora were 
plotted on grids showing reference lines for elevation and azimuth. When, due 
to rapid movement of the auroral form, measurements were not possible, the 
initial position of the aurora was estimated and the direction of movement noted. 
Longitudinal motion, that is, motion in the direction in which the aurora is 
aligned or extended, was recorded in 33 cases. 

Of these 33 cases, 6 were of east to west motion and 27 of west to east motion. 
The former occurred at or before 22” 36™, the latter between 22h 01™ and 065 01m 
local time. In 5 cases only was it possible to estimate speeds of longitudinal 
motion ; these estimates lay in the range 200-1000 m/sec. Comparable speeds 
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are reported by Kim and Currie (1958), who also find an increase in speed with 
increase in magnetic disturbance. 

Direction of longitudinal motion of auroral fo ms was compared at corres- 
ponding times with the H component magnetogram. In general it was found 
that positive departures (+AH) from the quiet-day curve corresponded with 
east to west motions of the aurora, while negative departures (—AH) corresponded 
with west to east motions. This trend is illustrated in Table 1. 


GEOMAGNETIC CO-LATITUDE 


Kp INDEX 


Fig. 1.—Co-latitude of the northernmost auroral lower border 

versus K,. The circle represents the mean, the line represents 

the total range, of co-latitude, and the numbers give the total of 

quarter-hourly observations, for each value of K p (Data from 
gg. longitude 152 °E.) 

It should be noted that in the are or band type of aurora present during a 
positive bay, longitudinal motion is virtually along itself and can be recognized 
only when a discontinuity can be discerned. Discontinuities of recognizable 
size or shape are not of frequent occurrence. Moreover, it has already been 
noted that a proportion of positive bays commence during daylight, and that 
positive bays are of shorter average duration than negative bays. It is pre- 
sumably for these reasons that the number of recorded cases of east to west 
motion is small. 

These limitations are not present when radio-echo methods are used to trace 
the motion of ionization at auroral heights. Unwin (1959) has reported such 
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motions over Macquarie Island and has shown that the echoes come from a layer, 
a few to 25 km thick, centred between a height of 110 km and a little more than 
120 km. Unwin found that east to west moving echoes occurred between about 
155 30™ and 015 30™ local time, while west to east moving echoes occurred 
between about 185 30™ and 09 30™ local time. On any one occasion the change 
from westerly to easterly motion was quite distinct. Velocities ranged from 
about 200 to 4000 m/sec. Unwin reports a distinct increase in velocity with 
increasing values of the Macquarie Island K-index. 

A summary of some observations of auroral motion and motion of auroral 


ionization in the northern hemisphere is given in Table 2. The marked similarity 
with the data of this paper and with those of Unwin is clear. 


TABLE 1 
MOTION OF THE AURORA, MACQUARIE ISLAND 1958 


Number of Cases 
Direction of 


Auroral Motion N 
fe) 
+AH —AH Departure 
East to west a ae 4 1* 1 
West to east oe me os 27 — 


* Value —l5y, a small value for this station. 


IV. DISCUSSION 

The results of this and quoted works may be summarized as follows. North- 
ward lateral drift or progression and westward longitudinal motion of auroral 
luminosity and auroral ionization are each associated with positive magnetic 
bays. Southward drift or progression and eastward longitudinal motion are 
associated with negative bays. Whether there is a direct association of lateral 
drift and longitudinal motion of individual auroral forms cannot be ascertained 
from the present data. These are not precise enough to enable detection of 
slow lateral drift during the relatively short time intervals in which longitudinal 
motion is usually discernible. Speeds of longitudinal motion are an order of 
magnitude greater than lateral drift speeds. 

Positive magnetic bays and associated auroral movements northwards and 
westwards are characteristic of the evening hours. Negative bays and associated 
movements are characteristic of the night and early morning hours. The 
transition is usually rapid and may occur within a time range of several hours 
(but usually within a range of about 1 hr) centred a little before midnight. 

The maximum co-latitude and the speed of longitudinal motion and possibly 
of lateral drift increase with geomagnetic disturbance. . 

It is apparent that the longitudinal motion takes place in the direction of an 
electron current which would produce the observed sign of magnetic bays. 
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Gadsden (1959b) has shown that the occurrence of auroral ionization (giving 
rise to radar echoes) and of visible auroral luminosity are statistically independent 
events. However, the morphologies of these phenomena are similar. 


The simplest hypothesis consistent with the observations is that the radio 
and optical methods detect different attributes (produced by different causes, 
either of which may be absent) of the one basic auroral phenomenon. This is 
characterized by a slow lateral drift and more rapid longitudinal motion. This 
motion of patterns of ionization and/or luminosity is identified with the mass 
motion of electrons which constitutes a magnetic bay producing electric current 
within the narrow confines of the auroral phenomenon. At any one time several 
such phenomena, generally roughly parallel to one another, may be observable 
from the one station. 


V. ACKNOWLEDGMENTS 
The writer acknowledges the benefit of discussions with Dr. F. Jacka, Chief 
Physicist, and Mr. K. D. Cole of the Antarctic Division, Department of External 
Affairs, and also with members of ANARE. The Macquarie Island magneto- 
grams were made available by the Bureau of Mineral Resources, Geology and 
Geophysics, Department of National Development. 


VI. REFERENCES 

Brzss, R. C., Gartrern, C. W., and Kimpaty, D. 8. (1955).— Astrophys. J. 122: 205. 

Bonn, F. R., and Jacxa, F. (1960).—Aust. J. Phys. 13: 610. 

ButtoucH, K., and Kaiser, T. R. (1955).— J. Atmos. Terr. Phys. 6: 198. 

Buxwiovues, K., Davipson, T. W., Katsmr, T. R., and Warxtns, C. (1957).— J. Atmos. Terr. Phys. 
Ti, BY ie 

Exvey, C. T. (1957).—Proc. Nat. Acad. Sci., Wash. 43: 63. 

GapspENn, M. (1959a).—Ann. Géophys. 15: 395. 

GapspeEn, M. (1959b).—Ann. Géophys. 15: 403. 

Heppner, J. P. (1954).—Thesis, California Institute of Technology. 

JacKA, F. (1953).—Aust. J. Phys. 6: 219. 

Kim, J. S., and Curriz, B. W. (1958).—Canad. J. Phys. 36: 160. 

Lyon, G. F., and Kavapas, A. (1958).—Canad. J. Phys. 36: 1661. 

Meex, J. H. (1954).—Astrophys. J. 120: 602. 

Metnet, A. B., and Scuurre, D. H. (1953).—Astrophys. J. 117: 454. 

Nicnots, B. (1957).—Sci. Rep. No. 1, Air Force Contract No. AF 19(604)—1859, Geophysical 
Institute, University of Alaska. 

Rosertson, OC. S. (1960).—Aust. J. Phys. 13: 470. 

Unwin, R. S. (1959).—Nature 183: 1044. 


A DYNAMO THEORY OF THE AURORA AND MAGNETIC 
DISTURBANCE* 


By K. D. CoLEt Boe 
[Manuscript received April 4, 1960] 


Summary 

A model of an aurora regarded as a plane slab of highly ionized air parallel to the 
geomagnetic field within the ionosphere is examined. The model is stable in the presence 
of a wind of neutral molecules which, blowing the slab across the geomagnetic field, 
generates an electric polarization field perpendicular to its faces and a current along 
its length. This current is concentrated in a small height range and is chiefly due 
to electron drift. 

The aurora moves normal to itself with a speed of the same order as the wind speed, 
while the drift carries its luminosity and ionization pattern along its length at a speed 
an order of magnitude greater. Measurement of these speeds will permit determinations 
of collision frequencies of ions and electrons within the aurora. 

Assuming, in the equatorial vicinity of the auroral zone, an equatorwards wind in the 
evening and a polewards wind in the night and morning hours it is possible to explain the 
major movements of the aurora and associated bay type magnetic disturbances. It is 
suggested that the magnetic disturbance indices K and K, are indicators of wind speed. 
A likely mechanism of maintenance of luminosity and ionization in the aurora is outlined. 
The general features of magnetic disturbance current systems within the auroral zone are 
briefly considered in relation to the theory and it is suggested that auroras are visible 
manifestations of the current flow. 

Qualitative speculation on a possible source of energy for winds in the polar dynamo 
region leads to a unified account of the Van Allen radiation belts, auroras, and some 
(low latitude) airglow. 


I. INTRODUCTION 

Two parameters are required to define the atmospheric dynamo problem, 
namely, the conductivity of the ionosphere and the velocity of the wind. The 
quiet daily variations in the geomagnetic field can be explained in terms of the 
conductivity of the quiet ionosphere in conjunction with the velocity of wind 
inferred from atmospheric tidal oscillations (Baker and Martyn 1953). To 
define a model for the explanation of geomagnetic disturbance one has, in the 
absence of other criteria, the choice of specifying a disturbed conductivity 
pattern (Obayashi and Jacobs 1957) and/or a disturbed wind velocity (Maeda 
1957). 

The present paper does not attempt to explain the world-wide distribution 
of magnetic disturbance, but rather considers a highly probable local situation. 
A model of the aurora is defined and its properties examined. These are compared 
with observational data, especially on large-scale movements of the aurora and 
associated magnetic disturbance. 


* Preliminary accounts of this work were presented at ANZAAS, Perth, August 1959, and at 
the Antarctic Symposium, Buenos Aires, November 1959. 
fj Antarctic Division, Department of External Affairs, Melbourne. 
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Although a one-to-one correspondence has not been established between the 
luminosity, magnetic, and radar observations of the aurora, these are here treated 
as aspects of the one phenomenon (Bond 1960) supposed as having the form of the 
common auroral “ arcs’’ and “ bands’’. These are typically long (thousands 
of kilometres) compared to their vertical extent (tens to hundreds of kilometres) 
and to their thicknesses (a few kilometres). 


Il. THEorY 
The current density at any point in an ionized gas is given by 


j=o,E,,+0,E, +o,(H xE,)/H, Rs Kos: op douse ation oneal (1) 


where E is the electric field and || and | indicate components parallel and per- 
pendicular to the magnetic field H. o 9, c,, co, are the direct (parallel to H), 
Pedersen (transverse), and Hall conductivities respectively. 

In a lightly ionized gas 


oo —Ne?(1/m,y; +1/m,»,), 
Oy =Ne(t, +t,), 
g=—WNer(h;—h,), 


where N =electron density, e=charge on electron, 


tei =Ve,il Me,i( VE; EF coe )s 


hej — We,i/Me,i( V2; = wo?) 3 


where v,,;, @,,;, and m,,; are the collision frequencies (with neutral particles), 
gyro frequencies, and masses of electrons and ions respectively. The assumption 
is made that the mobilities t and h are independent of ionization density NV. 


In general o, is much larger than oc, or o,; however, our first considerations 
are restricted to situations in which E,, is zero. Later (in Section IIT (d)) limited 
consideration is given to the effects of a non-zero Ej). 


(a) The Model 
(i) The Steady State—Clemmow, Jobnson, and Weekes (1955) have 
examined the motion of a cylindrical irregularity in an ionized medium. We 
investigate another elementary model now proposed. As with Clemmow, 
Johnson, and Weekes, dissipative effects such as diffusion and recombination 
are neglected here. 

Given a homogeneous, unbounded, weakly ionized neutral gas with ionization 
density N (ions) in the presence of uniform and steady orthogonal electric and 
magnetic fields E°, H, is there a steady state solution for which an infinite plane 
slab parallel to H and of ionization density N’ (ions) persists in the gas without 
change of form ? 

The unit normal (n) outward from the slab makes an angle « with E°. 
Rectangular axes 0(x, y, 2) with unit vectors e,, C2, 3, have the directions of 
E°xH, E°, and H respectively (see Fig. 1). The superscript 0 is used for 
quantities in the undisturbed state in the absence of the slab and the superscript p 
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for the incremental quantities resulting from polarization charge. The total 
electric field is EX +E”. The geometry of the model will make E?” orthogonal to H. 
Drift velocities of electrons (negative ions are neglected) and positive ions 
(assumed singly charged) are given, respectively, by 


v=—e(tEs +5 5) idiva oles sie Rie see (2a) 
HxE 
vine HE i wes =) rrr Se (2b) 


A prime is used to indicate quantities inside the slab. Then, if V denotes 
the velocity of the faces of the slab, the condition that the electronic and positive 


E°xH 


Fig. 1— Gas between the parallel planes differs in electron density from that 

of its surroundings. Magnetic field H is parallel to the planes and orthogonal 

to E® (an applied electric field). e,, €,, €; are unit vectors in the direction 

of E°x H, E°, and H, and n is an outward unit normal to the parallel 
planes, 


ion currents should separately be continuous yields two relations which must be 
satisfied on the surface of the slab, namely, 


n.(vg-tv3—V)N =nlve-y. VN, eee (3a) 
fh. (¥e--vi SVN =n AVeeve VN ee (3b) 
From equations (2) and (3) it follows that 


[w,/v,t+-o,/v,]eH® sin « 
n = a Mi Wien py WT na es e. 0 SesOiie eee We eis. 7 (4) 
{(ve+@e)/v.}m.+{(vi +3)/v;}m; 
It is thus shown that a steady state is possible in which the slab does not 
change form and that this state is independent of electron density either inside 
or outside the slab. 
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If E® is due to a wind of velocity —Ue, and 8 is the angle between — Ue; 
and n, then sin $=—cos « Then, in a frame of reference moving with the wind, 


site —(o,/v,+a,]v;) U cos § 
ie 
®,/V.+v,/@,+@;/V;+y;/0,; 


Now in the ionosphere w,/v,>«,/v;, hence* 


V.ties—U, Cos 8/(14-v,v)/0,0,). 0... 05.255. (6) 


In a frame of reference in which the wind velocity is U (i.e. a rest frame on the 
ground) the velocity of the faces is given by 


Vereen UN Wee ee Abe eae (7) 


The ratio V,/U.n is plotted in Figure 2 as a function of «,/v,, assuming with Baker 
and Martyn (1953) that (o,/v,).(v;/@;)=1650. 
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Fig. 2.—Curve (a), (63/6,)/(62/01)max, 28 a function of q,/v,. 
Note: (63/0;)max.%20. Curve (6), V,/U.n as a function of 
w,/v,- Height scale is approximate. 


(ii) Current Flow.—We now investigate the current flowing in the slab. 
The polarization field will be entirely internal and normal to the slab faces. The 
rate of accumulation of electrons on a face of the slab is 


an,/at=eN'n.(v9+v% —V)—eWn.(ve—V), 
whilst the rate of accumulation of positive charge on a face is 
On, | Ot —eN'n.(v9-+v¥ —V) —eNn.(vi—V). 
Thus the rate of charge (q) accumulation is 


0q/at —eN'n.(v$—ve sy? v3) —eNn.(v3—v9). 


* After completion of this work a recent paper by Clemmow and J obnson (1959) came to 
hand. Equation (6) is almost identical with their equation (61) for the velocity of propagation 


of an irregularity. 
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Now the polarization field E” is —4nqn, hence 


(—1]47) dH?’ /ot =(N’ —N )e2(t; +t,)E° cos a+ N'e?(t; +t,) EY 
+(N’—WN)e?(h; —h,)H° sin a. 
Thus, 


Ry hae (cos oa nee sin “| [1—exp { —47xWe?(t; +7,)¢}]. 


Let us resolve the total current density j into two components, j, normal to 
the slab and j, parallel to the slab in the direction nxH. It follows, since 
j=N.(v; —=V,), that 


jn=0,E° cos «—o,H® sin ao 


oO" x —N ipo (0s ao — ©? sin 2) exp (—47ojt), .... (8a) 
1 WN ony 
Gt 
jp = —03E° sin a—o,H° cos a+ ie sin % 
1 
N’—N Gye ; 

0, f° (cos a= a sin 2) exp (—4703t), os ASD) 
4 pCa? COS Ge BIN, aves! e ols arok ~ See ee eee (8c) 
jp=—6,H sin a—GgH® COB, «eee eee ee ee eee (8d) 


where 
63=6,+63/6). 
Thus the exponential time constant for approaching the steady state is 


1/426}. If E° arises from a wind with speed U of neutral particles in the direction 
opposite to e,, the internal current in the steady state is given by 


dn==0, UH sin 6--6,U0H GOs 3.) os 5.c5 eee ness (9a) 
o2 
jp=03UH cos 8—o,UH sin 8— ~UH £08. 8 cess (9b) 
1 
and the polarization field by 
Ev Bes Dinsial CH (sin 5+ cos s) mtn oie (10) 
N ony 


This polarization field has a maximum (in the case of N-+0) of 
| He’ nant = UNE A/G) shiSifa] (oleh ie) 0:4 Ve) 0) SMiarens. (11) 


when $=4n—tan-! o,/o,. Clearly | H® |max. may be much greater than UH, 
the maximum electrostatic field derivable from the movement of an isotropic 
conductor across a magnetic field (Alfvén 1950). 


(b) Acceleration of Charged Particles 
Under the influence of an electric field E and magnetic field H, charged 
particles move in trochoidal paths between collisions. The energies of the 
particles continually oscillate by virtue of movement to and fro along E. It 
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may be shown from the equations of Chapman and Cowling (1953, p. 324) that 
in the absence of collisions the distance d travelled to and fro by a particle of 
mass m and charge e is 2(v/@+eH/mw?) where w is its gyro frequency and v its 
speed prior to the application of H. In the field | H?’ |max, given by equation (11), 


a=2(7 Saal tae eke eae (12) 

@  M*,/ oy 

Thus the maximum change in energy of the particle over a half cycle is 
AE=2| +mvU 4/(o3/6,) +mU%(o,/6,)|, ........ (13) 


the sign + or — being used accordingly as the charge of the particle is positive 
or negative. Except for gas particles of initially very high velocities and for 
even moderate values of U, equation (13) becomes 


Nereis iae Maes oO ee Pla See Oe eae (14) 


In an isotropic gas o3;/o,—1, of course. In a lightly ionized tenuous gas 
(e.g. ionospheric F region) o3/o,~1. In the ionospheric # region o,/o, may 
reach a value of about 400 (Baker and Martyn 1953). In stellar atmospheres 
values greater than 10! are common (Piddington 1954). 

In the dynamo region of the Earth’s atmosphere significant values of U may 
arise through the action of a neutral particle wind; the effects of such winds 
are considered in later sections. In stellar atmospheres U may attain significant 
values through the action of radiation pressure of the Milne kind (Milne 1926). 

It appears likely that the movement of the Earth’s exosphere through the 
interplanetary magnetic field or the solar system through a galactic magnetic 
field may produce significant electrostatic fields which would manifest themselves 
as anisotropies in cosmic ray intensities. Such fields may also accelerate 
individual charged dust particles. 


III. APPLICATION TO THE AURORA AND MAGNETIC DISTURBANCE 

Baker and Martyn (1953) have calculated the conductivities as functions of 
height in a model ionosphere. It is seen from their work that the bulk of atmos- 
pheric dynamo current flows between 90 and 120 km height. Also between the 
same height limits Hall currents far outweigh Pedersen currents. Thus a good 
approximation in the dynamo region is 

GeO O Geen tatters Set oe ee (15) 

Weekes (1956) exhibits graphs of the mobilities of electrons and ions 
separately as functions of height. From these graphs it is seen that in the 
dynamo region the Hall current due to electrons dominates that due to ions. 
Hence throughout the greater part of the dynamo region we may say 

INGA, 


25 ee, OF 
M, VetWe 


O 9 O (electrons) — 


Since w,/v,>1 in the region considered, 
mae 0) LIINE IM, a eate en bea ore (16) 
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Since electrons dominate the flow of current, 
JR SHNC i a ee ew ee (17) 


where v, represents the drift velocity of electrons. 

We now identify the model described in Section II (a) above with an aurora 
immersed in the ionosphere and identify U with the velocity. of horizontal wind 
perpendicular to the (near vertical) polar geomagnetic field. The current flow 
within an aurora is given by equations (9). From equations (9b), (15), (16), 
(17), we obtain for the drift velocity of the electrons within the auroral mode 


Used —(Gy] G4) U0. COS O sien dee eas (18) 


in the direction +nxH. (It is assumed that N’SWN so that the first term of 
equation (9b) is the dominant one.) 

Figure 2 shows the variation of o,/o, a8 a function of w,/v, and height. The 
ratio o,/o, reaches a maximum value of about 20 in the dynamo region (cf. 
Baker and Martyn 1953). The maximum value depends on the value chosen 
for the ratio of collision frequencies of ions and of electrons. This point will be 
discussed below (Section III (e)). 

We now identify the speed V, given by (7) with the horizontal movement 
of an aurora normal to its length and the velocity of drift of electrons v, with 
that along the aurora of large-scale features of the luminosity pattern and of the 
ionization pattern producing radio echoes. This is justified when N’SN, for 
HY’ (equation (10)) and hence v, (equation (2a)) and v, become independent of WV’. 
The theory predicts that for the lower portion of the ionosphere V, is of order U 
and v, of order 10U (see Fig. 2). 


(a) Auroral Movements and Magnetic Disturbance 

If we now suppose that patterns of ionization established within the auroral 
zone are blown by a neutral particle wind it is possible to account for a number of 
observations on auroral morphology. Appropriate wind speeds, generally 
10-100 m/see and sometimes more, are observed in the dynamo region (Elford 
1959a, 1959b). 

In the vicinity of and equatorwards of the auroral zone, auroras are oriented 
on the average roughly east-west. If then the wind in this region has a diurnal 
variation, being equatorwards in the evening and polewards later in the night 
and morning, one would expect a similar movement of auroras. This is observed 
(cf. Bond 1960). Similarly we explain the diurnal variation of latitude of ares 
reported by Jacka (1953). Equations (17) and (18) also explain Robertson’s 
(1960) observation on the occurrence of positive bays in the evening and negative 
bays later in the night and morning. 

The drifts of auroral ionization (Unwin 1959) and luminosity patterns 
(Bond 1960) along the aurora with speeds an order of magnitude greater than 
typical wind speeds are explained by equation (18). The diurnal variation in 
ionization drift speed (Kaiser 1958) is explained by the cos § term in (18). 

Heppner (1954) shows that the bulk of the magnetic bay producing current 
does in fact flow in or near auroras. This indicates that N<N’. To estimate 
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the magnitude of a magnetic bay we take account now of only the current flowing 
in visible auroras (equation (9b)). The currents flowing outside an aurora 
indicated by equations (8c) and (8d) may contribute significantly to disturbance 
if they occupy a great enough region. At Macquarie Island when the magnetic 
K-index is 6 there is a magnetic disturbance of order 500 gamma and auroras 
typically fill the whole sky. Thus the total effective (latitude) thickness of 
auroras is of order 100 km. Recent rocket measurements of electron densities 
in auroras give values of 10°10? per cm® (Jackson and Seddon 1959). If the 
auroral dynamo be 10 km thickness in height and at 100 km above the ground, 
then a disturbance at the ground of over 500 gamma is easily explicable. 

It has been observed (Zaborshchikov and Fediakina 1958) that during 
positive magnetic bays hin, is generally above 100 km, whilst during negative 
bayS hm, is generally below 100km. Moreover, increased absorption of 
extraterrestrial radio waves is more associated with negative than with positive 
bays. These observations can be qualitatively explained by a modification of 
the theory which will not impair the foregoing conclusions. 

We consider the effect of allowing for the inclination of the Earth’s magnetic 
field to the vertical. Figure 3 illustrates a cross section of an aurora in the 
magnetic meridional plane with the wind blowing equatorwards. Positive 
charges will gather on the leading face of the aurora and negative ones on the 
trailing face. Thus in this situation electrons will tend to be lifted out of the 
lower regions of the ionosphere so that their effect on absorption is lessened. 
Apparently they are lifted to a region where their horizontal drift velocity is less 
(namely, above the peak in drift velocity near 100 km, cf. Fig. 2). However, 
when the wind is polewards electrons will be depressed into lower regions where 
their effect on absorption is increased yet the bulk of them still occupy a region 
of high drift speed. This also explains Robertson’s (1960) observation that 
negative bays are of larger amplitude in general than positive bays. 

The two physically distinct components of velocity of an aurora are V, 
and v,. Since v, is of order 10 and more times V,, the arbitrary resolution of the 
movements of auroras north-south and east-west will not in general give informa- 
tion about V, and v, unless the orientation is known. However, as mentioned 
above, near the auroral zone the average orientation of auroras is east-west, so 
that on the average we may associate north-south movement with V, and 
east-west movement with v,. Now Bond (1960) exhibits a graph for the latitude 
of auroras as a function of K,, whilst Unwin (1959) exhibits a graph for the drift 
velocity of auroral ionization as a function of A at Macquarie Island. The 
average values from these two sets of data are plotted in Figure 4. (The fit is 
achieved by choosing the same ordinate for the drift 3000 m/sec and 41° gm. 
co-lat. at K,, K=9 and K,=0 corresponding to 100 m/sec and 27° gm co-lat.) 
In terms of the present theory, the correspondence revealed in Figure 4 suggests 
that K and K, are indicators of wind speed. 

Thus, though changes in ionization density are important, it appears that 
the wind in the polar dynamo region is a prime factor in determining the mor- 
phology of auroras and magnetic disturbance. A source of energy for this wind 


will be discussed later (Section IV). 
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(b) Luminosity of Auroras 
In the model presented in Section II the only electric field available for the 
acceleration of charged particles is that given by (10), which is orthogonal to H. 
The maximum energy acquired by a charged particle which does not suffer 
collision is given by equation (13). Thus a 100 m/sec wind may accelerate 


Ux H 


Fig. 3.—Effect of inclined magnetic field, 


charged nitrogen and oxygen molecules to energies of approximately 1 eV 
Under favourable circumstances this field may contribute to the luminosity of 
auroras. 

However, in a real aurora significant electric fields may develop parallel to 
the geomagnetic field. For conductivity is a function of height and therefore 
the build-up of polarization at the faces of the aurora will not be uniform at all 
heights. Under these conditions current must flow along H—this is equivalent 
to creating an electric field parallel to H. The effect of this current flow along H 
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will be to delay the steady state suggested by equation (9). An upper limit to 
this electric field parallel to H may be obtained by assuming that in the initial 
stages of luminosity development all the current which flows in the direction 
of the normal n (see Fig. 1) is redirected along the geomagnetic field lines. Thus 
from equations (1) and (8a) (putting t=0, H®=UH, §=i4n+«a and N’SN) 


E,, ~7r(o,UH sin 8+02UH cos $)/o0,  ........ (19) 


where r is the ratio of the height extent and the latitude extent of the dynamo 
section of an aurora. An auroral dynamo may be 20 km in height extent and 
less than 1km thick. Since in the dynamo region (c,/69)max, and (G5) o5) max. 
are about 10~? (Baker and Martyn 1953), it may be shown that ,, values of order 
10~° V/em are obtainable from a wind of 30 m/sec. Chamberlain (1955) in a 
study of auroral rays concludes that an electric field of order 10-5 V/em is 
adequate to explain their luminosity. 
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Fig. 4.—Superposition of two sets of data: (OQ) the EWE average 

drift speed of auroral ionization as a function of K-index at Macquarie 

Island (Unwin 1959); (x) the average maximum co-latitude of auroras 
as a function of K b (Bond 1960). 


The luminosity of an aurora appears to move with the same velocity at all 
heights. This observation does not conflict with the theory. The electric field 
causing luminosity at great heights moves with a speed dictated by the electron 
flow at the base of the aurora in the dynamo region. The auroral dynamo may 
occupy only a small height range and so only a small range of values of V, and 
v4 may be effective in determining the speed of luminosity. It appears that the 
values of V, and v, associated with the bulk of “ dynamo ” electrons will dominate 

and that values of V, and , associated with other electrons may contribute to 
diffuseness at the edges of auroral luminosity. 

As mentioned above, the leakage of charge along the field lines increases the 
time to reach the steady state (maximum) electron drift (v,;) in an aurora. tt is 
thus expected that the luminosity peak of an aurora should precede the associated 
magnetic disturbance at the ground, as has been reported by Bless et al. (1959). 

Cc 
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(c) The Collision Frequencies of Ions and Electrons with Neutral Particles 
From equation (18) it follows that o,/o,~v,/U cos 5. It can be seen from 
the equations of Baker and Martyn (1953) that when v,/v,210 (as is presumably 
the case in the ionosphere), then (o3/6,)max.4(v;,/¥,) (to within 1 per cent.). 
Thus 
vy 


Me V4 : 
Yio 7 a atu PE Dy Pee (20) 


V y 


s e 


(o,/c,) U cos So ahs 


Also it may be shown that 


Therefore, from equation (18) 
ViGwy joe «is. Aeeeees sees (21) 


Thus the collision frequencies of electrons and ions may be determined by ground- 
based observations of wind and auroral movements. 


(d) World-wide Magnetic Disturbance 

Dynamo theory of world-wide magnetic disturbance (Maeda 1957 ; Obayashi 
and Jacobs 1957) in its usual form argues from a diffuse continuously (spatially) 
varying current. It does not envisage structural discontinuities of the kind 
discussed here. It may be capable of explaining on a world scale average magnetic 
disturbance but would fail to draw out the small-scale and short-term detail. 

It is suggested that auroras are just visible lines of the magnetic disturbance 
current system. This is a natural outcome of the above theory and is supported 
by observations as follows. At Macquarie Island (just outside the auroral zone, 
Bond and Jacka (1960)) the orientation of auroras is generally east-west (Bond 
1960); at Mawson (just inside the zone) the aurora often shows large-scale 
curvature (as distinct from relatively small-scale convolutions and kinks which 
appear as a general feature of auroras in most places) (ANARE records) ; at 
Scott Base (well inside the zone) the aurora often takes the form of complete 
loops in the sky (Carter Observatory 1957, 1958); at Dumont d’Urville (well 
inside the zone) the orientation of arcs in 1954 generally lagged the orientation 
of the Sun-Earth line by an acute angle (Weill 1958). These observations are 
qualitatively what could be expected following inspection of diagrams for 
disturbance current systems (see Figs. 227-240 of Vestine et al. 1947). 

The polar disturbance current system (Vestine et al. 1947) could be explained 
in terms of a wind system with a source at the centre of the “ afternoon ”’ current 
cells and a sink at the centre of the ‘‘ morning ”’ current cells. The lines of flow 
of the wind would be roughly orthogonal to the current lines, hence across the 
polar cap wind in the dynamo region would blow from the afternoon sector to 
the morning sector. 


IV. ENERGY SOURCE FOR THE WIND SysTEM 
Corpuscular bombardment of the Earth may be effected by a solar wind 
(Parker 1959) or hydromagnetic shock waves (Singer 1957) or weak hydromagnetic 
waves (Cole 1959). Let us consider the last mentioned. Hydromagnetic waves 
can deliver energy in two modes; the extraordinary and the ordinary mode. 
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The extraordinary mode travels with the Alfvén speed V, in all directions, 
whilst the ordinary mode travels with the phase velocity V, cos ), where ) is 
the angle between the wave normal and the magnetic field (Piddington 1955). 

In the case of hydromagnetic waves coming all the way from the Sun or 
being generated at the interface of the Karth’s outer atmosphere and the inter- 
planetary medium, the energy of the ordinary mode will accumulate principally 
in high magnetic latitudes and in the lines of force connected to them, whilst the 
energy of the extraordinary mode may be expected to accumulate in low magnetic 
latitudes and in the lines of force connected to them because of the curvature of 
the geomagnetic field lines. This would explain in a simple unified fashion the 
occurrence of the two maxima in the radiation intensity at heights above the 
Earth’s surface (Van Allen 1959). This mechanism is without the difficulties 
of the actual transmission of particles. For, whereas disruption of the Earth’s 
magnetic field must be invoked to explain how low energy extraterrestrial particles 
become trapped in the Earth’s magnetic field, no such disruption is necessary for 
the capture of energy from hydromagnetic waves. The latter is merely a function 
of the geometry of the capturing system. 

Let us consider the effect of hydromagnetic waves of all wavelengths coming 
either from the Sun or from interplanetary space near the Earth. Waves of 
length less than the scales of inhomogeneities (e.g. in density or magnetic field) 
in the Earth’s exosphere would retain their character as waves and suffer damping 
and loss of energy, as demonstrated by Piddington (1959). However, waves of 
length greater than the scales of inhomogeneities would soon degenerate into a 
system of gyrating particles in the Earth’s field. These particles will lose energy 
chiefly at their mirror points by collision ; in this way ionization is established 
and energy is delivered to the dynamo region. Also the polar dynamo region 
may receive energy from low energy cosmic rays. 

One would expect an “ auroral zone ’’ associated with each of the naturally 
occurring radiation belts, just as a temporary “‘ auroral zone ’’ was associated with 
the man-made radiation belts of the Argus experiments (Christofilos 1959). 
This would account for the maximum of intensity of the airglow observed near 
em. lat. 35° by Nakamura (1958) as the ‘“ auroral zone ”’ of the inner Van Allen 
radiation belt. Also Sandford (1959) has observed that ‘‘ the intensity of the 
O, lines 5577 A and 6300 A in the airglow and aurora vary as similar functions 
of the K,-index and therefore suggests that the increased O; emissions in the 
airglow during periods of increased magnetic activity are essentially the same 
phenomena as aurora’”’. This supports the above thesis. 


V. DISCUSSION 
The above theory correlates many widely dispersed observations. The 
_ application of the dynamo idea to the theory of auroras is not new; Vestine 
(1954) quotes references to work up to that time. However, the present approach 


is new. 
The theory leads to the inference that the “ weather ” of the polar dynamo 


region is under the control of corpuscular or hydromagnetic wave bombardment ; 
these effects dominate the effects of tidal oscillation. 
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The following picture of an aurora has evolved from consideration of the 
above theory. Some agency, probably charged particles from the outer Van Allen 
radiation belt or hydromagnetic waves directly from interplanetary space, 
establishes a wind system in the polar dynamo region. The wind blows an 
established pattern of ionization across the magnetic field. This generates 
electric fields, in the region of high electric conductivity, which drive currents 
and cause luminosity and maintain ionization for very much longer times than 
would be permitted by ‘ quiet ’’ ionospheric conditions. 

To test the theory simultaneous observations of wind velocity (by the meteor 
trail drift method) and auroral orientation, movement, and height are required. 
These are required to check that indeed the speed of an elongated aurora normal 
to itself is just V,. If this were established, the aurora would become a very 
valuable visible adjunct to the study of winds in the polar dynamo region. 
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RAYLEIGH WAVE DISPERSION FOR A SINGLE LAYER ON AN 
ELASTIC HALF SPACE 


By B. A. Boir* and J. C. BUTCHER* 
[Manuscript received May 3. 1960] 


Summary 


Numerical solutions of the period equation for Rayleigh waves in a single surface 
layer were calculated using the SILLIAC computer at the University of Sydney. Values 
of the phase and group velocities for both the fundamental and first higher mode are 
tabulated against period for eleven models. These related models allow a sensitivity 
analysis of the effect of variation in the seismic parameters. 


I. INTRODUCTION 

The complicated algebraic form of the period equation for Rayleigh wave 
propagation in a semi-infinite layered solid elastic medium discourages manual 
computation and only a few solutions for fairly simple cases have been calculated 
in this way (Ewing, Jardetsky, and Press 1957). Recently, solutions for the 
fundamental mode in three cases, each with three layers, have been obtained on 
EDSAC (Stoneley 1954); and IBM 650 machines have been used to evaluate 
the fundamental, first, and second higher mode velocities for a 3-layered model, 
and two 4-layered models (Oliver, Dorman, and Sutton 1959). 

In the present paper solutions for the fundamental and first higher mode 
have been obtained for 11 related 2-layered models. Values of the elastic 
constants were chosen to fit the seismic velocities determined for the Western 
Australian crust (Bolt, Doyle, and Sutton 1958), so that the solutions are relevant 
to studies of the crustal structure of Australia. Computations were carried out 
on SILLIAC at the Basser Computing Department, the University of Sydney. 


Il. THE PERIOD EQUATION 
Consider plane waves of length 27/x, phase velocity c, propagated in a 
homogeneous perfectly elastic medium, 1, of uniform thickness H, welded to the 
plane boundary of a semi-infinite similar medium, 2. 


In such media, the velocities «, of compressional and shear waves are 


O,” ==(2 2.) 6,7 Bet pees eee ean eee (1) 


where e is the density and A, wu the Lamé parameters. 


* Department of Applied Mathematics, University of Sydney. 


RAYLEIGH WAVE DISPERSION FOR A SINGLE LAYER 499 


As a convenient notation (Stoneley, loc. cit.) we define 
ri=W(1 —o8/ai), 
si=17(1 —c*/B3), 
b, =(2 +3) /x2? =2 —c2/B2, (1=1,2), 
cosh r,H=Cr,, sinh r,H=Sr,, cosh s,xH=Cs,, sinh $7, =88,, Pe) Us, 
where subscripts indicate the particular medium. 


From the equations of motion and boundary conditions it may be shown 
(Sezawa 1927) that Rayleigh surface waves exist so long as 


)Ne= 
(xb,/7,)Sr,+(xb,/72)Cr, 2Cr, +(27,/72)Sr, 8o/x—x/To 
—(2s,/%x)Ss;—(2x/r.)Cs, —b,Cs,—(x?b,|753,)Ss; 
2b,(1—p)Cr, + (44—26))Cs, (4r,/*)(1—p) Sry .(2—b.) 


+ (xb;/8;)(24—b,)Ss, 


(xb,?/17,)Sry+(%b,b2/72)Cr, 2b,(Cr,—Cs,) u(28,/x—b,%/r2) 
—(48,/x)Ss,;—(2xb,/r,)Cs, + (u.b2/72)[27S7,—(x7b,/81)Ssq] 
=0. 


The derivation implies that x>0, ¢<,<f.<a,<a,, but real solutions 
also exist for B,<c<f,. In this case Cs, is real and Ss, imaginary; but in 
(2) s, always occurs as s,Ss, or Ss,/s,, so that A remains real. When c=8,, 
s,=0 and there is a limiting case but no analytical difficulty arises. 


III. NUMERICAL CALCULATIONS 
For each given value of xH equation (2) was solved for c using the Newton- 
Raphson iteration process. That is to say, if ¢, is the nth iterate, 


oA 
Cn4+1 =0,—A/ 5, 


The process was terminated as soon as two successive values of ¢ agreed 
to five decimal places, and generally only two or three iterations were necessary. 
This speed of convergence is comparable with that for more complicated matrix 
iterative methods (e.g. Haskell 1953). 

In (2) ¢ is a function of xH so that the waves are dispersive. The group 
velocity C, given by 

C=c+xde/dx, 


was calculated directly using 


de _ a oA 
ee wi sof 3o 


and the calculation time for each step was about 5 sec. 
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For each new value of xH the first iterate to be used was found from the 
formula 


¢,(xH +3(xH)) =¢e(xH) — BE seam | ge" 


where 3(xH) is the difference between the two values of xH. The partial 
derivatives required here and in the calculation of C were found by a numerical 
differentiation. 

The value of ¢, for the first value of xH used in any case was found by a 
subsidiary program which examined A for a set of values of ¢ in a specified 
arithmetical progression for specified values of xH, printing whenever a change of 
sign occurred in A as ¢ was increased. 

In the main program xH was printed at each step, followed by c, C, and the 
number of iterations required. All printing was in floating decimal form and the 
calculation time for each step was about 5 sec. Program tapes are available at 
the Department of Applied Mathematics, University of Sydney, for future caleu- 
lations. The elastic parameters can be varied to suit different crustal models 
and still higher modes may be calculated if required. 


IV. DISCUSSION OF THE SOLUTIONS 
Velocities for the fundamental and first higher mode (M,, or Sezawa mode) 
for the models A to K defined in Table 1 are tabulated in Tables 2 and 3. 
A solution for the last case in Table 1, previously obtained using a desk machine 
(Haskell 1953) was used as a program check. Typical period-group velocity 
curves for the case H=35 km are shown in Figure 1. 


TABLE 1 
VELOCITIES (IN KM/SEC) AND THE DENSITY RATIOS DEFINING THE ELEVEN MODELS 


Model A B C D E F G H I J K Haskell 
Oy dS ep O20) 6c 6:0" 1620670902 6-0lgNG 0 eee O 0 mG OmemNG OD 6-14 
By 3°69 3"6, (3:6 324 7 3-8, 3°60 3*6) 83:63 bo bee OHO 3-39 
oy 822 98229) 822) 8e20 852 82) ea oO sod ho eae 8-26 
Bo 4-8 4-8 4-8 4:8 4-8 4-8 4-8 4:8 4-8 4-6 5:0 4-65 


Pol Px 1-296 1-296 1-296 1-296 1-296 1-400 1-207 1-296 1-296 1-296 1-296 1-111 


Fundamental mode solutions exist for 0<xH<oco and as xH-0, c+§,. The 
higher mode M,, has a cut-off period corresponding to c=, and as xH—+oo, 
cf. 

There is always a well-defined minimum to the group velocity for both 


modes ; waves associated with these minima are designated Airy phases (Pekeris 
1948). 
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The models indicate that the Airy wave velocity is sensitive mainly to 8,, 
the shear velocity in the upper layer. Of the 11 models, D has the minimum 
and E the maximum Airy wave velocity for both models. In general M,, is 
affected more by changes in §, than is the fundamental mode. 


TABLE 2 
PHASE AND GROUP VELOCITIES (IN KM/SEC) FOR THE ELEVEN MODELS : FUNDAMENTAL MODE 


Model A Model B Model C Model D 
x os 
c C Cc C c C c Cc 

0-10 4-36925 4-33570 | 4:37124 4-33995 | 4-37303 4-34376 | 4-36905 4-33573 
0-20 4:33757 4-27786 | 4:34205 4-28772 | 4-34603 4-29643 | 4-33799 4-28025 
0-40 4-28428 4-18898 | 4-:29465 4-21185 | 4:30374 4-23155 | 4-28780 4-20051 
0-60 4-24051 4-11804 | 4-25717 4-15332 | 4-27151 4-18304 | 4-24841 4-13952 
0-80 4-20082 4:04330 | 4-22346 4-08879 | 4:24270 4:-12649 | 4-21299 4-07070 
1:00 4-16026 3-94839 | 4-18829 4-00162 | 4-21185 4-04537 | 4:17544 3-97379 
1-50 4-03376 3-58409 | 4-07180 3-64295 | 4:103844 3-69182 | 4-04407 3-54222 
2-00 3-86631 3-15691 | 3-90744 3-19540 | 3-94174 3-22796 | 3-84666 2-98114 
2-50 3-69583 2-91036 | 3-73392 2-92586 | 3:76571 2:93846 | 3-63735 2-67553 
3-00 3.55877 2-86055 | 3-59217 2-86729 | 3-61998 2-87221 | 3-47116 2-63439 
3°50 3-46159 2-90469 | 3-49105 2-91063 | 3-51551 2-91493 | 3-35576 2-70060 
4-00 3.39628 2-97465 | 3-42290 2-98242 | 3-44494 2-98834 | 3:27938 2-78922 
4-50 3-35327 3-04225 | 3-37791 3-05217 | 3-39828 3-06021 | 3-22953 2-87016 
6-00 3-29462 3-17795 | 3-31625 3-19259 | 3-33404 3-20454 | 3-16169 3-02553 
7:00 3-28146 3-22227 | 3-30223 3-23867 | 3-31929 3-25238 | 3-14619 3-07506 
8-00 3-27573 3-24613 | 3-29605 3-26435 | 3-31271 3-27934 | 3-13926 3-10288 
10-00 3°271 Bora 3:292 3-283 3-308 3:298 3-135 3:123 
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TaBLE 2 (Continued) 


Model E Model F Model G Model H 
XH 
c Cc c C c C c Cc 
0-10 4-37287 4-34302 | 4-37381 4-34502 | 4-36865 4-33483 | 4-35471 4-32365 
0-20 4-34492 4-29268 | 4-34705 4-29744 | 4-33698 4-27781 | 4:32589 4-27255 
0:40 4-29901 4-21820 | 4:30415 4-22995 | 4-28495 4-19322 | 4-27984 4-20023 
0:60 4-26230 4-16058 | 4-27085 4-17933 | 4:24311 4-12648 | 4-24405 4-14543 
0-80 4-22951 4-00978 | 4-24122 4-12273 | 4-20521 4-05402 | 4-21206 4-08405 
1-00 4-19615 4-02208 | 4-21006 4-04316 | 4-16601 3-95981 | 4-17849 3-99952 
1-50 4-09273 3-72703 | 4-10196 3-69011 | 4-04177 3-59920 | 4-06552 3-64642 
2-00 3-95656 3-37662 | 3-93835 3-20921 | 3-87758 3-18431 | 3-90405 3-20160 
2-50 8-81521 3-15273 | 3-75765 2-90646 | 3-71114 2-94398 | 3-73234 2-93056 
3:00 3-69775 3-08952 | 3-60786 2-84067 | 3:57702 2-89221 | 3-59147 2-87000 
3°50 3-61207 3-11527 | 3-50093 2-88809 | 3-48143 2-93213 | 3-49074 2-91196 
4-00 3-55334 3-17110 | 3-42907 2-96536 | 3-41684 2-99865 | 3-42275 2-98304 
4-50 3-51416 3-22942 | 3-38180 3-04015 | 3-37408 3-06412 | 3-37784 3-05251 
6-00 3:-46014 3-35244 | 3-31728 3-18843 | 3-31523 3-19682 | 3-31624 3-19291 
7-00 3-44796 3-39308 | 3-30267 3-23694 | 3-30180 3-24109 | 3-30223 3-23855 
8-00 3-44269 3-41526 | 3-29625 3-26243 | 3-29586 3-26470 | 3-29605 3-26449 
10-00 3:°439 3-432 3: 292 3-285 3: 292 3°286 3-292 3-283 
Model I Model J Model K 
“Hf 
c C c C c C 
0-10 4-38563 4-35422 | 4-21602 4-18953 | 4-51961 4-48313 
0-20 4-35618 4-30110 | 4-19144 4-14599 | 4-48544 4-42160 
0-40 4-30771 4:22227 | 4-15220 4-08436 | 4-42939 4-33089 
0-60 4-26880 4-16050 | 4-12169 4-03758 | 4-38469 4-26072 
0-80 4-23363 4-09317 | 4:09441 3-98522 | 4-34454 4-18460 
1-00 4-19705 4-00364 | 4:06580 3-91355 | 4:30292 4-08273 
1-50 4-07744 3-63987 | 3-97039 3-61853 | 4-16585 3-66131 
2-00 3°91049 3-18984 | 3-83482 3-24242 | 3-97327 3-14549 
2-50 3:73534 2-92143 | 3-68830 2-99467 | 3-77421 2-86045 
3-00 3°59281 2:-86488 | 3-56491 2-92496 | 3-61585 2-81537 
3-50 3°49134 2-90944 | 3-47475 2-95220 | 3-50507 2-87404 
4°00 3°42304 2-98161 | 3:41300 3-01126 |.3-43136 2-95'729 
4-50 3°37798 3-05213 | 3-37179 3-07212 | 3-38313 3-03522 
6-00 3-31626 3-19289 | 3-31466 3-19934 | 3-31759 3-18703 
7:00 3:°30224 3-23951 | 3-30156 3-24254 | 3-30280 3-23606 
8-00 3°296 3°263 3-29576 3:26585 | 3-29630 3-26285 
10-00 3° 292 3:284 3°292 3°285 3:292 3°284 


Variation in the density contrast between the media (as in Models F and G) 
has only a slight effect on the M,, curve but for 3:7 <e<4-1 km/sec a change in 
P2/e, of 0-2 causes 15 per cent. displacement in the fundamental mode curve 
relative to M,,. The seismological implications of this point have been discussed 
by Oliver, Dorman, and Sutton (1959). 
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An important further inference from the solutions is that different combina- 
tions of elastic parameter values give very similar dispersion curves in both modes 


(e.g. models A and G). 


The main conclusion from a comparison is that observed 


Rayleigh wave dispersion alone cannot be used to give precise values to individual 
elastic parameters. 


PHASE AND GROUP VELOCITIES (IN 


TABLE 3 


KM/SEC) FOR THE ELEVEN MODELS: FIRST HIGHER MODE 


Model A Model B Model C Model D 
yah 
c C c C Cc C c Cc 
3-00 4-73368 4-33063 | 4-73399 4-33471 | 4-73427 4-33834 | 4-65065 4-19180 
3-50 4-67200 4-27672 | 4-67321 4-28524 | 4-67426 4-29263 | 4-58331 4-15835 
4-00 4-61882 4-21087 | 4-62110 4-22184 | 4-62309 4-23134 | 4-52516 4-06725 
4-50 4-56798 4-10360 | 4-57133 4-11615 | 4-57423 4-12704 | 4-46595 3-90390 
5-00 4-51395 3-94197 | 4-51824 3-95487 | 4-52198 3-96610 | 4-39814 3-65893 
5-50 4-45254 3-72817 | 4-45755 3-73920 | 4-46191 3-74861 | 4-31771 3-36665 
6-00 4-38253 3-49883 | 4-38787 3-50566 | 4-39251 3-51138 | 4-22733 3-11285 
6-50 4-30691 3-31110 | 4-31220 3-31401 | 4-31678 3-31590 | 4-:13487 2-95405 
7-00 4-23119 3-19455 | 4-23619 3-19442 | 4-24049 3-19406 | 4-04753 2-88196 
7-50 4-15999 3-14068 | 4-16459 3-13878 | 4-16854 3-13750 | 3:96901 2-86468 
8-00 4-09571 3-12762 | 4-09990 3-12515 | 4-10350 3-12336 | 3:90025 2-87681 
8-50 4-03901 3-13851 | 4-04281 3-13635 | 4-04607 3-13334 | 3-84082 2-90484 
9-00 3-98959 3-16128 | 3-99303 3-15841 | 3:99598 3-15726 | 3-78972 2-93813 
9-50 3-94677 3-19045 | 3-94987 3-18651 | 3-95254 3-18601 | 3-:74585 2-97296 
10-00 3-909 3°225 3-912 3°219 3:914 3: 220 3:708 3-009 
15-00 3-721 3°4 3-722 3°5 3° 724 3°4 3-519 3°2 
Model E Model F Model G Model H 
c C (2 Gi Cc Gs c C 
3-00 4-72245 4-31980 | 4-74448 4-35152 | 4:73147 4-32609 
3-50 4:-74510 4-43259 | 4-66313 4-29467 | 4-68249 4-27572 | 4-67004 4-27940 
4-00 4-70198 4-36900 | 4-61445 4-24651 | 4-62697 4-19597 | 4-61772 4-21799 
4-50 4:66109 4-29493 | 4-56892 4-15446 | 4-57289 4-07662 | 4-56802 4-11449 
5-00 4-61956 4-19071 | 4-52051 4:00431 | 4-51517 3-90599 | 4-51523 3-95598 
5-50 4-57439 4-04808 | 4-46427 3-78974 | 4-45028 3-69187 | 4-45505 3-74317 
6-00 4-52340 3-87328 | 4-39773 3-54141 | 4-37783 3-47400 | 4-38600 3-51136 
6-50 4-46642 3-69478 | 4-32320 3-32741 | 4-30129 3-30201 | 4-31091 3-31930 
7-00 4-40583 3-54929 | 4-24674 3-19175 | 4-22578 3-19715 | 4-23533 3-19863 
7:50 4-34534 3-45600 | 4-17397 3-12796 | 4-15533 3-14823 | 4-16401 3-14120 
8-00 4-28810 3-40963 | 4-10795 3-11294 | 4-09192 3-13726 | 4:09951 3-12704 
8-50 4-23587 3:39548 | 4-04962 3-12319 | 4:03602 3-14822 | 4-04253 3:- 13763 
9-00 4-18927 3-40119 | 3-:99878 3-14714 | 3-98728 3-16862 3-99283 3-15998 
9-50 4-14820 3-41636 | 3-95473 3-17593 | 3-94500 3-19656 3-94972 3:-18784 
10-00 4-111 3-443 3-916 3-215 3°908 3: 232 3°912 3-223 
15-00 3-924 3°6 3°722 3-4 3°721 3°5 3-722 3-5 
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TABLE 3 (Continued) 


Model I Model J Model K 
“x 
c C c C c C 
3-00 4-73621 4-34247 | 4-58280 4-33042 | 4-86340 4-37757 
3°50 4-67602 4-29048 | 4-53799 4-22874 | 4-79202 4-34387 
4-00 4-62410 4:22530 | 4-49579 4-17135 | 4-73148 4-26216 
4-50 4-57428 4-11765 | 4-45582 4-09594 | 4-67191 4-11676 
5-00 4-52093 3-95383 | 4-41464 3-98539 | 4-60604 3-89748 
5-50 4-45978 3-73546 | 4:36897 3-83207 | 4-52942 3-62456 
6-00 4-38953 3-50057 | 4-31661 3-64609 | 4-44289 3-36702 
6-50 4-31335 3-30895 | 4-25783 3-46294 | 4-35269 3-18887 
7-00 4-23695 3-19085 | 4-:19572 3-32289 | 4-26584 3-09738 
7°50 4-16510 3-13632 | 4-13442 3-23896 | 4-18668 3-06728 
8-00 4-10025 3-12350 | 4-07703 3-19963 | 4-11674 3-07172 
8-50 4-04305 3-13463 | 4-02511 3-19233 | 4-05593 3-09614 
9-00 3:99321 3-15767 | 3-97907 3-20278 | 4-00346 3-12762 
9-50 3:95001 3:18687 | 3-93866 3-22118 | 3-95830 3-16210 
10-00 3:912 3:223 3-903 3°244 3-918 3-205 
15-00 3-722 3°5 3-720 3:4 3°724 3-4 
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THE PHOTODISINTEGRATION OF RARE EARTH ELEMENTS 
By H. H. Tures*} and B. M. Spicer* 
[Manuscript received March 10, 1960] 


Summary 


The neutron yields from praseodymium, terbium, and holmium have been measured 
as a function of bremsstrahlung energy, and from these the shapes of the respective 
cross sections have been deduced. For terbium and holmium, the cross section for 
photon absorption shows two peaks, in agreement with the predictions of both the 
hydrodynamic model and the independent particle model of the nuclear photoeffect in 
non-spherical nuclei. The case of praseodymium was treated as a control experiment, 
since this nucleus is nearly spherical. Its photon absorption cross section has only 
one peak, in agreement with expectation for a spherical nucleus. The results are 
discussed in terms of both the hydrodynamic model and the independent particle model. 
Values of the deformation parameter are deduced from the photodisintegration data 
and are compared with values obtained by other methods. 


I. INTRODUCTION 

Both collective theory (Danos 1956, 1958; Okamoto 1956, 1958) and shell 
model theory (Wilkinson 1956, 1958) predict that in non-spherical heavy nuclei 
the photonuclear giant resonance is broadened and will show two maxima if the 
nuclear deformation is large enough. Known nuclear electric quadrupole 
moments indicate that the nuclear deformation increases from a negligible amount 
for elements of atomic mass near 140 to a maximum for A near 170, and decreases 
again to small values for A near 200. Thus, by systematically investigating the 
shape of the giant resonance for nuclei over this mass range, it should be possible 
to demonstrate experimentally any broadening and splitting effect in the giant 
resonance. This would then lead to the testing of theories on the subject, and 
possibly to the use of photodisintegration data for the evaluation of nuclear 
deformations and electric quadrupole moments. 

In the heavy elements it is not practicable to observe the photonuclear 
absorption cross sections directly. However, to a fair degree of approximation, 
they can be derived from photoneutron cross sections for these elements. The 
approximation here arises because of the need to correct for multiple neutron 
emission. Because of the very high Coulomb barrier, the emission of charged 
particles gives negligible contribution. 

If the photoneutrons are to be observed directly, only mono-isotopic or 
almost mono-isotopic elements may be used as targets if separated isotopes are 
not available. In the region of interest here these are all odd-mass nuclei, 
namely: 1%La, 141Pr, 159Th, 165Ho, 1° Tm, 1757 y, 18'Ta, 97Au, Bi, 
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In a previous paper (Spicer et al. 1958) an experiment on the photonuclear 
giant resonance in La and 18!Ta was reported. This experiment used direct 
neutron detection (BF -counters in a paraffin moderator). Simultaneously a 
similar very accurate experiment on 1°Tb, 1§'Ta, %7Au was performed inde- 
pendently by Fuller and Weiss (1958). In the present paper further direct 
neutron detection experiments on ™!Pr, 1°°Tb, and 16°Ho are reported. 


This type of experiment measures the total neutron production cross section 
Sy,In= Syn t25y,2n +30y,3n ++ + - 


The computation of the absorption cross section from o, >, above the (y,2n) 
threshold requires a knowledge of the ratio S=oy,,/oy,2. S may be estimated 
from statistical theory of nuclear reactions (Blatt and Weisskopf 1952). For 
the nucleus 165Ho it was possible to determine S experimentally by combining 
the results from the direct neutron detection experiment with those from an 
independent activation experiment on the reaction 1®°Ho(y,n). 

The available experimental results on 1%La, 141Pr, 15°Tb, 1®°Ho, 1%!Ta, and 
197Au will be discussed in the light of theoretical predictions. 


IJ. THEORETICAL DESCRIPTION OF THE PHOTONUCLEAR EFFECT 
The predictions of the hydrodynamic theory have been discussed previously 
(Spicer et al. 1958) and lead to the result that the two eigenvalues are related by 


the equation 
LH ,/E, =hw,/io,=0-911a/b +0 -089. 


Here a refers to the semi-axis of rotational symmetry and b to the semi-axis 
perpendicular to a. If (H,—EH,) is large enough, the positions of H, and #, 
will coincide approximately with the positions of the experimental maxima 
(#, and H,). Frictional coupling between the two eigen-modes will tend to 
decrease (#,—H,) and flatten the minimum between #, and £,. 

Wilkinson (1958) used the Nilsson level scheme and computed photonuclear 
absorption cross sections for 140<A<210. His results show that cross sections 
in non-deformed nuclei have one narrow maximum. This is broadened in 
deformed nuclei, and for A—165 the computed cross section is split into two 
separated maxima. 


Wave functions and transition probabilities corresponding to the Nilsson 
Hamiltonian can be computed (see Nilsson 1955). They do, however, depend 
much more sensitively on the form of the chosen Hamiltonian than do the 
transition energies, and thus the computation of transition probabilities is not 
significant for excitation energies in the region 10-25 MeV. It is possible, however, 
to predict semiquantitatively the manner in which the nuclear deformation 
will affect electric dipole absorption cross sections, without actually knowing the 
Strength of the #1 transitions individually. In the spheroidal nucleus, nuclear 
States are specified by different quantum numbers than in the spherical case. 
Thus, new selection rules (Alaga 1955, 1957) become increasingly effective with 
increasing nuclear deformation. In deformed nuclei, these selection rules 
separate the allowed electric dipole transitions into two groups. In transitions 
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of group I, the excitation energy is used predominantly to add energy to the 
component of motion of one nucleon in the direction of a, the axis of rotational 
symmetry. In group II, the energy is added predominantly to the component 
of motion in a direction perpendicular to a. 


The systematic shift of levels with deformation of a prolate type causes 
excitation energies for group I transitions to be smaller than for group II 
transitions, in general. (For an oblate spheroid, the opposite condition tends 
to hold.) Also, as a consequence of the selection rules, there are approximately 
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Fig. 1—Cross section for photon absorption in 1°5Ho computed from the 
Nilsson level scheme. The three cross sections are for different assumed. 
nuclear deformations. 


twice as many allowed £1 transitions in group TI as in group I. Thus the #1 
absorption cross section will separate into two giant resonance components. 
These components will have mean energies of H, and Hy, respectively. This does 
not mean that the maxima in the composite cross section will occur at H, and Ey, 
(see Fig. 1). For deformations insufficient to give rise to separation of the peaks, 


a single broadened peak will result. as 
For particular nuclei, we investigated the broadening and splitting effect 
From the Nilsson level diagram (Nilsson 1955) the 


in the following manner. 
i rie di iti the selection rules 

itati 3 of the electric dipole transitions allowed by 
Sale ae it is assumed that all the allowed 


are found. In the first approximation, 
transitions have equal strength. The corresponding absorption lines are given 
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half-width yey, and a Gaussian shape. Contributions from the various 
transitions are summed at a series of intervals. The value of Vicye: may be taken 
as 3 MeV, since this is an estimated value for non-deformed nuclei given by 
Wilkinson (1956). However, it is possible that this value may be less than 3 MeV 
for strongly deformed nuclei, as in this case the level degeneracies are reduced 
considerably. Component H1 cross sections o, and” oj, corresponding to 
transitions in groups I and II, may be derived separately. Although the strengths 
of transitions may vary considerably from equality, these variations will not 
affect the shape and position of o,; and oj; significantly, as in this region of A 
the number of transitions is large. For example, for A near 160, there are 
14 proton transitions in group I and 28 in group II. The cross sections o; and oj;, 
and their sum oaps, for the nucleus 1®°Ho are shown on Figure1. The calculations 
were made for three values of the deformation. 

It was found empirically from a number of diagrams similar to Figure 1 
that the positions of H, and H,,, the mean energies of co; and o;;, follow the 
relationship 


Ei fg LEO 1 Co ew sin ane ew Ce ee (1) 
Here ¢ is the deformation parameter (Danos 1958) defined by 
6==(G* =O") (ts Ss pss was sn cr (2) 
It is related to the deformation parameter 5 of Nilsson (1955) by the relationship 
©=20(1/-- 02666). fe. we ee eee (3) 


R is the radius of a sphere having the same volume as the spheroid. 
(R3=ab?=r3A. Throughout this paper, the value r,=1-2 fermi is used (Ford 
and Hill 1955).) Thus the independent particle model predicts that the giant 
resonance is broadened in non-spherical nuclei, and, for large deformations, 
splits into two maxima in a manner similar to that predicted by the hydro- 
dynamic model. 

Using equation (2), the prediction of the hydrodynamic theory may be 
expressed in terms of an expansion in ¢ as follows ; 


L,/E,=1+0-46e+terms of order (e?). 


In view of the crudeness of the assumptions which led to (1) and (4), the pre- 
dictions on the splitting effect from hydrodynamic theory and from independent 
particle model theory are in substantial agreement. This is so in spite of the 
very different assumptions in the two theories. 


III. EXPERIMENTS ON PHOTONEUTRONS FROM 
id St cers AND 165Ho 

Bremsstrahlung from the Melbourne 18 MeV electron synchrotron was used 
to irradiate the targets. Each target consisted of 10 g of a rare earth oxide in 
powder form (Pr,O,,, Tb,0,, Ho,O,). The oxide was contained as a tablet in a 
polystyrene container. Sample impurities consisted of neighbouring rare earths, 
and were smaller than 1 per cent. in each case. The absorption of y-rays in 
the target and the container lid was less than 3 per cent. in all cases, and was 
taken into account in the yield analysis. 
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With the exception of an activation experiment on the reaction 
'S>Ho(y,n), the photoneutrons were detected by BF,-counters embedded in a 
paraffin moderator. The arrangement of target, BF,-counters, paraffin 
moderator, and concrete ‘ house ’’ was essentially the same as that described 
by Spicer et al. (1958). Four BF,-counters, separated from the target by 12-5 em 
of paraffin, were grouped as two pairs, and the counts from these were recorded 
through two independent electronic counting systems. The electronic counting 
loss in each system was determined experimentally and was kept less than 
1-2 per cent. by limiting the yield rate of the synchrotron when necessary. 
Corrections were made for the counting loss. The neutron detection efficiency 
of each counting system was checked before and after each run by inserting a 
Ra-Be neutron source into an accurately reproducible position in the paraffin 
moderator and recording the counts obtained per minute. The number of 
neutrons counted in a given run was thereby corrected for the average deviation 
in neutron detection efficiency from a standard average value. 

The y-ray dose was measured by using a 25 r Victoreen thimble as a trans- 
mission ionization chamber. The leakage of this instrument was measured 
periodically and found to be negligible over several days. Hach run was arranged 
so as to make the recorded dose approximately 10 r, thus effectively eliminating 
any relative error that would arise from a calibration error in the Victoreen 
scale. Temperature and pressure were recorded before and after each run and a 
temperature-pressure correction was applied to each dose reading. 

The energy scale of the synchrotron had been calibrated previously by means 
of photoneutron threshold measurements (Spicer et al. 1958). This calibration 
was checked periodically by measuring the **Mn(y,n) threshold. 

For proper interpretation of the experimental results, the standard error 
in each yield point must be known. This error may be substantially greater 
than the purely statistical counting error. It can be estimated only inadequately 
from the reproducibility of neutron counts, y-ray dose measurements, and 
energy settings of the accelerator, but can be obtained accurately if a yield 
ordinate has been measured many times at random time intervals during a 
complete experiment. 

During each experiment, individual yield ordinates were measured in some 
instances up to 14 times, and it was found that the spread of values obtained at 
each energy setting followed a normal distribution within statistical accuracy. 
Experiments were performed with praseodymium, terbium, and. holmium targets 
in position and background runs were made with the empty containers in the 
target position. The background remained unaltered when the empty poly- 
styrene container was removed. 

The contribution of photoneutrons from the oxygen content of the rare earth 


oxide samples was estimated. Under all circumstances its contribution to the 


total neutron production was less than 1 per cent. 

Yield ordinates were measured at 0°25 MeV intervals from 7 to 18 MeV. 
These results are shown in Figures 2, 3, 4, for praseodymium, terbium, and 
holmium respectively. Smooth yield curves were drawn through each discrete 
set of ordinates, fitting by eye with minimum curvature curves. For the 
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Fig. 2.—Total neutron yield from praseodymium. 
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Fig. 3.—Total neutron yield from terbium. 
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transformation of yield data to cross sections for photoneutron production the 
tables of Penfold and Leiss (1958), having AE=1 MeV, were used. Each set 
of measured yield ordinates was separated into four subsets, each subset having 
AE=1 MeV. These subsets were then analysed separately, and the resultant 
cross sections were then combined to make a composite cross section. Thus, the 
cross section ordinates making up one ‘unsmoothed”’ cross section curve come 
from four interlaced and independent sets of ‘ unsmoothed ”’ cross section 
ordinates. Standard errors of unsmoothed cross section values were computed 
from standard errors of yield ordinates. 
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Fig. 4.—Total neutron yield from holmium. 


The smooth cross sections shown were obtained from analysis of the smoothed 
yield curves as drawn through the experimental points. Values of the yield 
were read off this smooth curve and these were treated in the same manner as 
the unsmoothed yield points. The cross sections are shown in Figures 5 (1“1Pr), 
6 (°Tb), and 7 (16°Ho). The results on terbium and holmium are the averages 
from two independent experiments. These results agreed within the standard 
errors and demonstrated the reproducibility of results. Note that these cross 


sections are the total neutron production cross sections. 


To reconstruct the photon absorption cross section Gaps=Sy,n+y,2n from 
the total neutron production cross section oy,0n=Syn+26y,2, allowance has to 
be made for neutron multiplicity above the (y,2n) threshold. 
thresholds of some rare earths have recently been measured by Geller, 


ae Muirhead, personal communication 1959), and. agree 


Muirhead, and Halpern ( 
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within less than 0:3 MeV with those calculated from the semi-empirical mass 
formula of Cameron (1957). Therefore this formula was used to calculate the 
values of (y,2) thresholds. The values obtained were 17-6 MeV for 141Pr(y,2n), 


CROSS SECTION (MBARN) 
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Fig. 5.—Total neutron production cross section for 41Pr. The points 

are obtained by analysis using the experimental points. The smooth 

curve is obtained by analysis of a smooth yield curve drawn through the 
experimental points. 
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Fig. 6.—Total neutron production cross section for 15°Tb. 


14-8 MeV for *°Tb(y,2n), and 14-4 MeV for 1®°Ho(y,2n). Thus, in the present 
experiment on 141Pr, the neutron multiplicity need not be considered. 

Blatt and Weisskopf estimated the neutron multiplicity using statistical 
theory of nuclear reactions. Applied to the present instance, this does not take 
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account of the neutrons emitted in a direct photoeffect. This can be built in ‘to 
the statistical theory estimate, which then becomes lis * 


Syn +2 0y,9n 
1+(1 —a#){1 —(1+k/@)e-/O} 


Here x is the fraction of direct events to total number of events, k is the energy 
difference between the incident photon energy and the (y,2n) threshold value, 


Oy,n + Oy,2n = 
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Fig. 7.—Total neutron production cross section for 1°°Ho. 
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Fig. 8.—Cross section for photon absorption in 1°°Tb, using statistical theory 
of nuclear reactions to correct for neutron multiplicity. The factor x is a measure 
of the probability that a direct photoeffect will occur. 


and © MeV is the nuclear temperature of the target nucleus less one neutron. 
Values of @ were obtained from Blatt and Weisskopf (1952). 

Equation (5) can be brought into agreement with measured data on neutron 
multiplicities if a value of # between 0-3 and 0:4 is assumed. Experimental 
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data considered are by Montalbetti, Katz, and Goldemberg (1953), Whalin and 
Hanson (1953), Berman and Brown (1954), Silva, Goldemberg, and Smith (1958), 
and Carver and Turchinetz (1959a, 1959b). 

The cross sections in 1°Tb and 16°Ho, corrected for neutron multiplicity, 
are shown in Figures 8 and 9. Curves are drawn for values of # of zero and 0°3. 
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Fig. 9.—Cross section for photon absorption in 1°°Ho. The 
correction for neutron multiplicity is applied using statistical 
theory. 


IV. ACTIVATION EXPERIMENT 

The (y,) cross section in 1®°Ho was determined from a separate experiment 
by counting the ®--activity of the daughter nucleus Ho. Brown and Becker 
(1954) give the half-life of this nuclide as 36-7 min. The $--spectrum has a 
maximum energy of 0:99 MeV. 

The Ho,O, used in the BF,-counter experiment was subdivided into nine 
samples of 1-050 ¢ each. These were each pressed into tablet form (1-4 ecm 
diameter by 0-18cm) in identical polystyrene containers, and protected 
by 0-32 cm thick ‘ Perspex ”’ lids. The containers could be located accurately 
and reproducibly in a jig 30 cm from the synchrotron target. Dose measurements 
were made in the manner previously described. In this experiment the y-ray 
beam traversed the target before striking the Victoreen thimble. The absorption 
of y-rays in the target was therefore allowed for in the analysis. 

The container lid was removed for counting the activity. This was com- 
menced 5 min after the cessation of irradiation (to allow the 15O 2-min activity 
to die out to some extent) and continued for 30 min. Only pure 37-min activity 
of 1**Ho could be detected. Allowance was made for the finite dead-time of the 
Geiger counting system, and the appropriate correction applied to the results. 

A procedure was developed to allow for the decay of the activity during 
irradiation, and for variations in the yield rate of the synchrotron. The y-ray 
beam, after traversing the dosemeter, was allowed to strike a 500 g lead sample. 
This was located in the “ sample position”? in the paraffin moderator plus 
BF;-counter assembly that had been used for the direct neutron detection 
experiments. The counters therefore received copious photoneutrons from the 
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lead sample. The number of neutrons recorded in any time interval is pro- 
portional to the y-ray dose given in that time. The neutron counts were recorded 
at 2-min intervals and used as a relative measure of the dose which the holmium 
sample had received during that interval. The decay of the activity of the 
holmium sample during the 30 min irradiation period could then be calculated 
with due allowance for variations in dose rate. 

Activation points were taken at four energies below the (7,2n) threshold, 
and at 0-25 MeV intervals above 14-5 MeV. Each point was repeated at least 
twice. A smooth curve was fitted to the measured points (see Fig. 10). This 
curve was normalized to the total neutron yield curve at the energies below the 
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Fig. 10.—Yield curve for the 165Ho(y,n)!°4Ho reaction, measured 
by counting the 1*Ho activity. It was normalized on to the total 
neutron yield curve for holmium below the (y,2n) threshold. 


(y,2”) threshold at which the activity was measured. The (v.72) cross section 
at energies above the (7,2n) threshold was computed from the yield curve in 
the manner previously described, and is shown a8 oy, On Figure 11. On the 
same figure the total neutron cross section oy,En is shown. It is to be noted 
that the yield curve for the 1°*Ho activity has not been measured well enough for 
analysis using the experimental points to be useful. Consequently, all SESE 
of this cross section refers to the result of analysis of the smooth yield curve. 
Using the total neutron cross section oy, Xn given by the BF,-counter 
experiment and the (y,n) cross section of the activation experiment, the photon 


absorption cross section is given by 
1 
GOabs = Oy,n + Oy,2n = 3{Oy,En + Oy,n5- 
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The absorption cross section cap; computed with equation (6) is compared in 
Figure 12 with the absorption cross section computed using neutron multiplicity 
theory, using equation (5) with ~=0-3. The two cross sections are in agreement 
within statistical expectations, although a value of x as high as 0-4 would be 
equally suitable. 
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Fig. 11.—Result of analysis of the holmium activation 

experiment. Also plotted are the total neutron production cross 

section and the photon absorption cross section derived using 
the activation results. 


800 


600 


400 


a(MBARN) 


200 


8 10 12 14 16 18 
ENERGY (MEV) 


Fig. 12.—Comparison of the calculated cross section for photon absorption in 
16°Ho with experimental measurements of the same quantity. Curve I has the 
neutron multiplicity correction applied using the activation experiment. Curve II 
has neutron multiplicity correction applied using statistical theory of nuclear 
reactions, with an assumed 30 per cent. direct photoeffect. 


V. RESULTS 
The derived y-absorption cross section for the undeformed nucleus 141Pr 
showed one narrow giant resonance maximum, whereas in the strongly deformed 
nuclei +°T'b and 16°Ho the predicted broadening and splitting effect was observed. 
These results are consistent with those obtained for the spherical nuclei 13%La 
and "Au, and the deformed 1°Tb and 18!Ta, by Fuller and Weiss (1958) and 
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Spicer et al. (1958). In both sets of experiments essentially the same experi- 
mental technique was used. 

When comparing the experimental results on photonuclear absorption cross 
sections with theoretical predictions, the notation given below will be used. 
The component of the giant resonance will have maximum cross section (c)max. 
at energy H. The width of the resonance will be discussed in terms of [_; 
and I’,,;, which give the energy difference between H and the energy at which 
the cross section is half its maximum value below and above the energy H 
respectively. The total width of the resonance, I’, is given by ((C-y+/'43). 
The integrated cross section J is in good approximation proportional to the 
product of (c)max, and [. Thus the ratio of the integrated cross sections under 
two component peaks of a split giant resonance is given by 


I, _(69)max. ib 


ae ee (62) max. T,, +4 
I, (01)max. |i (6;)max. ee 


The total width of a double-peaked giant resonance shall be defined by 
Tern=l1-3~tPe4g they. cece ee ceeee (8) 


The giant resonance parameters for La, 1Pr, 7Au, A592 Elo, ands 228s 
are collected in Tables 1 and 2. 


TABLE 1 
GIANT RESONANCE PARAMETERS FOR NON-DEFORMED NUCLEI 


Nucleus. . 1397, AEP SOA 
Reference Spicer et al. Present Paper | Fuller and Weiss 
(1958) (1958) 
E (MeV) 14-9 15-1 13-5 
(C)max, (mbarn) 460 480 590 
2xI_; (MeV) 4-2 4-4 4 
Ter. (MeV) 4-2 — 4-2 


TABLE 2 
GIANT RESONANCE PARAMETERS IN DEFORMED NUCLEI 


Nucleus SIAM 8) 165Ho 181T a, 
i L 
Reference Fuller and Fuller and Spicer et al. 
Weiss (1958) Present Paper | Present Paper Weiss (1958) (1958) 
Tqr. (MeV) 6:5 6°7+0°-5 740-5 6-2+0-2 6-1 
E, (MeV) 12-5 12-440-2 12-140-2 12-45 12-6 
(6;)max. (mbarn) 260 410 420 308 500 
IT, (MeV) 2-4 3°3 2-8 2-3 2-0 
E, (MeV) 16-3 16-0-40-2 16-20-2 15*45 15-3 
(62)max, (mbarn) 310 460 510 348 a ; 
T, (MeV) 4-0 4-5 4-7 4-4 
I,/I, 2-0 2-0 2-0 2-16 1:8 
#,/E 1-30 1-29-L0-03 1-34+0-03 1-25+0-01 1-2140-03 
2 1 
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VI. INTERPRETATION OF RESULTS 

If the experimentally determined values of H, and F, are equated to the HE, 
and H, respectively of the Danos-Okamoto theory, the nuclear deformation 
parameter may be calculated from equation (4) by inserting the relevant values 
from Table 2. These values are entered in Table 3. 

The independent particle model theory may also be used to obtain values of 
the nuclear deformation. The shapes of the cross sections were calculated 
from the energy level diagram of Nilsson (1955), and the value of ¢ was estimated 
by comparing the computed cross sections with the experimental ones. The 
comparison for the case of 1°°>Ho is shown in Figure 12. 


TABLE 3 
NUCLEAR QUADRUPOLE MOMENTS AND DEFORMATION IN }°%rb, 16°Ho, AND 181Ta 


Method Nucleus € Qo Reference 
(barn) 
Coulomb excitation ome D4 Bs) 0:69 8-1 Martin, Marmier, and de Boer 
(1958) 
Coulomb excitation ne 0-63 6-9 Alder et al. (1956) 
Survey, low-lying levels .. 0-74 ~- Mottelson and Nilsson (1959) 
Photodisintegration (Danos) 0-63 7:0+0°7 Present paper 
Photodisintegration (IPM) 0-72 8-0+0-8 Present paper 
Paramagnetic resonance .. EEO 0-4 +4-3 Baker and Bleaney (1955) 
Coulomb excitation ie 0-69 8-0 Martin, Marmier, and de Boer 
(1958) 
Coulomb excitation 67 7:8 Alder et al. (1956) 
Survey, low-lying levels .. “72 Mottelson and Nilsson (1959) 


Photodisintegration (Danos) 
Photodisintegration (IPM) 


74 8:6+0-9 Present paper 
“72 8:4+0-9 | Present paper 


ooo & 


Atomic spectroscopy in ae a 0-43 |+5-8+0-5 Murakawa and Kamei (1957) 

Coulomb excitation 0-41 5-6 Martin, Marmier, and de Boer 
(1958) 

Coulomb excitation 0-50 6-8 Alder et al. (1956) 

Survey, low-lying levels .. 0-53 Mottelson and Nilsson (1959) 

Photodisintegration (Danos) 0:54 7-0 Fuller and Weiss (1958) 

Photodisintegration 0-46 6-0 Spicer et al. (1958) 


If the nuclear charge is assumed to be uniformly distributed throughout the 
nucleus, the values of ¢ obtained may be used to compute a value of Q,, the 
intrinsic nuclear electric quadrupole moment. The quadrupole moment of a 
charge Ze uniformly distributed throughout a spheroid is given by 


Qa =0 4216202). eee ee ain ee (9) 
The relation between < and Q, is 
c=2 BO (ree ate eet ee (10) 


using equation (2) and R’=r34. In Table 3 the values of ¢ and Q,) obtained 
from the photodisintegration data are set out. Also shown in Table 3 are values 


the strong coupling approximation in the case of the nucle 
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for « and @,) obtained from measurements of atomic spectra, from Coulomb 
excitation measurements, and from a survey of the low-lying levels of deformed 
nuclei. 

The theory of fast electric quadrupole transitions between low-lying levels 
(Coulomb excitation) and of the spacings of low-lying levels in deformed nuclei 
are based on the collective theory of nuclear rotational states (Bohr and Mottelson 
1953). Both of these types of study give values of wr directly in terms of the 
collective theory, and therefore do not specify the sign of Q. The accuracy of 
such Q-values is limited to some 20 per cent. still, due to the inherent difficulties 
in the absolute measurements necessary. To compute the deformation parameter 
¢, again some sort of assumption has to be made concerning the nuclear charge 
distribution. Again, we assume a uniform charge distribution. 
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Fig. 13.—Plot of [gp and « as a function of neutron number over 
the region of interest. 


The atomic spectroscopy measurements do not give the value of Q 
immediately. Rather they give a value of the quadrupole moment with respect 
to the nuclear angular momentum direction, say Q, which is less than the value of 
Qo. Qo is the quadrupole moment relative to the direction of the nuclear axis 
of symmetry, which in strongly deformed nuclei is not necessarily in the same 
direction as the angular momentum vector. In the extreme case of strong 
coupling (Bohr and Mottelson 1953) Q and Q, are related by Q =();./,, where 


~ 1(2—1) 
1 (£+1)(21 +8) 


i i 1958) has argued for the validity of 
and J is the ground state spin. Okamoto (1958) 2 seer ie aeaneTe 


yg 


It must be pointed out that the comparison of the above data with photo- 
disintegration data should be treated with some caution. The three alternative 
methods of determining nuclear deformations which were mentioned are all 
concerned with nuclei in their ground state or in states very near it. In contrast, 
photodisintegration data concern transitions between the ground state and 
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excited states 15-20 MeV above the ground state. Again, Okamoto (1958) has 
argued that the comparison made here is worth while, and the close agreement 
exhibited in Table 3 indicates that his argument may be correct. 


It is to be noted that the ratios J,/Z, noted in Table 2 do approximately 
fulfil the theoretical prediction of the value 2 for that ratio. 


On Figure 13 the values of « and of Tg, are plotted as a function of the 
neutron number N. It is seen that the broadening effect in the photonuclear 
giant resonance is dependent on nuclear deformation. 


VII. CONCLUSION 
The quantitative investigation of the broadening and splitting effect in the 
photonuclear giant resonance, in particular the latter effect, requires extremely 
accurate experiments. The splitting effect remained unobserved in earlier 
experiments because in those cases the experimental resolution was not sufficient. 


The splitting effect was not observed in the (y,n) cross sections of 1%La, 
141Pr, and 7Au. This is in agreement with the prediction that the effect will 
not be present for non-deformed nuclei. This does not exclude the possibility 
of some fine structure in these cross sections, but merely states that the experi- 
mental resolution at present attainable is insufficient to detect it. 


The experimental resolution in the experiments of Fuller and Weiss (1958), 
Spicer e¢ al. (1958), and of the present work was sufficient to indicate that the 
splitting effect observed is real. A discussion of the resolution of bremsstrahlung 
experiments is given by Thies (to be published). The predictions on photo- 
nuclear absorption cross sections from hydrodynamic theory and from 
independent particle model theory are very much alike, and are in excellent 
agreement with observation. It was not possible to decide whether experimental 
cross sections favour either one of the two theories. 
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SOUTHERN HEMISPHERE METEOR SHOWER ACTIVITY IN 
JULY AND AUGUST 


By A. A. WEISS* 
[Manuscript received March 21, 1960] 


Summary 
Observations of the radiant positions and the activities of the $-Aquarids and 
three minor showers of July and August, obtained with high resolution radio equipment 
at Adelaide from 1957 to 1959, are given in detail. Data for the §-Aquarids are combined 
with similar observations in 1953 to obtain the radiant ephemeris 
a= 342-3°+0-76 (@—126-0°) 
d5=—16-4°+0-13 (Q—126-0°). 
The §-Aquarid activity extends from July 20 to August 14, with a maximum near 
July 29. Minor shower activity, extending over and beyond the whole duration of 
$-Aquarid activity, is identified as proceeding from three main centres already known 
visually : the Capricornids, the Pisces Australids, and the Cetids. 


I. INTRODUCTION 

The period from the beginning of July to the middle of August is one of 
high meteor shower activity in the southern hemisphere. The principal shower 
is the 6-Aquarids. Numerous minor showers with radiants south of the ecliptic 
have been detected by visual observations (McIntosh 1935). In addition, a 
previously unreported day-time shower, the Phoenicids, has been discovered by 
the Adelaide radio surveys. 

With the exception of 1955, observations over July and August have been 
made each year from 1953 to 1959. The observations prior to 1956 have already 
been published (Weiss 1955, 1957), and all available data on the Phoenicid 
shower have been summarized in Weiss (1960). The present paper describes the 
remaining observations for 1957 to 1959, and summarizes the measurements, 
with high resolution equipment, of the radiant coordinates and echo rates for 
the $-Aquarids and three of the more prominent of the minor showers. The 
observations for 1954 and 1956, made with the 27 Mc/s ‘ wind ”’ equipment, 
lack the resolution necessary for the separation of the $-Aquarid from the minor 
shower activity, and have been excluded from this summary. 


Il. EQUIPMENT 
The bulk of the observations have been made with the 67 Me/s radiant 
survey equipment (Weiss 1955, 1960). This is a high resolution radar system 
with two narrow-beam aerials directed at low elevation along azimuths 35-5° 
S. of H. (termed the N aerial) and 14° N. of E. (S aerial). The line densities of 
electrons in the faintest trails detected by this equipment are close to 11012 
electrons/em, corresponding to visual magnitude M »~tT. The method of 
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measuring radiant coordinates has been described adequately in the papers cited. 
The equipment operated at the same high sensitivity in 1957 and 1958, and at 
much lower sensitivity in 1953. 

During 1959 the 67 Mc/s equipment was replaced by a new equipment 
operating at 61 Mc/s, the transmitter power of 50 kW remaining unchanged. 
The observations were made during the period of reconstruction with a temporary 
aerial system with only one receiving channel; the beam axis was directed due 
west at an elevation of 29°. Because of the temporary nature of the installation 
the geometry of this aerial system has not been investigated in detail. The 
resolving power of the system, however, is inferior to that of the 67 Me/s array, 
and the detection threshold is M,~-+6. 


Ill. THE PHOENICID SHOWER 
For the sake of completeness the following details of this day-time shower 
are given. The activity extends from July 5 to July 19, with a fairly sharp 
maximum at July 14, © =112°, when the echo rate is 30/hr. Radiant coordinates 
for 1958, which are the most reliable determinations, are 
a= 31-1°+1-05 (©—109-6°), 
d= —47-9°+0-54 (© —109-6°). 
Data for the earlier years will be found in Weiss (1960). 


TV. THE 5-AQUARID SHOWER 
(a) Observations for 1957, 1958, and 1959 

The daily positions of the radiant, and the echo rates, are listed in Table 1. 
The times of detection of echoes at 700 km slant range, from which the radiants 
were determined, are plotted in Figure 1. 

The shower appears to persist until about August 14, but near the emergence 
of the Earth from the stream the range-time plots show considerable scatter, 
presumably associated with a diffuse radiant. This, coupled with the low echo 
rate, precludes accurate measurement of the radiant. In 1957 echoes were 
received from N and § aerial channels on alternate days, and the radiants for 
individual days were determined by interpolation between times of passage 
separated by two days. Since only one channel was available in 1959, the radiant 
could not be measured to the same order of accuracy as in the previous years, 
and no data have been given. 

The times for the 700 km echoes, plotted in Figure 1, were determined by 
maximizing the echo count within theoretical range-time envelopes. The main 
source of error is random grouping amongst both shower and sporadic meteors. 
We see from Figure 1 that except near the fringes of the stream the error in the 
measurement of these times probably does not exceed +5 min; the consequent 


errors in radiant coordinates are: R.A. -+1°, declination +3°. It is also evident 


from Figure 2, by comparison of the 1957 and 1958 data, that the interpolations 
which are necessary to derive radiant coordinates for the individual days of 
1957 do not introduce appreciable errors. The overall uncertainty in an 
individual radiant determination should not exceed +2° in R.A. and +5° in 


declination. 
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The ‘ envelope”? echo rates of Table 1 for 1957 and 1958 are the actual 
numbers of shower echoes detected inside theoretical range-time envelopes for 
the appropriate radiant positions on the corresponding days. No attempt 
has been made to reduce these envelope rates to hourly rates. Theoretical 
envelopes have not been calculated for the 1959 aerial system, and an empirical 


TABLE 1 
RADIANTS AND ECHO RATES FOR THE §-AQUARID SHOWER 1957-1959 


Radiant Coordinates Echo Rates 

Date 1957 1958 Envelope Total 
ot 8 x 8 1957 | 1958 | 1959 | 1957 | 1958 | 1959 

July 20... $32. (i238 17 97 
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envelope has been determined from the echo rates themselves. The “ total ” 
echo rate is the number of shower echoes detected each day over the following 
periods : 1957 and 1958, N aerial, 01-04 hr ; S aerial, 02-05 hr ; 1959, 22-02 hr. 
Small corrections have been applied to both sets of rates to allow for known 
variations in sensitivity over each run, assuming a mass-distribution parameter* 
s~1-6 for the 5-Aquarids and s=2-0 for sporadics. The value of s~1-6 has 


* This parameter is the exponent in the differential law for the distribution of meteor masses 
namely, v_ dm=const. m-%dm. : 
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been determined (in 1959) from a comparison of the echo rates in two receiver 
channels held at different sensitivities. The high echo rates in 1959, relative 
to 1957 and 1958, are due to the increased aperture of the aerial system used in 
1959, and do not necessarily indicate a higher shower activity. 
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Fig. 1.—Times of detection of echoes at 700km slant range, 
8-Aquarids. 1953, m; 1957, x ; 1958, @. 


(b) Radiant 
Radiant}positions on individual days are plotted in Figure 2. The radiant 
ephemeris estimated by the method of least squares is i 


a= 342-2°40°77 (O©—126-0°), 
$= —16-6°-+40-19 (© —126-0°). 


These daily motions in « and 5 are considerably smaller than those, namely, 
0-96° in a, 0-41° in 8, determined by a different method by McIntosh (1934) 
from the New Zealand visual observations. It is evident from Figure 2 that 
McIntosh’s daily motion in 6 is not incompatible with the radio determination. 
This is, however, not the case for the daily motion in «, nor can the discrepancy 
be attributed to the two different methods used in determining the ephemeris. 
“The ephemeris was formed by McIntosh by reducing all radiant coordinates to 
apparent ecliptical latitude and the difference in longitude, ~A—A, between the 
apparent radiant A and the apex of the Earth’s way, A, taking the mean position 
of the whole and translating it back to the «,5 coordinate system. 6B and A—A 
have been termed apparent coordinates because corrections for zenith attraction 
and diurnal aberration have not been applied. These corrections are small, 
E 
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however, and barely exceed 0-5°. The data used by McIntosh are plotted in 
Figure 3, and it is seen that there is, for the visual observations, no systematic 
dependence of }—A and of f on the longitude of the apex, A. Similar treatment 
of the radio data, on the other hand, reveals systematic changes in both A—A 
and in @ with the age of the shower. These, when translated back to the a,3 
system, yield the daily motions sketched in Figure 2. The corresponding 
ephemeris, which differs only slightly from the least squares ephemeris, is 


a= 342-3°40-76 (@—126-07), 
§=—16-4°40-13. (© ~126-0°). 


The Adelaide radiant coordinates have not been corrected to epoch 1950, 
to which the longitude of the Sun is referred, because the probable errors far 
exceed the corrections. The corrections required to reduce the radio and visual 
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Fig. 2.—Right Ascension and declination of the §-Aquarid radiant. Determinations for individual 

days: 1953, m; 1957, x ; 1958, @. The full lines are the best fits to the radio data determined 

through ecliptical coordinates in the manner indicated in Section IV (b) of the text. The dashed 
lines are similar fits to visual data for 1926 to 1933 (after McIntosh). 


observations, separated by over 20 years, to common epoch will be appreciable 
but as the reference epoch for the visual observations is not stated, no conreetinns 
have been applied. An error not exceeding 0-5° may thus have to be taken into 
account in comparing mean radiants ; the daily motions, of course, will not be 
affected. 

That these systematic changes in A—A and § are real can be seen from the 
following argument. A constant value of A—A implies that the time of transit 
of the radiant changes by only a few minutes over the duration of the shower : 
any change in the time of passage of the shower radiant through the 5 
collecting sectors is thus produced almost entirely by a shift in the declination of 
the radiant. Now, for radiants whose declinations lie in the range 0°> Ss —20° 
the time of passage for the S aerial is extremely insensitive to the dodina gat 
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(see, e.g. Figure 1 of Weiss 1955), and if A—A is constant the time of passage 
through the 8 aerial collecting sector should change by a few minutes at most. 
In actual fact we find (Fig. 1) for the S aerial an advance of 18 min from July 22 
to August 6, which is identical for each of the three years for which observations 
are available. Since a change of this magnitude is quite outside the experimental 
error in the measurement of times of passage, we may conclude that A—A, at 
least, cannot be constant for the duration of the shower. Further consideration 
of the times of passage for the N aerial then indicates that 6 also is not constant 
for the duration of the shower. 
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Fig. 3.—Comparison in ecliptical coordinates of the New Zealand visual radiant and the Adelaide 
radio radiant for the 8-Aquarids. Corrections for zenith attraction and diurnal aberration have 
not been applied, and coordinates are not reduced to common epoch. 


Although it is possible that the discrepancy between the visual and the radio 
ephemerides could be due to incomplete separation of the minor showers from the 
5-Aquarids in the Adelaide observations, this explanation is unlikely in view of 
the excellent agreement between the mean positions of the radiant as found from 
the two sets of data. 


(c) Radiant Diameter 
It is not possible to measure the radiant diameter because two of the minor 
showers cannot be fully resolved from the §-Aquarid radiant. The large varia- 
tions in radiant diameter attributed to the d-Aquarid stream on the basis of the 
1953 observations (Weiss 1955) are now recognized as largely due to the inter- 
mittent nature of the minor shower activity. 
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(d) Activity 
The shower is first detectable on July 20 and persists until August 14. This 
duration agrees with the visual observations. 

‘ An overall picture of the dependence of shower activity on solar longitude 
has been determined by multiplying the envelope rates for 1953 and 1959 by 
factors of 3-2 and 0-30 respectively, and then combining these with the 1957 and 
1958 data. Echo rates so obtained are plotted in Figure 4. The activity is 
slightly asymmetrical about a flat maximum between © =124° and © =125°. 

Variations in equipment sensitivity have been so large that no conclusions 
can be drawn regarding variations from year to year in the activity of the 
d-Aquarid stream. 


ECHO RATE 
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Fig. 4.—The rate of shower echoes falling inside range-time envelopes for the 
5-Aquarids. Echo rates for 1953 and 1959 have been reduced to the 1957/1958 level 
of activity. 1953, m; 1957, x ; 1958, @; 1959, a. 


V. THE MINOR SHOWERS 

The period from mid July to mid August is one of remarkably high minor 
shower activity in the southern hemisphere. It can be seen from Table 1 that 
after excluding sporadic activity the total echo rates for 1957 and 1958 in 3 hr 
centred on the time of peak $-Aquarid activity exceed the envelope rates by a 
factor usually exceeding 2. A similar excess was noted in 1953 when the 
equipment operated at much lower sensitivity. Some at least of the excess echoes 
are to be attributed to 5-Aquarid meteors which fall outside the theoretical 
envelopes, which are computed for a point radiant, but the majority are associated 
with minor stream activity. 
At the higher sensitivity of 1957 and 1958, three main centres of activity 
could be recognized. The times of occurrence of echoes at 700 km range in 
each of the two aerial channels are plotted in Figure 5. The large scatter in each 
set of points is due to the low echo rates of these minor showers, which rarely 
ever exceed the sporadic rate. Radiants have been determined by drawing 
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straight lnes through these plots. In addition, day-by-day radiants were 
computed from the 1958 data, in which information was available on most days 
(see Table 2). In general, the quality of the individual radiant determinations 
is poor, again because of the low echo rates ; a radiant classified as A is probably 
accurate to 3°, but the quality of the radiant deteriorates rapidly for the B and C 
classifications. The radiant ephemerides, and the 1958 measurements for 
individual days, are plotted in Figure 6. 
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Fig. 5.—Times of detection of echoes at 700 km slant range, minor showers 1957 and 1958. Large 
dots represent determinations of good accuracy, and small dots, determinations of low accuracy. 
I, Capricornids ; II, Pisces Australids ; Til, Cetids. 


The identification of the Cetid radiant is beyond doubt. It corresponds 
closely with a visual radiant at 352°/11°, —18° active from July 19 to August 1, 
reported by McIntosh (see, ¢.2., southern hemisphere finding list compiled by 
Ellyett and Roth (1955)). 

Although the times of passage of the other two centres through the N and 8 
collecting areas are fixed with reasonable precision, there is some uncertainty 
in the way in which they are to be paired for the determination of the radiants. 
The criterion adopted was the relative strengths of the activities in each aerial 
channel on the same days, which indicated that the 00-39 passage through the 
N aerial is linked with the 01-31 passage through the S aerial, and 01-22 (N 
aerial) with 02-31 to 02-05 (S aerial). This grouping also gave the best agreement 
with the visual results, according to which the minor shower activity is confined to 
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TABLE 2 
RADIANTS OF THE MINOR SHOWERS 1958 


Capricornids Pisces Australids Cetids 
Date 

4 Oy Class a 8 Class a - ) Class 

July 20 307 —26 B 322, —l7 C 341 —20 C 

21 306 —22 C 326 —10 C 342 —I16 B 

22 313 —17 B 324 —16 A 343 —17 C 

24 314 —18 Cc 325 —18 C 347 —I15 Cc 

20 313 —23 Cc 327 —18 C 346 —l17 Cc 

26 315 —21 B 330 —18 C 348 —19 A 

30 320 —19 C 332 —13 B 350 —20 B 
Aug. 1 321 —12 C 333 —19 B — — — 

2 321 —20 Cc 335 —14 C 353 —18 C 

3 322 —l7 Cc 333 —19 A — — — 

4 319 —26 C 337 —18 C 355 —19 A 


only two radiants other than the Cetids—see Table 3. Although there are some 
discrepancies between the visual and radio radiants, the visual nomenclature 
has been retained to avoid confusion. 

The strengths of all three showers vary irregularly from day to day. The 
Pisces Australids are the most active of the three; and the Cetids, which are 
resolved from the 6-Aquarid shower only with difficulty, are the weakest. The 
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Fig. 6.—Ephemerides for the minor streams of July and August. The three straight lines, which 
are averages for 1957 and 1958, indicate the mean positions of the respective radiants in Right 
Ascension and declination as functions of the longitude of the Sun, whose value is marked along 
the lines. Individual radiants shown are for 1958; each is joined to the point on the mean 
radiant path corresponding to the longitude of the Sun at which it was determined. 


onset of all three showers occurs prior to July 20, the date from which the detailed 
analysis was made; the Capricornids and the Pisces Australids persist beyond 
the middle of August, but the Cetid activity was not detected after August 8. 
The comparatively small excess of the total rates (in 4 hr) over the envelope 
rates for 1959 is thought to be misleading and not a true measure of the minor 
shower activity. As already mentioned, the resolving power of the 1959 aerial 


METEOR SHOWER ACTIVITY IN JULY AND AUGUST d31 


system is lower than that used in previous years and the activity from only one 
of the three recognized minor radiants can be resolved completely from the main 
$-Aquarid activity. The envelopes used in the analysis of the 1959 data were 
determined empirically from the observations, and it is almost certain that a 
considerable fraction of the Cetid and Pisces Australid activity has been included 


TABLE 3 
COMPARISON OF RADIO AND VISUAL RADIANTS 


Adelaide Radio Visualt 
Shower ee : ; 
ties at Radiant Tata ae Radiant 
Activity Piciten 
ctivity 5 Activity zs 5 
Capricornids | July 20*—Aug. 18*| 311/337 —20 July 10-Aug. 5 | 300/325 | —10/—19 
Pisces July 20*—Aug. 18* | 323/346 | —16/—23 July 14-Aug. 22 | 330/339 —30 
Australids J July 26-Aug. 8 | 337/350 | —30/—33 

6-Aquarids July 20 —Aug. 14 | 335/353 | —14/—19| July 22-Aug. 15 | 331/351 0/—17 
Cetids .. | July 20*-Aug. 8 | 342/0 —18 July 19-Aug. 1 | 352/11 —18 


* Shower may commence before, or terminate after, these dates. 
+ From Ellyett and Roth (1955). 


in the envelope activity for the 3-Aquarids. There are, however, in the range- 
time plots some excess echoes whose range changes with time in the manner 
expected of shower activity and which precede the main $-Aquarid activity by 
130 min. These are identified as proceeding from the Capricornid radiant. In 
fact, this time difference of 130 min gives for the Capricornid radiant in 1959 a 
Right Ascension which agrees better with the visual radiant than the 1957 and 


1958 determinations. 
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METEOR HEIGHT DISTRIBUTIONS AND THE FRAGMENTATION 
HYPOTHESIS 


By A. A. WEISS* 
[Manuscript received April 19, 1960] 


Summary 


This paper is concerned with a comparative study of the incidence of fragmentation 
amongst bright and faint meteors, as disclosed by photographic measurements of heights 
and light curves of brighter meteors (visual magnitudes <-+4), and by measured radio- 
echo height distributions of faint meteors (visual magnitudes ~+6). The discussion 
is based on photometric and height data for sporadic and shower meteors obtained by 
Harvard and Canadian observers, and on radio-echo height distributions measured 
at Jodrell Bank. Mean photographic light curves are shorter than predicted by evapora- 
tion theory, and there is a large residual scatter in the heights of maximum brightness 
after reduction to standard meteor velocity and brightness. ‘Trail length is independent 
of reduced height, and also of velocity if this exceeds 25 km/sec ; for slower meteors 
trail lengths shorten rapidly. The radio-echo height distribution as a function of trail 
length is calculated for a simple model which incorporates these features of the ablation 
process for brighter meteors. Comparison with radio height measurements indicates 
that meteors below the photographic threshold show a closer approach to the predictions 
of the evaporation theory than meteors bright enough to record photographically. 


I. INTRODUCTION 

The photographic recording of meteor trails has demonstrated that the 
process of trail formation amongst visual meteors is not correctly described by 
the evaporation theory (Herlofsen 1948), which treats the interaction of a solid 
particle with the atmosphere. Thus, Hawkins and Southworth (1958), in a 
discussion of basic data for sporadic meteors, conclude that the observed lengths 
of meteor trails are considerably shorter than predicted by the evaporation theory 
and appear to depart from the theory in a random fashion. They refer to a 
suggestion by Jacchia (1955) that this discrepancy can be accounted for by the 
phenomenon of fragmentation of the meteoroid. Jacchia, Kopal, and Millman 
(1950) had already advanced this explanation to account for their observation of 


exceptionally short trajectories and abnormally great heights of the 1946 Draconid 
(Giacobinid) meteors. 


These are direct measurements, and the inapplicability of the evaporation 
theory is inescapable. It is more difficult to ascertain the extent to which the 
evaporation theory fails to describe the ionization curves of the fainter meteors 
accessible to radio equipments, because a radio equipment responds to only a 
short segment of the individual meteor trail. This question, however, is of some 
importance, not only in its own right but also because radio measurements of 
the heights of echoing points, of the echo rate, and of echo durations have been 
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applied by Browne ez al. (1956), and other workers, to the determination of the 
distribution of the masses of meteors whose brightnesses extend down to visual 
magnitude M,,~-+7, and Evans (1954, 1955) has derived values of the atmospheric 
pressure and scale height from the height distributions. In all these applications 
the assumption has been made that the ablation of a meteoroid is correctly 
described by the evaporation theory. Unless the nature of the ablation process 
changes significantly over a range of meteor brightness of only two or three 
magnitudes, this assumption would appear to be unwarranted. 


An indirect approach to the measurement of ionization curves of radio 
meteors has been made by Greenhow and Neufeld (1957). These workers 
compared echo amplitudes at two rather ill-defined points on an individual meteor 
trail, and from a statistical analysis of a large number of echo amplitude ratios 
they concluded that the mean ionization curve of these fainter meteors 
(+6<M,<-+8) is shorter than predicted by evaporation theory. 

Apart from second-order effects introduced chiefly by the finite aperture of 
the aerial system and the celestial spread in radiant points of sporadic meteors, 
the width (r.m.s. deviation) of the sporadic radio-echo height distribution depends 
on four main factors. These are (a) the mass distribution of the meteoroids ; 
(b) the atmospheric scale height and its variation with height ; (c) the mean 
ionization curve of the meteor trails; (d) the dependence of the height of 
maximum ionization on the mass of the meteoroid. The mass distribution for 
sporadic meteors is known (Kaiser 19546). If one assumes the evaporation 
theory, (¢c) and (d) are prescribed and hence the width of the radio-echo height 
distribution can be computed as a function of the atmospheric scale height 
Kaiser (1954a, 1954b ; Weiss 1959). It was by comparing with his actual height 
measurements theoretical height distributions computed under this set of 
assumptions that Evans (1954, 1955) was able to study the properties of the 
atmosphere. The atmospheric scale height, however, is sufficiently well known 
to justify an alternative treatment of the height measurements, which is rendered. 
more appropriate by the doubts which have arisen over the validity of the 
evaporation theory. In the absence of a physical theory of fragmentation the 
mean ionization curve (as regards both length and profile) and the dependence 
of the height of the trail on meteoroid mass may be regarded as free parameters. 
A comparison of the measured width of the radio-echo height distribution with 
the width computed with the above-mentioned trail characteristics (c) and (d) 
as free parameters should then lead to an estimate of the extent of the fragmenta- 
tion amongst the sample of meteors measured. 

The necessity for retention of two free parameters, and the possibility of 
random departures of individual meteors from the mean behaviour, are complica- 
tions which effectively preclude the direct application of this procedure to the 
radio meteors. The random variations are important simply because a radio 
equipment responds only to a short segment of the meteor trail. A measurement 
of the height of a single echoing point tells us nothing about the height of the 
point of maximum brightness of the trail, the shape of the ionization curve, or 
the maximum brightness of the trail. This compels us to analyse the radio height 
data on a statistical basis, with the possibility of erroneous conclusions unless the 
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causes of random variations and errors are scrutinized. Random departures 
of individual meteors from the mean height of maximum ionization may well be 
decisive ; random fluctuations in the shape of the ionization curve are of much 
less consequence. 

For these reasons, it would obviously be desirable to have prior knowledge 
of the behaviour of these characteristics amongst the brighter photographic 
meteors, where the incidence of fragmentation is beyond dispute. Neither the 
Harvard nor the Canadian workers carried the analysis of their photographic 
data to this point; they went little further than establishing that, for their 
samples of meteors, the trails were shorter than predicted by the evaporation 
theory, occurred at abnormally great heights, and were subject to random 
fluctuations. As already mentioned, Evans’s approach to the radio-echo data 
was to assume the validity of the evaporation theory and then to proceed to the 
evaluation of the atmospheric parameters ; his numerical values were qualified 
by the remark that if fragmentation were widespread it seemed surprising that 
scale height results which were at all reasonable should be obtained. 


The present paper is a comparative study of the incidence of fragmentation 
amongst bright and faint meteors, based on a more detailed analysis and rediscus- 
sion of the observational material obtained by Jacchia, Kopal, and Millman 
(1950), Evans (1954, 1955), and Hawkins and Southworth (1958). The photo- 
graphic data are first analysed in detail, as regards light curves, trail lengths and 
heights, with particular emphasis on scatter about the mean behaviour. This is 
followed by the theoretical derivation of the height distribution which would be 
observed, by radio techniques, for a model group of meteors which exhibited these 
characteristics, including the scatter. This model height distribution is then 
compared with the actual observations of radio-echo height distributions of faint 
meteors. This comparison suggests that fragmentation is less marked amongst 
the fainter meteors. 


II. PHENOMENOLOGICAL DESCRIPTION OF THE ABLATION PROCESS FOR 
PHOTOGRAPHIC METEORS 

In this section we derive a tentative phenomenological description of the 
observed lengths, light curves, and heights of photographed meteor trails. The 
analysis is based on photometric and height data for 286 sporadic and 74 shower 
meteors, with brightnesses extending down to a little below M,=-+-4, as tabulated 
by Hawkins and Southworth (1958); and for 150 Draconid shower meteors, 
with brightnesses down to M,~+1, measured by Jacchia, Kopal, and Millman 
(1950). These samples are taken to be representative of the brighter meteors. 


It should be emphasized at this juncture that this description may not be 
suitable as the basis for a theory of fragmentation of meteoroids. Rather, 
the analysis is intended to enumerate those characteristics of photographed 
meteor trails which are relevant to the calculation of the corresponding radio- 
echo height distribution. These characteristics are not necessarily those which 
one would choose to examine preliminary to a discussion of the physical causes 
and consequences of fragmentation of the meteoroid. 
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(a) Mean Light and Ionization Curves 
The data given are insufficient to define light curves for individual meteors, 
but by reducing all meteors to the same arbitrary height at maximum brightness 
a mean light curve may be determined. Photographic magnitude M, is converted 
into intensity I by the simple relation 


M,=constant—2-5logl.  .............. (1) 


Mean light intensity curves so obtained are sketched in Figure 1; the profiles 
are normalized so that plate limit corresponds to I/Imax.=1. The light curve 
for the sporadic meteors is designated BB, and that for the Draconid shower 
meteors, CO. For the sporadic meteors the level of maximum brightness above 
plate limit has been tabulated by Hawkins and Southworth. This quantity was 
not given for the Draconids, but since the observations extended over only 5 hr, 
plate threshold has been assumed constant and some 4 magnitude below the 
faintest measured meteor. 


1-0 


DIRECTION OF FLIGHT——> 


O-8 


L/Imax,OR @/Omax 


HEIGHT (KM) 


Fig. 1.—Ionization and light curves for sporadic and shower meteors. Curve AA, 

ionization curve proposed by Greenhow and Neufeld for faint sporadic meteors ; 

BB, mean light curve for photographed sporadic meteors ; CC, mean light curve 

for Draconid meteors; dashed curve, light and ionization curves predicted by 

evaporation theory with atmospheric scale height 6-4km. In each case the height 
scale has been adjusted to height of maximum brightness as zero. 


The mean ionization curve proposed by Greenhow and Neufeld (1957) for 
faint sporadic meteors (+6 <M,< +8) is also reproduced in Figure 1 (curve AA). 
It will be seen that the ionization curve is somewhat longer than the light curve 
for the photographic sporadic meteors, which are about three magnitudes brighter, 

put that both profiles are a good deal shorter than the prediction of the evapora- 


tion theory, namely, 

T/Imax. =4] max. =(9/4)(@/Pmax.)(1 — $¢/Pmax.)” CONG O50. (2) 
(Herlofsen 1948). 9 is the atmospheric density and « the electron (ion) line 
density. The subscript “ max.” refers to the point of maximum electron density 
(and light intensity) on the trail. 


536 A, A. WEISS 


The small sample of shower meteors tabulated by Hawkins and Southworth 
also have light curves which are shorter than predicted by (2). There also seem 
to be systematic differences between the light curves of meteors belonging to 
different showers. In particular, the Taurid meteors show the closest appreach 
to the evaporation theory and the Quadrantids the largest deviations from it. 
However, the samples are far too small for these to be other than hints of 
important differences in the characteristics of meteoroids belonging to different 


showers. 


MAXIMUM PHOTOGRAPHIC BRIGHTNESS, Mz (MAG) 


10 20 30 40 sO 60 7O 
VELOCITY, v (KM/SEC) 


Fig. 2.—The relation between meteor velocity and photographic absolute magnitude 

at maximum brightness, for sporadic meteors. Vertical lines give the r.m.s. devia- 

tions of the mean brightnesses for velocity groups. The full line is the regression 

line for the weighted means, and the dashed curve the relation (11) predicted by 
the evaporation theory, to an arbitrary zero. 


(b) Heights of Maximum Brightness 
On the assumption that corrections to the observed heights of maximum 
brightness, Amax., for velocity and for brightness are independent, these heights 
have been reduced to standard velocity and standard brightness. 
For sporadics the standards are: velocity, 40 km/sec ; absolute brightness 
0-O0mag. The following empirical relations were determined during the course 
of the reduction to reduced height hyax, in kilometres, 


hiner 247 419576 eee meee ee ee (3) 
M,= 2-04— 0909670) eee (4) 
heise gp 2930+ 1-20 enn ee (5) 


Here v is in km/sec, a prime indicates an observed quantity reduced to standard 
velocity and brightness, and the subscript v denotes reduction to standard 
velocity only. M,is the absolute photographic magnitude at maximum intensity. 
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Relation (3) is given by Hawkins and Southworth. The data leading to (4) are 
summarized in Figure 2, in which individual points are the mean brightnesses 
for narrow velocity ranges and the regression line is determined from the weighted 
means. The regression coefficient in (5) was derived in a similar manner from 
the data presented in Figure 3; the individual points are now mean heights for 
narrow ranges of M,,. 


(b) SPORADICS 


=5 


REDUCED BRIGHTNESS, Mp, (MAG) 
REDUCED BRIGHTNESS, Mp, (MAG) 
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REDUCED HEIGHT, Humax. v (KM) 


Fig. 3.—The relation between height and maximum brightness. For sporadic meteors both 

parameters have been reduced to the standard velocity of 40 km/sec ; for the Draconids reduction 

to standard velocity is unnecessary. Horizontal lines are r.m.s. deviations of the means. Full 

lines are regression lines for the weighted means, and the dashed lines are the relation (12) 
predicted by the evaporation theory, to an arbitrary zero. 


For the Draconids (Fig. 3), 
Rona e000 +0-48M,,. Faces sveltetah © euaene 07/0 18 (6) 


The above notation is retained although, in this case, no velocity corrections are 
necessary. As before, the standard brightness is 0-0 mag, with the proviso 
that the different type of plate used here may give rise to a small difference 
between standard brightnesses for sporadics and for the Draconids. 
It is well to stress here that the empirical relations (3)-(6) apply only to 
' the samples of photographic meteors analysed in this paper, but they are expected 
to be typical. After reduction to standard velocity and brightness there remains 
a considerable scatter in the heights of maximum brightness, whose distribution 
is illustrated in Figure 4. There is no dependence of this residual scatter on 
velocity or maximum brightness. Although such a spread is to be expected 
as a consequence of fragmentation, three possible causes of a spurious scatter 
immediately come to mind. They are: (a) imperfect reduction to the standard, 
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arising from uncertainty in the regression lines in Figures 2 and 3; (b) the 
adoption of M, as a basis for reduction to a standard meteoroid. Physically it 
would be more acceptable to reduce all heights to a standard mass. The 
distribution of trail lengths found in Section II (¢) implies that there is no one-to- 
one correspondence between M, and meteoroid mass; (¢) a systematic change 
in the form of the light curve with trail length. (a) is provéd to be negligible 
by the fact that the scatter in Mnax. for the Draconids is unchanged if the theoretical 
relation (12), instead of the empirical relation (6), is used in the reductions. 
(b) is probably responsible for the largest errors in reduction. These errors 
would be minimized if the integrated absolute magnitude were to be used instead 
of M,. However, the Draconid data (Fig. 7) suggest that the reduction in the 
scatter in Rinax.. Which would ensue if this substitution in basic data were 


40 
(a) SPORADICS (b) DRACONIDS 


NUMBER 


80 90 100 110 70 80 30 100 10 
REDUCED HEIGHT, bh’ max.(KM) 


Fig. 4.—The distribution of the heights of maximum brightness, after reduction 

to standard velocity and standard brightness (M,=0). For sporadics the 

standard velocity is 40 km/sec ; for the Draconids reduction to standard velocity is 
unnecessary. 


possible, cannot be very large. The spread in M, for given integrated absolute 
magnitude only slightly exceeds 1 mag, which is to be compared with a total 
range in M, of 8mag. If the spread in M, for sporadic meteors is no larger 
than for the Draconids, most of the residual scatter in hmax, is undoubtedly 
real, and as such is characteristic of the ablation process. Figure 7 of Jacchia 
Kopal, and Millman (1950), in which beginning, maximum, and end fieichts 
for individual Draconid meteors are plotted, leaves one in no doubt as to the 
reality of the scatter in this case. 

The mean reduced height for the Draconids, 95-0 km, is close to the 93-6 km 
found for sporadics. The Draconid geocentric velocity (23 km/sec) is lower 
than the sporadic standard velocity (40 km/sec). We find using (3) that 
Draconid trails are formed, on the average, 9km higher than sporadics of 
comparable brightness. 
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(c) Lengths of Trails 

Observed distributions of trail lengths, defined as beginning height —end 
height,* are plotted in Figure 5. The abscissa is the ratio of observed length to 
the length expected, from the mean photographic light curves BB and CC of 
Figure 1, for the appropriate brightness above plate limit. This quantity will 
be termed the relative length of the trail. The light curve has been assumed to 
be independent of velocity, although this is known not to hold for the brightest 
meteors. Distributions for sporadics and for the Draconids are very similar, if 
the low velocity tail of the former is disregarded. Scatter diagrams reveal that 
relative trail length and reduced height hmax, are only weakly correlated. There 
is a sight tendency for relatively shorter trails to form at greater heights, but for 
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Fig. 5.—The distribution of the observed relative lengths of meteor trails. 

The observed relative length is the ratio of the observed length to the length 

expected, from the mean light curves BB or CC of Figure 1, for the appropriate 

brightness above plate limit. Sporadic meteors with velocities <20 km/sec 
are indicated by cross-hatching. 


practical purposes it is sufficient to regard these two properties of the trail as 
independent. For velocities exceeding 25 km/sec the distribution of relative 
trail length is independent of velocity, but as the velocity decreases below 
25 km/sec the trails shorten rapidly (see Fig. 6). This may imply that parone 
departures from the evaporation theory are a consequence of a long time of 
flight of the meteoroid, and that, on the dustball hypothesis, severe fragmentation 
occurs only after an appreciable time-lag, perhaps associated with thermal 
effects. 

For the Draconids there is no correlation between relative trail length and 
the ratio (h,—Mmax.)/(Mmax.—h,) (0 =beginning, e=end of trail), which suggests 
that there is no systematic relation between relative trail length and the form of 
the light curve. For sporadics, on the other hand, there is a tendency for 
relatively short trails to have the comparatively short, steep rise to maximum 
brightness which is to be expected in cases of extreme fragmentation. This 


* Assuming that the light curve profile is independent of the zenith angle 7 of the radiant, 
the true length of the trail is found by multiplication by sec x. 
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discrepancy between the two groups of meteors may be real or it could be due 
to the small number of segments in the Draconid trails, which in many cases 
makes the position of maximum brightness rather uncertain. 


(d) Integrated Absolute Brightness of Draconid Meteors 
In addition to apparent photographic magnitude at maximum brightness 
(which is readily converted to absolute magnitude M,), Jacchia, Kopal, and 
Millman have tabulated for the Draconids the integrated absolute magnitude, 
M,=—2-5 logfIdt, to an arbitrary zero. Figure 7 is a plot of M, versus M,, 
in which trails longer and shorter than the mean length have been distinguished. 
Since we are dealing with a homogeneous velocity group of meteors, the influence 
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Fig. 6.—Scatter diagram of the relative trail length v. meteor velocity, for sporadic 
meteors. 


of velocity on brightness does not come into consideration. If the integrated 
absolute magnitude M; is taken as a measure of meteoroid mass, this diagram 
implies that, for given mass and velocity, a short trail is accompanied by a larger 
maximum brightness (smaller M,) than a long trail. 

This relation can be given a definite numerical form. Neglecting decelera- 
tion, the above definition of M; implies that for given M,, {Idh—=constant. If 
the shape of the light curve varies with the relative trail length 7 in such a way 
that the difference in height between any two points of equal light intensity is 
proportional to J, it follows that Umax.=constant for given M ; This restriction 
on the shape of the light curve is artificial, but we may still expect that in general 
Umax.~constant for given M,. We see from Figure 5 (b) that for the oreat 
majority of trails 0-4</<2-0, i.e. 1 varies by a factor of 5. Now a ratio of 5 
in Imax. corresponds to a difference of 1-75 mag in M,. This agrees closely with 
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the spread actually found in M, (Fig. 7). We conclude that almost the entire 
spread in M, for given M; is to be attributed to variations in the trail length, 
and that, for given mass and velocity, fIdh is independent of Imax. 

The spread in M, amongst the two classes of meteor, those with J>1 and 
those with /<1, is rather larger than expected, especially amongst the brightest 
meteors. This scatter could well be due to irregularities (flares) in the individual 
light curves, and to minor variations in the shape of the light curves for given 
land M;. It is not, however, the controlling feature in the spread of M, for 
given MV ,. 
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Fig. 7.—The relation between the absolute magnitude at maximum brightness, M,, 

and the integrated absolute magnitude, M; (to an arbitrary zero), Draconid meteors. 

Xx, Trail lengths longer than the mean value for given Mp; @, trail lengths shorter 
than the mean. 


(e) Remarks on the Analysis 
In their original discussions of the photographic data, both Jacchia, Kopal, 
and Millman, and Hawkins and Southworth stressed two aspects of the light 
curves, namely, the shortness of the trails and the random fluctuations which 
occur. The present analysis supplements the earlier conclusions and indicates 
the extent of the departures by individual meteors from the mean behaviour. 
Perhaps the most interesting new feature is the large scatter in the reduced 
height at maximum brightness ; this is so large that it assumes considerable 
importance in the computation of the radio-echo height distributions for frag- 
i oroids. 
Beet fa features of the ablation process, described above, apply equally 
nd to the Draconids, despite the different physical 


dic meteors a ; E 
 aierige wee gs involved. The major limitation arises from the 


structure of the meteoroid 
F 
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lack of light curves of individual meteors, with the result that it has been possible 
to obtain only a rough qualitative estimate of the dependence of the form of the 
light curve on trail length. Nevertheless, it is reasonable to conelude that this 
description should not be modified in any essential detail by more refined 
observational data. 


III. CoMPARISON WITH EVAPORATION THEORY 
The empirical relations (3)-(6) are now compared with the predictions of 
the evaporation theory. If the deceleration of the meteoroid is neglected, 
the equations describing the relevant results of the theory for an isothermal 
atmosphere (see e.g. Hawkins and Southworth for a summary) are 


Cmax (CONG. 8, wets ae ee (7) 


Omax, is the atmospheric density and Imax, the luminous intensity at maximum 
brightness, m and v are the initial mass and velocity of the meteoroid, and 7 is 
the luminous efficiency of an evaporated meteor atom. In an isothermal atmos- 
phere of scale height H, 


Omax. —P9 EXP {(Amax. —ho)/H}, 


and (7) becomes, for constant mass, 
hae CONStanhe ——2 Helin ee ere eee (9) 


Hawkins and Southworth, who adopted the value of H=6-38 km established 
by the Rocket Panel (1952), found that (9) gave good agreement with their 
observational data. 

Writing in (8) t=7t)v" where 7~—0-3 (Weiss 1957) and using Opik’s relation 
between visual magnitude M, and intensity J, 


M ==24 6-2-0 lool. 7 nce ueaiee eee (10) 
we obtain, again for constant mass, 
M,=constant*—6-75 logy. =.........-. (11) 


This is to be compared with (4), with M, replaced by M,, subject to the 
uncertainties associated with colour index. According to Jacchia (1957) colour 
index increases from —1-:8 mag for Mp<—3 to —1-0 mag for M,=+2, but 
the change in colour index is largely if not entirely physiological in origin. In 
any case, the discrepancy between (4) and (11) is so wide (see Fig. 2) that there 
is little doubt that the evaporation theory fails completely to describe the 
observed relation between meteor brightness and velocity. 


* These constants can be evaluated from the evaporation theory if the parameters of the 
atmosphere and the meteoroid are specified. Their absolute values do not concern us here 
With the exception of (12), these constant terms are a function of the mass. Since meteor aos 
and velocity are independent, the relations (9 and (11) will still apply to the mean values of 


hmax, and M,, provided the constant is replaced by an appropriate average over all detectable 
meteor masses. : 
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Eliminating the meteoroid mass from (7) and (8) gives, for constant velocity, 
In Imax, =constant* —3hmax./H, 
which in the notation of (5) and (6) becomes 
hax. »=constant* +1-96 Mp, .........05- (12) 


The observed gradients are less than required by (12), but comparison of this 
relation with Figure 3 shows that, for both sporadic and Draconid meteors, the 
central more populous points are not inconsistent with (12). 


The basic equations (7) and (8) are unchanged if an atmosphere in which 
the scale height increases linearly with height is substituted for the assumed 
isothermal atmosphere (Weiss 1959). The former is a good approximation to 
the actual atmosphere in the meteor region. Hence the only modification 
required to obtain, for the actual atmosphere, the relations corresponding to 
(9), (11), and (12) is to replace H in (9) and (12) by the local scale height. 
Although significant, these alterations are not so large as to invalidate the 
conclusions drawn from the comparisons already made between the observations 
and the predictions of the evaporation theory in an isothermal atmosphere. 


IV. RApIo HEIGHT MEASUREMENTS AND THE FRAGMENTATION OF 
FAINT METEORS 

We now proceed to evaluate the width of the radio-echo height distribution 
which would be measured for the sample of brighter sporadic meteors analysed 
in Section II, and to compare it with actual radio measurements of the heights of 
echoing points of fainter meteors. The observational data used in this com- 
parison are the radio measurements, segregated into velocity groups, made at 
Jodrell Bank by Evans (1954, 1955) on sporadic meteors with M,~-+6. Evans 
has applied corrections for an expected error of +1 km in the individual height 
measurements, and also for the finite range of velocities included in each group. 


The sample of meteors examined by Evans is very different from the samples 
of photographic meteors which have hitherto been considered. They are a good 
deal fainter and are restricted to those meteors for which both velocity and 
height of the echoing point could be measured. Assuming the theory of evapora- 
tion in an isothermal atmosphere, Evans deduced values of the scale height 
increasing from 6-2 km at a height of 85 km to 7-8 km at a height of 107 km, and 
atmospheric densities which agreed well with the density profile from rocket and 
photographic measurements of meteors obtained in New Mexico. 

Since in the actual atmosphere the scale height increases with height across 
‘the meteor region, the scale heights deduced by Evans are too large. This comes 
about because the scale height determined for each velocity group of meteors is 
really an average taken over the whole spread of the heights of the echoing 
points of the meteors within the group. ‘The scale height gradient, which is not 
considered in the isothermal theory, is thus a cause of a small additional spread 
in the height distributions. 


* See previous footnote. 
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The theory of the radio-echo height distribution in a non-isothermal 
atmosphere, in which the evaporation theory and a linear variation of scale 
height with height in the atmosphere have been assumed, has been given by the 
author (Weiss 1959). With the aid of this theory, a reinterpretation of Evans’s 
measurements leads to the scale height 


H 2671-0 -06(,—95 0) ki 5s eee (13) 


for the range of heights h from 85 to 105km. This agrees well with the 
ARDC 1956 Model Atmosphere, which in the height range 90-105 km is repre- 
sented approximately by the linear relation 


H—6'4-0-09(h—95-0) kom... ee ee (14) 


The agreement is probably better than this, since in deriving (13) no allowance 
has been made for loss of echoes due to diffusion of the meteor trail. Diffusion 
will limit the number of echoes measurable in the higher height ranges, thereby 
reducing the measured r.m.s. deviations and mean heights for the higher velocity 
groups. Both mean height and height gradient in (13) are therefore lower limits. 


TABLE 1 
R.M.S. DEVIATIONS OF THE RADIO-ECHO HEIGHT DISTRIBUTIONS 
R.M.S. Dependence of Y 
Deviation Pmax, OF 

Meteoroid Mass 0-0* 0-2 1-0 2-0 cot 

Elementary, 32, .. mt 1-01 0-72 0-54 0-48 0-33 
m® 1-17 0-65 0-31 0-19 0-00 

Complete, da oe m* 1-10 1-07 0-96 0-92 0-86 
m°® — 0-93 0-78 0:74 0-71 


* Evaporation theory. 
} Trails of zero length. 


Were it not for the scatter in the reduced height hinax, for the photographic 
meteors, there would be no reason to doubt that the samples of fainter meteors 
examined by Evans obey the evaporation theory. This scatter is, however, 
So large (r.m.s. deviation +4-2km) that if it is present amongst the samples 
measured by Evans it cannot fail to make a significant contribution to the 
observed r.m.s. deviations. The conclusion that the ablation of the fainter 
radio meteors is correctly described by the evaporation theory therefore cannot 
be accepted without further scrutiny. 

Details of the calculation of the r.m.s. deviation of the radio-echo height 
distribution as a function of the length of the ionized trail will be found in 
Appendix I. The definition of trail length already given in Section II, namely, 
beginning height —end height, is retained. A parameter y is introduced, in terms 
of which the length of the trail between the two heights where the line density 
a=0 is Hln(1+3/y). The complete r.m.s. deviations $x, with which the 
observations are compared directly, are listed in the last two lines of Table ie 
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for several values of y. dx as tabulated includes the scatter in the reduced 
height hinax.. 

For a trail with y=1-0, which corresponds to the mean photographic light 
curve BB of Figure 1, and an assumed scale height of 6-4 km, we see from 
Table 1 that the calculated r.m.s. deviation of the radio height distribution, 
dz, is 6-1km. This is some 15 per cent. lower than the 7-2 km found by Evans 
(1955). This is for the case where pmax.ccm?; if height is independent of mass, 
the deficiency of the calculated below the measured r.m.s. deviation is increased 
to 30 per cent. Since the photographic data suggest that the real dependence 
of height on mass falls between these two extremes, we may expect that the 
calculated r.m.s. deviation of the radio-echo height distribution corresponding 
to meteors brighter than M,~-+-4 will fall some 20 per cent. below that actually 
measured for meteors with M,~-+6. 

There are, however, two deficiencies in the representation of the actual 
photographic light curves by simplified models which will tend to reduce this 
discrepancy. Firstly, the value of y=1-0 was arrived at by equalizing the 
integrals {(a/o%max.)dd corresponding to the mean light curve BB of Figure 1 and. 
to the adopted model (18). The differences in shape of the two ionization curves 
are not insignificant, and this may result in a slight understatement of the 
calculated r.m.s. deviation. Secondly, the radio height measurements do not 
include meteors with speeds less than 20 km/sec. It is seen from Figure 6 that 
about half of these very slow meteors have abnormally short trails, and conse- 
quently the adopted mean value of y=1-0 may be a little high. Against these 
two factors must be offset any contribution to the residual scatter in the reduced 
height hinax, for photographic meteors which arises from irregularities and random 
variations in the shapes of the light curves; and also any loss of echoes by 
diffusion. 

One other factor requires consideration. The values of dx, given in Table I 
were computed for an isothermal atmosphere. In fact, as we have seen, the 
meteor region is one in which the scale height increases quite rapidly with height. 
If the real gradient is approximated by the linear function (14), we find (Weiss. 
1959) that for ~=0 (evaporation theory) dz, should be increased from 1-01 to 
1:07, ie. by 6 per cent. For y>0, da, will require increases of approximately 
this same magnitude, but the increases in dx will be much less, 1 or 2 per cent. 
at most, because of the large correction terms in (21), all of which are insensitive 


to the scale height gradient. 


Although the precise numerical evaluation of some of these corrections is 
difficult, not one of them is large. Collectively they seem quite inadequate 
to influence the first impression that the width of the radio-echo height distribu- 
tion computed for sporadic meteors brighter than M,~-+4 falls considerably 


‘below that measured for fainter meteors, M,~-+6. Measured radio-echo r.m.s. 


deviations can be brought into agreement with predictions based on the photo- 
graphic model only by assuming that the errors in the measurement of heights: 
of individual echoing points have been grossly underestimated, by a factor of 
at least 4. Since this seems improbable, we may conclude that the ablation 
process which is typical of bright (photographic) meteors is not applicable to 
fainter (radio) meteors. 
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V. DISCUSSION 

The conclusion that the consequences of fragmentation are less marked 
amongst the fainter meteors is supported by other radio evidence. It is 
unfortunate, however, that the only radio data which might confirm the photo- 
graphic evidence for fragmentation amongst bright meteors are limited and 
inconclusive. 

On the reasonable assumption that a grouping of meteors by line densities 
at the reflection points is also, on the average, a grouping by absolute brightness, 
the variation of the height distribution with meteor brightness has been ascer- 
tained for a sample of meteors measured at Adelaide by radio techniques. In 
September 1953 heights and line densities were measured for 482 meteors. These 
have been divided into three groups, with line densities at the reflection points of 
2-9 x104, 1-0-1-9 1012, and >2.%x10!% electrons/em. Mean heights for the 
three groups are respectively 90-6, 91-7, and 92-6km; r.m.s. deviations from 
the means are 5-9, 6-8,and6-9km. The increase of mean height with increasing 
brightness is not inconsistent with the hypothesis of an increasing severity of 
fragmentation, but the increase in the r.m.s. deviation is unexpected. 


Evans (1955) has pointed out that estimates, by radio techniques, of the 
deceleration for meteoroids in the same range of brightness as those to which his 
height distributions apply seem to be little, if any, greater than the theoretically 
expected values for a solid particle. 


Finally, the ionization curve of Greenhow and Neufeld (Fig. 1), although 
shorter than predicted by the evaporation theory, is almost twice as long as the 
mean light curves deduced from the samples of photographic meteors under 
discussion. 


A difference in physical structure between fainter and brighter (and therefore 
more massive) meteoroids is an obvious explanation of the closer approach to 
the predictions of the evaporation theory shown by the former. If this is so, 
small meteoroids must be regarded essentially as solid particles. But this is 
not the only possible interpretation, and indeed it is less attractive than the 
alternative which is now considered. 


If it is postulated that the physical structure of large and small meteoroids 
is identical and that they differ only in the number of particles of which they are 
aggregates, the difference between the ablation processes for the two classes of 
meteoroids may be apparent rather than real. The crumbling of a small 
meteoroid, composed of relatively few fragments, will be completed at a much 
earlier point on the flight path through the atmosphere than the fragmentation 
of a large meteoroid. The time-lag in the onset of fragmentation, which we have 
noticed as a characteristic of very slow meteors, supports this view. 


The fragments into which the meteoroid is divided are individually expected 
to obey the evaporation theory. As is shown in Appendix I, this still implies a 
shortening of the ionization curve, which is the more marked the higher the 
density of the atmosphere at the point on the flight path where the evaporation 
of the fragment commences. Consequently, a small meteoroid, which fragments 
early in its flight, will produce a light and ionization curve conforming closely to 
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the evaporation theory for a solid particle; the large departures from this 
theory shown by a massive meteoroid follow from the longer time which is 
required for fragmentation to be completed. 


This alternative explanation requires that the individual particles composing 
the meteoroid be sufficiently minute. 
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APPENDIX I 


The Radio-Echo Height Distribution as a Function of Trail Length 

In this appendix we derive the characteristics of the radio-echo height 
distribution as a function of the length of the ionized trail. Since light and 
ionization curves are proportional, it would be desirable to carry through these 
calculations using a functional representation of the mean light curves of Figure 1. 
This has been attempted, but the resulting mathematics is intractable. Instead, 
anticipating that the height distribution will be much more sensitive to the 
length than the shape of the ionization curve, the height distributions have been 
calculated for a simple model of a fragmenting meteoroid. Since our object 
is to compute the width of the radio-echo height distribution corresponding to 
the picture of a fragmenting meteoroid derived in Section II, the scatter in the 
reduced height hnax, has been incorporated into our model. 


It can be shown that a particle of mass m, inserted at speed into an isothermal 
atmosphere, at a level where the density i8 9, produces by evaporation a trail 
whose maximum electron line density is 


max. =5(6/uH)mo(1 +7/3)8 COSY. sere eee eee (15) 


This maximum occurs at the level emax, specified by 


pmax. 01 tPol3 =k +13). esses eter ee (16) 
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Here we have written y= /e;. 91 is the density at the height of maximum 
ionization for the same particle entering the Earth’s atmosphere from outside 
(e=0). It is given by 
8/3 \* efi cosy 
oe (2) iia PE en ee ary ee 
for a spherical particle ; if the particle is not spherical, a shape factor of order 
unity must be included in (17). The meaning of the symbols employed is as 
follows : 
On=density of meteor particle, 
v=meteor velocity, 
u=mass of individual meteor atom, 
y=zenith angle of meteor radiant, 
H=scale height of isothermal atmosphere, 
68=ionizing efficiency of an evaporated meteor atom, 
l=latent heat of evaporation of the meteoroid, 
A=heat transfer coefficient. 


The usual assumptions, that the deceleration of the particle is negligible, and that 
v?>121, have been made. These are valid for all but the slowest meteors. 


The ionization curve generated by this particle is 

a=0, e<YP1 
ALB; 

391 

Provided that the particle mass m, is replaced by the meteoroid mass m, expres- 
sions (15)-(18) will describe the simple case of a meteoroid for which ablation is 
negligible until the meteoroid fragments, at the level specified by 99, into a number 
of particles, all of mass mp), which individually obey the evaporation theory. 

The length of the meteor trail, defined as the difference in height between the 
two points where «=0, is H In (1+3/y), independent of the zenith angle y. The 


total number of ions generated by the meteoroid is independent of the length 
of trail, as may be verified by evaluating the integral 


—Hiny 
| edhe 
—Hin(3+y) Ue 


2 
=(8/uH) cos 1m (145 ~ ao) » Ye1<e<(3-+y)ex- 


The ionization curve (18) exhibits a progressive shift of the position of maximum 
ionization towards the end of the trail as the trail shortens. These last two 
properties are characteristic of the observed photographic light curves. 

Consider now the case of a homogeneous velocity group of meteors pro- 
ceeding from a point radiant. Using the ionization curve (18) and following 
through the analysis of Kaiser (1954a), we find for the elementary height distribu- 
tion v*, measured within a narrow angular sector of the echo plane, 


of 3 3(s —1) 3Y 3(s—1) 
Stree = (earl yo Oo Ooo (19) 


#=e/e1>(2/3)8. 


with 
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s is the mass distribution parameter, assuming a differential law n,,dmocm-sdm 
for the distribution of meteoroid masses. The first term of (19) is identical 
with the elementary height distribution according to the evaporation theory. 
The second term is introduced by the finite length of the trail, since for any 
given height there is now a maximum detectable mass as well as a minimum 
detectable mass. 

Alternatively, we may consider the case in which the ionization curve is 
described by (18), but with pmax. constant replacing pmax,cmt. The elementary 
height distribution then becomes 


zag Y\ _ 282-0 
: =jeu(1+3) — soy rs Te, eee (20) 


oa 
Let us now introduce the new height variable 
v,=—In z=(h—h,)/H, 


h, being the height corresponding to the density ;. R.M.S. deviations dx, 
of the height distributions (19) and (20) have been computed for s=2-0, which 
is the accepted value for fainter sporadic meteors (Kaiser 1954). This value 
of s is probably a lower limit, since Hawkins and Upton (1958) found s=2-34 
for the sample of photographic meteors analysed in Section II. These values of 
5x, are listed in Table 1 for several values of y. y=0-2 and y=2-0 correspond 
roughly to the longest and shortest trails measured photographically for meteors 
with velocities greater than 20 km/sec. -This cut-off in velocity has been imposed 
to conform to the lowest velocity group studied by Evans. 


Before these r.m.s. deviations can be compared with the observations, 
corrections are necessary for (a) the spread in heights introduced by the finite 
aperture of the aerial, (b) the distribution of radiants over the celestial sphere, 
and (c) the scatter in the reduced heights (Fig. 4), which is approximately 
Gaussian. The complete r.m.s. deviation dx is found from 


(Snj2= (Say)? +0 4+0,+0, +R “aches cee: (21) 


C, and C, are the corrections for aerial polar diagram in azimuth and in elevation. 
They are given respectively by expressions (16) and (20) of Kaiser (19540). 
In the present case they have been computed for a double Gaussian aerial beam ; 
with s=2-0, C,=—C,=0-06. For the height distribution (20), of course, the 
aerial aperture can contribute no additional spread, and these corrections vanish. 
O,, is the correction from a point radiant to the uniform distribution of radiants 
which is simulated by observation over an extended period. It is given by 
expression (9) of Kaiser (19545) ; again with s=2-0, C,=0-08. The r.m.8. 
deviation of the residual scatter in reduced height (Fig. 4 (a)) is 4-2 km, and with 
H=6:4km, R=0-66. 

Complete r.m.s. deviations 5x are also listed in Table 1. These are the 
values with which the observations are to be compared. The corrections for 
aerial polar diagram and for the radiant distribution are sufficiently small that 


uncertainty in these geometrical factors is of little account. 


ON THE IDENTIFICATION OF EXTRAGALACTIC RADIO SOURCES 
By B. Y. Mi1s* 
[Manuscript received March 22, 1960] 


Summary 

Identifications of radio sources with galaxies and clusters of galaxies have been 
sought systematically in a limited region of the sky. The optical data have been taken 
principally from the National Geographic Society—Palomar Observatory Sky Atlas 
and a catalogue of clusters of galaxies prepared by Abell (1958) from the same Atlas. 
The radio data are taken from results obtained with the Sydney cross-type radio 
telescope, supplemented at times with additional information from a recent Cambridge 
Catalogue (3C). A total of 46 possible identifications with galaxies are listed and 55 
possible identifications with clusters of galaxies, the great majority of which are new. 
Most of these galaxies are double systems, but no other common features could be 
recognized: it seems probable that many galaxies of completely normal appearance 
are very strong radio emitters. The possible nature of the double galaxies is discussed 
briefly and attention is drawn to a corresponding duplicity in the radio brightness 
distribution observed in some strong radio sources. Im many cases the emission from 
clusters appears to be associated with a single galaxy or pair of galaxies in the cluster ; 
evidence for the existence of radiation of intergalactic origin is inconclusive. A selection 
of the possible galaxy identifications has been used to derive a provisional radio luminosity 
function ; it appears probable that there is no distinct class of “ radio galaxies ’’, but 
all radio luminosities appear to be represented, the numbers of radio sources in a given 
volume of space falling rapidly with increasing luminosity. A radio luminosity function 
of the form derived is capable of accounting for a substantial number of the radio sources 
of small size observed at high galactic latitudes ; the remainder could be accommodated 
by the uncertainty in the data. It is found that, at the levels to which source counts 


can currently be taken, cosmological effects are likely to be small, although the most 
distant sources are markedly affected. 


I. INTRODUCTION 

To understand the physical significance of radio astronomical observations 
it is necessary to identify the types of objects in which the radiation originates. 
A number of such objects are known; at galactic distances these include 
hydrogen emission nebulae, supernova remnants, and the whole Galaxy itself ; 
extragalactic objects include other spiral and irregular galaxies radiating in much 
the same way as our own, large clusters of galaxies, and very rare types of galaxies 
radiating many orders of magnitude more than the average. However, the 
great majority of detectable radio sources remain unidentified. 

The galactic radio sources are rare objects concentrated very close to the 
plane ; they may be identified only when there is fortuitously a ‘¢ window ” 
through the general dust clouds in their direction, or else when they are very 
close. A high proportion of identifications is therefore not expected. At high 
galactic latitudes, however, obscuration is low and the lack of identification 


must be related to the types of objects causing the emission ; if these are galaxies 
they must be very distant. 


* Division of Radiophysies, 0.S.I.R.O., University Grounds, Chippendale, N.S.W. ; present 
address : School of Physics, University of Sydney. 
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The need is evident for a systematic comparison of a large and homogeneous 
sample of radio sources with an equally homogeneous sample of galaxies. In the 
present investigation such a comparison has been made and many new possible 
identifications of ‘“‘ radio galaxies’? and clusters found. It has also proved 
feasible to take the first steps in the derivation of a radio luminosity function for 
galaxies which is consistent with the existing observational results. 


II. OBSERVATIONAL DATA 

The radio data have been taken principally from a recent catalogue of sources 
between declinations +10° and —20° observed with the Sydney cross-type radio 
telescope at a wavelength of 3-5 m (Mills, Slee, and Hill 1958). For the optical 
data, the basic source has been the original plates of the National Geographical 
Society—Palomar Observatory Sky Atlas in the same region. These have been 
supplemented by additional radio data at more southerly declinations and a 
recent Cambridge catalogue of radio sources (3C) which has been made available 
in advance of publication ;* also many plates taken with the larger telescopes 
of the Mount Wilson and Palomar Observatories have been examined. 


The basic catalogue of radio sources includes 1159 within the declination 
limits of +10° and —20°. The Cambridge 3C catalogue used for supplementary 
information includes 112 sources in the same area, of which about 85 can be 
recognized in the Sydney catalogue. The great majority of the 3C sources 
which could not be recognized were weak sources close to the galactic plane 
where the situation is confused by the existence of a high concentration of 
extended non-thermal and thermal sources, so that different types of instruments 
might be expected to yield different results. When a source was listed in both 
catalogues, it was found that the Right Ascension listed in the Sydney catalogue 
could usually be improved substantially by use of the 3C data. This possibility 
arose because the latter catalogue was prepared using an interferometer which 
gave quite precise values for the time of transit across a fictitious meridian 
inclined at a few degrees to the true meridian, but subject to possible lobe 
ambiguities. Once the correct lobe had been identified (about 30 per cent. 
required lobe shifts to bring them into coincidence with the Sydney position) 
and the Sydney declination adopted (because very large apparent errors in some 
of the 3C declinations caused uncertainty as to how they should be applied), 
the most probable Right Ascension was then determined giving the 3C data 
some three or four times the weight of the Sydney values. For the stronger 
sources it is estimated that the probable error in Right Ascension due to this 
procedure is of the order of 4 min of are, compared with an average probable 
error of about 3 min in declination. These errors would be roughly doubled for 
the weaker sources included in the 3C catalogue. When the 3C data are not 
available, the accuracy in Right Ascension is only slightly better than in 
declination. 

The angular sizes of radio sources having possible identifications were also 
measured at the same wavelength when this was feasible. The instrument 


* This catalogue has now been published with some slight modifications (Edge et al. 1959). 
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used is described elsewhere (Goddard, Watkinson, and Mills 1960) ; briefly, 
it consists of two interferometers of identical sensitivity but different spacings 
between aerials, one of 30 and the other 30000. By measuring the ratio of 
responses on each spacing a crude estimate of the angular size of a radio source 
is possible. The sizes quoted later in this paper were derived on the assumption 
that the brightness distribution across a radio source is circularly symmetric 
and has a Gaussian distribution. On this assumption the angle between half- 
brightness points is given by 
§, 37 -4(In X)* seconds of are, 


where X is the ratio, response at 30A/response at 3000A. When X is large, 
greater than 3 or 4, fine structure in the distribution is likely to dominate and 
lead to a low estimate for the size. As X approaches unity, measurement errors 
become important and limit the size which can be measured to something greater 
than 10-30 sec of arc, depending on the intensity of the source. Because the 
sensitivity of the angular size interferometer is less than that of the radio telescope 
used in preparing the original catalogue, not all the selected radio sources could 
be measured. 

Finally, comparisons were made between the radio sources and a very 
homogeneous collection of data on clusters of galaxies prepared by Abell (1958). 
Abell’s catalogue, being prepared from the plates of the Palomar Sky Survey, is 
based on the same material as the other optical data used. 


III. ANALYSIS OF THE DATA 

A preliminary inspection of the Sky Survey plates showed the futility of 
any attempt to identify all the radio sources. The positional uncertainties of 
the weakest sources are so large that numerous galaxies are contained within 
their limits, and one can only consider seriously coincidences with the very 
bright galaxies which had been investigated earlier (Mills, Slee, and Hill 1958). 
Also, near the galactic plane, the obscuration is so heavy and the star density 
so high that the search appeared hopeless. Accordingly, only a limited sample 
of sources was inspected ; this sample was chosen from the Sydney catalogue as 
follows. 

(i) All Class IT sources (i.e. sources having | b |>124°) with flux densities 
2x10-*° Wm? (c/s)-! or higher. The catalogue is reasonably complete 
in this category. 

(ii) Fainter sources, or sources with lower latitude, recognized in the 30 
catalogue. 

(iii) All sources, not included above, in two areas at high latitudes selected 
because of their low radio background emission and freedom from any 
obvious clustering of bright galaxies. These areas are defined by the 
Sky Survey plates centred on 094 36™ to 112 12m, —18° to 0°; and 
235 12™ to 235 36m, —18° to 0°. 

A total of about 400 positions was examined, made up of about 230 in 

category (i), 20 in category (ii), and the remaining 150 in category (iii). 

Initially, the source positions were plotted approximately on the Sky Survey 

prints using offsets from catalogued stars. The area within three probable errors 
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of the position was searched for the presence of single galaxies brighter than 
m,~16, or double systems brighter than about magnitude 19. All such 
coincidences were then investigated more closely on the original plates. Here 
positions were plotted to an accuracy of about 4 min of arc and coincidences 
taken to be significant only when the discordance was less than two probable 
errors (unless the object showed great peculiarity). The fainter sources in 
category (iii) received slightly different treatment as the listed errors were ignored 
and coincidences sought within 10 min of are of their positions: only galaxies 
showing some peculiarity, or double systems were noted. 


A total of 46 coincidences were considered to be worth recording. The 
angular sizes of 38 of these have subsequently been measured. 


In categories (i) and (ii) (containing 250 sources) there were 39 coincidences, 
of which 19 were with single galaxies and 20 with apparently double or multiple 
systems. To check the possibility of chance coincidences, a number of randomly 
selected positions was examined. Out of 200 positions, each with an associated 
area of 50 square minutes (corresponding roughly to the mean positional uncer- 
tainties), 4 coincidences were found with single galaxies of magnitude 16 or 
brighter, and 2 with very faint double systems of about magnitude 19. The 
number of chance coincidences with galaxies brighter than magnitude 16 to be 
expected from the data of Hubble (1936) is 6, in agreement with the above 
result. Thus it appears safe to conclude that the majority of coincidences, 
particularly with the double galaxies, are significant. The angular size data 
have proved very useful in making individual decisions: they also indicate that 
a very high proportion of the double galaxies are likely to be real identifications. 


Among the radio sources in category (iii), seven coincidences were found, 
comprising six double galaxies and one peculiar single system. On the basis of 
the above statistical tests, about half a dozen coincidences might be expected 
with faint double systems due to chance, so that little weight can be given to the 
coincidences. 

A different procedure has been adopted in making comparisons with Abell’s 
catalogue of clusters. In this case the clusters are often very extended objects 
so that precise positions are of little use. Accordingly, all the 1159 radio sources 
in the Sydney catalogue have been compared with the 877 clusters listed in the 
area of the radio catalogue. Coincidences have been noted for positions agreeing 
within 1 min in Right Ascension and 20 min of are in declination, a compromise 
based on the desire to exclude as many chance coincidences as possible without 
eliminating many real identifications. A total of 55 coincidences were noted. 
The chance coincidences were estimated by comparing the radio positions of 
1950 epoch with the cluster positions of 1855 epoch given in Abell’s catalogue. 
These totalled 16 so that the majority of the 55 coincidences are likely to be real 
identifications. 

LV. PosstsLE IDENTIFICATIONS 

The data relating to observed coincidences have been collected in Tables 1, 
2, and 3. Tables 1 and 2 include coincidences observed with single or double 
galaxies during examination of the Sky Survey plates, and Table 3, the 
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TABLE 3 
COINCIDENCES BETWEEN RADIO SOURCES AND THE CLUSTERS OF GALAXIES IN ABELL’S CATALOGUE 
The photographic magnitude given is that of the tenth brightest cluster member 


Radio Source Cluster 

Position Flux Position 

(1950) Det (1950) 
Ref. No. a Ref. mp 

R.A. Dec. Ww Ss ¥ No. TRolé\e Dec. 
h m ie i et h m Ee 

00—14 00 05-68 —19 58% i 2 00 58:8 —19 54 17:3 
00+03 00 14-23 +06 487 12 16 00 14-1 +06 30 17-0 
00—111 00 25-34 —13 10’ 13 36 00 25-1 —13 03 17:6 
00+07 00 30-88 +05 535 25 55 00 31:3 +06 05 17°2 
00—017 00 54-5! —Ol 39? 90 (72) 119+ 00 53-8 —Ol1 3l 15:0 
o1—14 O1 08-24 —14 33° 16 157 01 08:5 —14 40 16:9 
01—05 01 23-51 —Ol1 35? 88 194} 01 23-0 —Ol 45 13-9 
01—113 01 36-94 —17 49° 10 235 01 37-7 —17 40 175) 
02—09 02 29-48 —04 55% 12 362 02 29-1 —05 05 17°7 
02+010 02 55-1! +05 534 51 4007 02 54-9 +06 O1 13-9 
03—11 03 03-53 —12 215 18 415 03 03:4 —12 14 16:3 
03—12 03 05-44 —16 44° v7 416 03 04-9 —16 55 17:7 
04—05 04 09-53 —Ol1 505 15 477 04 09-6 —02 00 17-5 
04—112 04 32-02: —13 265 38 496 04 31:3 —13 21 15:3 
04—016 04 46-83 —09 555 16 513 04 45-8 —O09 48 17-5 
04+0/1 04 51-54 —02 33° 10 520 04 51-6 +02 53 17:4 
044013 04 56-33 +05 208 10 526 04 57-1 +05 24 16:4 
05—14 05 13-0? —15 568 il 538 05 12-9 —15 45 17:4 
0s8—o02 08 03-13 —00 30° 15 6237 08 03:1 —00 49 16-9 
08 +010 08 55:04 +03 368 9-0 732 08 55:2 +03 22 etal 
09—02 09 06-53 —09 38° 17 754 09 06-4 —09 27 15-2 
o9—l14 09 15:71 —I1I1 53? 690 7807 09 16:0 —12 04 16:6 
09 +05 09 41-83 +09 57° 32 862 09 42-1 +09 44 set 
09+06 09 43-0 +02 21° 7 869 09 43-6 +02 36 17-4 
10—021 10 59-82 —00 52% 23 1148 10 59-8 —00 48 17°6 
11+09 11 38-44 +05 438 8-2: 1346 11 38-5 +05 58 16°8 
11—113 11 47-13 —I1 47° iby 1391 1147-3 —12 02 18-0 
12—0/ 12 01:83 —04 36° 12 1453 12 01:0 —04 23 17°8 
1458 12 01-5 —04 48 17:3 
12—08 12 15-93 —04 47’ 16 ais) l/ 12 16-5 —04 45 16:6 
12—1/1 12 34-0 —14 13’ 9-6 1572 12 34-1 —14 26 17:3 
12—113 19 37-35 —15 38* 14+ 1583 12 37-8 —15 42 17-8 
12—019 12 48-02 —0Ol1 36° 14 1620 12 47-1 —0O1 20 Vie 2 
13+09 13 40-34 +02 20° 13 1773 13 39:6 +02 29 15:6 
13—117 13 59-12 —I11 355 13 1836 13 58:9 —ll 24 17-7 
14+01 14 01:68 +09 22? 28 1850 14 01-0 +09 22 17°5 
14—018 14 53:52 —05 444 16 1994 14 53-6 —05 39 17:7 
15+03 15 08:62 +06 08° 24 2029 15 08-4 +05 57 16-0 
15+06 15 14-45 +00 188 16 2050 15 13-7 +00 17 17-1 
15+05 15 14-22 +07 118 140 2052+ 15 14-2 +07 11 15:0 


+ Coincidences with single or double galaxies in these clusters are noted in Tables 1 and 2. 
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TaBLE 3 (Continued) 


ee en ee es GS ae 


Radio Source Cluster 

Position Flux a 

(A800) Densit (2080) 

y Ref... |———_ m 
Ref. No. (10-26 cen os Aa ‘p 
heres 0 awa Salar Priel, ii = 

15+0/1 15 37-48 +06 088 22 2109 15 38-1 +06 11 17-4 
15—012 15 38-15 —Ol 54° 37 (21) 2103 17 37-3 —02 Ol 17-1 
20—111 20 45-03 —18 207 15 2328 20 45-4 —18 00 16-4 
21—120 21 48-73 —I15 547 8-8 2382 21 49-3 —15 53 16-0 
21—121 21 48-93 —19 537 18+ 2384 21 49-5 —19 47 15-9 
22—012 22 29-24 —08 336 15 2448 22 29-1 —08 42 16-0 
22115 22 35-82 —I17 366 17 2462 22 36-4 —17 38 16-2 
22—014 22 36-98 W—04 138 16 2464 22 36-6 —04 14 17-8 
224015 22 57-24 +09 438 15 2512 22 57-0 +09 50 17-1 
23—13 23 07-73 —10 457 7:6 2544 23 07-6 —I11 05 17-2 
23-402 23 08-23 +07 288 22 2551 23 09:0 +07 38 17°5 
23—14 23 09-63 —12 546 11 2549 23 08-7 —13 05 17-0 
23—08 23 63 —02 297 9-8 25717 23 16-0 —02 33 17-6 
23-—112 23 22-62 —12 295 30 25977 23 22-7 —12 24 16-6 
23—116 23, 27-6? —18 478 13 2612 23 28-3 —18 56 Lig 
23—015 23 38-03 —00 086 11 2644+ 23 38-5 .—00 12 16-6 


coincidences with clusters in Abell’s catalogue. Photographic magnitudes 
given in Tables 1 and 2 have been taken, where possible, from Humason, Mayall, 
and Sandage (1956), de Vaucouleurs (1957), or Holmberg (1958). Where the 
galaxy is too faint for inclusion by these authors, estimates have been made 
using a comparison sequence of images of galaxies calibrated in photo-red 
magnitudes and kindly made available by Dr. G. O. Abell. The red magnitudes 
have then been transformed to photographic, using the data of Abell (1958) 
and Holmberg (1958); it is expected that the accuracy is about 1 magnitude. 
In several cases photographs with the 200-in. Hale telescope have been referred 
to ; this unpublished material was kindly made available by Dr. R. L. Minkowski. 
Also, several of the coincidences noted in the tables had been recognized by 
Dr. Minkowski in an earlier unpublished investigation of a portion of the area 
examined here, using essentially the same data. 


In Table 1 are collected those coincidences for which the available angular 
size data are consistent with an identification. Criteria for such consistency are 
by no means well established and, in general, considerable latitude has been 
allowed. From an examination of known and well-established identifications 
it has been concluded that, with the present aerial spacings, a radio galaxy of 
photographic magnitude 19 would have an expected equivalent angular size of 
2 min of are, with an uncertainty of about a factor of 2 or 3. It has also been 
borne in mind that intrinsically strong radio sources are likely to have a greater 
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extent in space than intrinsically weak sources. A number of coincidences 
excluded from Table 1 because of inappropriate sizes are given in Table 2, along 
with those coincidences in which the radio source is too faint for a size measure- 
ment. Since all the faint sources and all those with incompatible sizes are 
included in Table 2, it is possible that not many of the coincidences are significant. 
In all, there are 30 sources in Table 1 and 16 sources in Table 2. In recording 
the difference of radio and photographic magnitudes the estimated radio 
magnitude at 1-9m wavelength is quoted, based on the magnitude system of 
Hanbury Brown and Hazard, together with scale and spectrum adjustments de- 
scribed before (Mills 1959). The relation used is m,.,.=—53-4—2-5 log S3.;+0°-8. 
In this scale a flux density of 10-?° Wm~-? (¢/s)-? at 3-5 m wavelength corresponds 
to a magnitude of 9-9. 


A brief discussion of each coincidence is given in the ‘‘ notes ’’ column and 
some of their common features are examined in Section VI. One coincidence, 
however, that with the source 16+010 (Hercules A) needs detailed discussion. 
This source had earlier been identified by Minkowski (1957) with the highly 
peculiar double galaxy noted in the table. This identification, together with 
others, has been used to derive a collimation error for the 3-5m Cross; the 
position of the source given by Mills, Slee, and Hill (1958) includes the derived 
correction of +48 Right Ascension and agrees quite well with the galaxy position. 
However, the more accurate Right Ascension of the 3C position appears incom- 
patible with an identification agreeing precisely with the uncorrected Sydney 
position. A new determination has been attempted at Sydney, basing the 
collimation correction on the positions of two well-established identifications, 
Virgo A and Hydra A, closely straddling the source Hercules A in declination. 
However, the result was essentially the same, the indicated correction being 
+35; this position is quite compatible either with the 3C position or the galaxy 
position, and doubt concerning the identification arises because of the claimed 
accuracy of the former. However, inspection of the Sky Survey and 200-in. 
telescope plates shows that there are no galaxies in the immediate vicinity bright 
enough to suggest a radio size as great as that measured, except the system noted, 
and even this is somewhat faint. One possibility is that here we have a close blend 
of radio sources, the dominant one being associated with the double system, the 
weaker with one of numerous faint galaxies in a north preceding direction ; 
examination with an instrument of very high resolving power in both coordinates 
seems necessary to clarify the situation. 


Coincidences between radio sources and clusters are listed in Table 3. Of 
the 55 coincidences noted, it was shown that the chance expectation is 16 so that 
the probability of a selected coincidence being significant is not so high as in 
Table 1. In every case the apparent radio emission is considerably more than 
anticipated from the integrated emission of “ normal ” ovalaxies in the cluster. 
Several of these radio sources have, however, also been listed in Tables 1 and 2, 
indicating that, in these cases at least, the radiation is more likely to be originating 
in a single peculiar system within the cluster than in the cluster as a whole ; 
such cases are marked by a dagger in Table 3. It is not possible to decide at 
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present whether the remaining significant coincidences are similar or whether 
one has examples of ‘radio clusters ’’ where radiation from the intergalactic 
material reaches abnormally high values. 


V. THe Rapio LUMINOSITY FUNCTION 

Interpretation of statistical analyses connected with the distribution of 
radio sources requires, as a starting point, some knowledge of the radio luminosity 
function, that is, knowledge of the number of radio sources in a given volume of 
space as a function of their emission. Since there is, at present, no reliable way 
of determining directly the distance of a radio source by radio means, one must 
proceed by identification of the sources and the application of optical criteria 
of distance. In the following analysis identifications with galaxies only are 
considered. 

There are sufficient possible identifications resulting from the present 
investigation to take the first step in determining the radio luminosity function. 
Thus, if it be assumed that this function is constant throughout the volume 
of space accessible to the present observations, then we determine at the outset 
the probability of a galaxy, chosen at random, having a specified ratio of radio to 
optical emission. The result obtained with the present material can only be 
regarded as provisional, to be improved when better determined positions and 
angular sizes become available ; in particular, the operation of selection effects, 
necessarily imposed by the quality of the data, allows the existence of distant 
unidentified sources of high radio emission. 

In the next step, absolute values for radio emission and space density may 
be obtained from optically derived properties of the galaxies. However, the 
latter step is not taken explicitly here because the optical data are at present 
inadequate. Instead, some statistical properties of the radio sources are deduced 
from the probability function above, and compared directly with observation, 
making some appropriate simplifying assumptions where necessary. 

For simplicity, also, the initial assumption will be made that red-shift 
effects are negligible over the range in which identifications have been sought. 
A later analysis in which these are considered indicates that errors will be small 


compared to those arising from the uncertainties in identification and small- 
number statistics. 


(a) Derivation 
We shall for brevity define the ratio of radio to optical emission of a galaxy 
by m=m..—m,. Then we may write P,,-Am as the probability that a galaxy 
chosen at random has m between m—1Am and m+i4Am. Now the area searched 
systematically for identifications is very close to 104 square degrees so that the 
number of galaxies brighter than m, in the chosen area can be obtained from 
Hubble’s (1936) relation, on ignoring the red-shift correction, as 


log N=0-6m,—9-2 +4-0=0:6m,—5-2, 
whence 


log P,,-Am=log a” Att —0 <6 aD chs maaan itis stalk (1) 


where M»*Am is the number of identifications with galaxies brighter than m ( 
having m lying between m—1Am and m+-4Am. 
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Restricting attention to identifications with m,. brighter than 

9-2 (S3.,>2 10-25 Wm-? (¢/s)-!), which have been looked for systematically, 

it is evident that we can have a complete set of identifications only with galaxies 
brighter than 

Mm ,=9-°2—m. 


Substituting this in (1) we have for this set of identifications 


log P,,-Am=log n,,Am+0-6m—0°3. «6... ee eee (2) 


For single galaxies which have been investigated systematically only to magnitude 
16 this expression is correct for m greater than —6-8. For double galaxies 
investigated to m,=19 the limiting value of m is about —9-:8. It is evident 
that, for single galaxies with m<—6-8 and for double galaxies with m<—9:°8, 
N is a constant defined by the limiting photographic magnitude. Thus for 
single galaxies with m<—6°8 


log P,,, Am=log Nm Am+0-6 xX —6°8—0°3 
=log 1," Am—4 EG Eiehcrgl tc DO ee (3) 


and similarly for double galaxies with m<—9°8 
log P,," Am =log M»* Am —6°2. 6 ccs oie ge Hoe Peco ee (4) 


To construct the probability function the data of Table 1 have been used. 
It is not suggested that all coincidences listed in this Table are real identifications, 
or that all identifications in the area are included, but it appears to be an adequate 
sample for making a preliminary analysis. The photographic magnitudes of 
two galaxies in heavily obscured regions (coinciding with the sources 17 —06 
and 17—23) have been corrected by 4 magnitudes and the sources 14—019 
and 22—09 have been omitted because they do not fit into the scheme 
used in deriving equations (2), (3), and (4) ; the former has too low a flux density, 
the latter an associated galaxy which is too faint. In order to obtain sufficient 
numbers in each magnitude interval, Am has been taken as 2 magnitudes ; 
because of this large interval the equations are no longer exactly correct, but the 
uncertainties in the data do not warrant alteration. Likewise, the numbers are 
too few to differentiate between single and double galaxies, and in Table 4, where 
the results are tabulated, they are lumped together. When equations (3) and 
(4) had to be applied the probabilities were computed separately and then 
added. 

The probability function is plotted in Figure 1, where it is evident that, 
between m=+1 and m=—13, the straight line defined by 


log P,=—1-040-42m eee eee ee eee ees (5) 


. gives a reasonable approximation to the observed points. The point marked 
by a cross at m=+1 has been derived from the data on normal galaxies (Mills 
1959) ; about 50 per cent. of spirals representing about 25 per cent. of all galaxies 
had values for m between +0-5 and 41-5, leading to P,,~t#. 
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(b) Applications of the Luminosity Function 
Let us now consider in a general way the flux densities and numbers of the 
radio sources predicted by a law of this form and compare them with the actual 
observations. Initially we will ignore cosmological effects ; later they will be 
taken into account. 
TABLE 4 
DERIVATION OF THE PROBABILITY FUNCTION FROM THE 


OBSERVED COINCIDENCES 
The interval Am is taken as 2 magnitudes 


mt n,,°Am | log P,, 
fo a = 
=| 1 —1-2 
—- 3 4 —1:-8 
— 5 3 —3-1 
—i7 9 —3-8 
— 9 7 —5-1 
—l1 3 | —6-0 
—13 1 | —6°5 


Substituting equation (5) in equation (1), we find the following relationship 
between v,, m, and m, 


log n= —6 -2 +-0-42m +0+6m,, 


or, if we take /,, as the total number of sources over the whole sky in an interval 
Am=2, we have 
log WV ,= —5°3+0-42m+0-6m, 
= —5-3—0-18m-+0-6m,... «2.2... es (6) 


In defining a family of uniformly distributed radio sources, the most 
convenient parameter is the flux density of the first ‘“‘ average ’’ source, S, (Mills 
and Slee 1957). This may be derived for each value of m by putting W,,=1 
in equation (6), and converting m,., to flux density at 3-5 m using the formula 
in the preceding section, thus ; 

Som =antilog (—24-60 —0-12m). 
Thus we have ns 
So,-13 ==9)-OSC10 225i (e/s)-3 
So, === iris Gl Oem a VV ne (Cia) 
So,-9 =3 -0 X10-** Wm (e/s)-!, ete. 


To obtain the flux density S, of the first ‘“ average ’’ source of the whole 


population, we note that the total number of sources with flux densities greater 
than 8 is given by 


N = (So, 41/8)8/2 +(So,-1/8)3/2-+(8o,-3/S)2/2+. . . 
=(§,/8)*/2, 
whence 


3/2 3/2 3/2 
Se =8971+8924 +982 54+. cere 
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If the population is limited to m=—13, which represents the greatest radio 
emission recognized, we may perform the sum and find that 


Det 10> Wm (c/s) =" 


From the direct counting of observed sources Mills and Slee (1957) concluded 
that S,=2-1x10-?3 Wm-~-? (¢c/s)-1. Since it was also found that a substantial 
proportion of the sources were of large apparent size and therefore not repre- 
sented by the present class of identifications, it would appear that, within the 
uncertainty of the data, the agreement is good. However, cosmological effects 
have been ignored and before a definite conclusion can be reached it is necessary 
to examine the results in a more realistic universe. 


fo) — rr 
[ 


= 


a —L 


4 
+1 -1 -3 -5 7 =. =) 1d 


m 


Fig. 1—The probability that a galaxy chosen at random 

has radio emission defined by m, where m is the difference 

between the radio and photographic magnitudes of the 
galaxy (m,..—mp). 


No attempt will be made to decide between different cosmologies, but one 
particular type will be worked out in detail to illustrate then une Ofs aiteet: 
expected. The Hinstein-de Sitter universe is chosen as a model ; it is the simplest 
of the “exploding ”? models and is characterized by the absence of curvature 
and cosmological constant. If the assumption is made that the total number of 
radio sources of any given emission has remained unaltered during the evolution 
of the universe it turns out that the number counts may be expressed’ by. a very 
simple relation, which is derived in Appendix I. Thus, for a class of source of 
fixed emission we have the relation 


41] 4+27 
gaBif-(A) YO es 3) 


0 
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where is the number of sources with flux density greater than S, 4 is the 
total number of such sources in the observable universe, and the spectrum of 
the sources is of the form 

SOG: 

The correction to the normal source count relation is given by the bracketed 
term, M, is derived fundamentally from the measured red-shifts of identified 
galaxies, but for clarity we will assume a particular value of Hubble’s constant, 
namely, H =100 km/s/megaparsec, and calculate the space density of the various 
source classes and their absolute luminosities. The final source counts, however, 
are independent of the value chosen for H. 

It is first assumed that all identified galaxies have the same absolute photo- 
graphic magnitude. Some such assumption is necessary to proceed because of 
the present inadequate data; it is not unreasonable as shown by Dewhirst 
(1959), at least for those galaxies of the highest radio emission which dominate 
the source counts and for which cosmological corrections are required. The 
mean, absolute photographic magnitude of six identified galaxies listed by 
Dewhirst is M,——20-5, also using H=100 km/s/megaparsec. This value will 
be adopted in the present calculations although it may eventually prove to be 
too luminous, because of selection effects. Thus the absolute radio luminosity 
of a source is given by 


L,,=4n(3 -09 x 107")? antilog —0-4(52-6+M,-+m) W (c/s) 
=1:7 x10”? antilog —0-4m W (c/s)-! at 3-5 m wavelength, 


for example, 
Dy73 2-1 107® W (e/s)-* 


It is readily shown, by consideration of the number of sources in a given 
spherical volume, that the space density, Om is given by the relation 


m=O) Teo Li) a « iin ee (9) 
where So, is obtained as in (7). It is also shown in Appendix I that 
MN) =367 p,(cT')>,* 
where T is the reciprocal of Hubble’s constant, whence 
MN o.m =2167°?(eL)>(So.m/Ln)° = 
=1"3 x0 antilog 042m) i ee (10) 


for example, 
AN jis Ont 10% 


In Figure 2 the source counts expected from equations (8) and (10) have been 
plotted for values of m between +1 and —13. The assumption is made that 
y=+1:0 for all sources. Also shown are the combined counts and those obtained 


iS The symbol Pz is used in Appendix I to denote the space density of objects at the present 
time, Le. it is equivalent to op, above. > | 
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from direct observation. The latter have been obtained from the known strong 
sources in the whole sky (open circles) and the catalogue of Mills, Slee, and Hill 
(1958) covering about 3 steradians (solid circles). The standard amrbhe in the 
numbers due to statistical fluctuations are indicated by vertical wings. To 
obtain the best estimate of the true source count, corrections are ane IG at 
low flux densities; these have been calculated at 10-?5 Wm-? (c/s)-! and 
210-25 Wm? (¢/s)-1 using the data of Mills and Slee (1957), plus an additional 
small correction for a type of error discovered subsequently. This error arises 


2 = as ae T =e 


LoG ‘Nl 


10723 10724 10725 1072° 
FLUX DENSITY (W™-2(c/s)~') 


Fig. 2.—Comparison of the observed source counts with those predicted by the 
derived radio luminosity function. The observational points are shown as 
dots or crosses, depending on their origin, having vertical wings indicating the 
standard error due to statistical fluctuation. Two points corrected for the 
instrumental effects which occur at low flux densities are indicated by crosses. 
Expected counts for each class of source are labelled with the corresponding 
value of m; the total expected count which is the sum of these is also shown. 


from the systematic over-estimation of flux densities not very much above the 
sensitivity limit; it is a partially subjective effect, but very similar between 
different observers. It will be discussed in a later publication. The corrected 
points are indicated by crosses, and a straight line drawn through these represents 
the true source counts as closely as possible with the present data. 

It appears that the numbers derived from Table 1 are too low by a factor of 
about 3 to account for all the observed sources. However, substantial numbers 
of the sources appear to be of large size, possibly irregularities in the galactic 
corona or relatively close blends (e.g. two or more radio galaxies in a cluster) ; 
these would not appear in Table 1 and would not contribute to the radio luminosity 
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function. They comprise about 30 per cent. of listed sources having flux density 
near 5 X10-25 Wm~? (c/s)~! and, perhaps, as indicated by preliminary results with 
the Sydney angular size interferometer, an even greater proportion of fainter 
sources. One is therefore left with a discrepancy of a factor of 2 which, in view 
of the uncertainty in the data, is barely significant. However, there are several 
possible additional identifications between strong sources and galaxies of 
apparently normal appearance which are too faint for inclusion in Table 1 although 
the radio sizes are admissible. These could quite well raise the probability 
function at high emissivities sufficiently to counteract any deficiency. They are 
clearly a subject for future investigation both at radio and optical wavelengths. 
A further possibility, that there is a significant number of sources having m< —13, 
does not, at present, appear very likely, since these would drastically alter the 
form of the source counts. 

One of the most striking features of the computed curve is the small 
influence of cosmological effects. There is no appreciable curvature above 
about 10-25 Wm-? (¢c/s)-! where the mean slope is about —1-3, and between 
10-25 and 10-26 Wm-? (c/s)! the slope decreases to only about —1-0, although 
the sources of highest emission show very marked effects in this region. The 
reason, of course, is the dilution by the intrinsically weaker sources, which being 
closer, show little effects and indeed initially tend to increase the numerical value 
of the slope. The difference between various cosmologies would also be masked 
in this way and, since this difference is unlikely to be large in the first place (e.g. 
Hoyle 1959), the simple counting of observed sources is unlikely to be a very 
sensitive test. However, suitable angular size data could prove very valuable 
by allowing the separation of the various classes of source. 

Finally, one may compute the integrated extragalactic emission from the 
class of radio sources defined by Figure 1. It is shown in Appendix I that the 
equivalent brightness temperature of the integrated background in an Hinstein- 
de Sitter universe is given by 

Lin P CT 
2 dk(B +24) Soy ee 


Making the assumptions as before we find 


L,, =3 +1 antilog 0-02m degK 
whence 
m=—13 
T= > T,=20°K. 
pie oak = 
This is a lower limit for two reasons: (a) it is probable that the derived 
radio luminosity function does not include all extragalactic sources, as. shown 
above, and (b) in deriving the luminosity function it was assumed that all galaxies 
identifiable as radio emitters had an absolute magnitude of —20-5. While this 
value is probably appropriate for deriving the cosmological correction to the 
counts of sources of high emission, for which it. was used originally, it is 
inappropriate for, calculating the integrated emission of the “ normal ”’ Spiral 
galaxies which have the greatest weight in determining 7; these would have a 
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smaller emission, a higher space density, and a greater integrated temperature. 
As an example, if we put M,——16-5 (Hubble’s mean value corrected for the new 
value of H) instead of —20-5, the integrated emission is increased by 2 magnitudes 
(equations (9) and (11)) to 125 °K. Shain (1959) has concluded that the most 
probable value of the extragalactic contribution to the observed temperatures 
near the Galactic poles is about one-third of the total, or about 250 °K, but he 
stated that it could be very much less. 


VI. Discussion 

It follows from the results presented in the last section that the Class II 
sources of small angular size are likely to be extragalactic. Analysis suggests 
that these sources can be substantially represented by a class defined by the 
possible identifications in Table 1, possibly with the addition of a class of 
apparently normal galaxies of relatively high radio emission. 

There is less information about the sources of large size. In a few cases 
these may be identified with clusters of galaxies, either representing radiation 
from the intergalactic material in the cluster or from two or more strongly 
emitting galaxies. However, few of the sources coinciding with clusters in 
Table 3 are listed as “‘ extended ”’ or ‘‘ perhaps extended ”’ in the Sydney catalogue 
and, conversely, extended Class II sources do not show a significantly higher 
correlation with clusters. It is therefore concluded that rich clusters of the type 
listed by Abell are not associated with the majority of the extended radio sources. 
The most likely explanation appears to be that these sources are mainly galactic 
irregularities and/or blends of closely spaced “radio galaxies ”, Angular size 
data which are being collected at present should throw more light on this problem 
and it will not be discussed here any further. 

It is perhaps appropriate to consider a possibility which arises because of 
the apparently continuous nature of the radio luminosity function. If radio 
galaxies are formed only with the highest emission by a rare cataclysmic process 
and thereafter the emission slowly decays to a normal value, the number of 
galaxies having a luminosity between L and L+dZ will be inversely proportional 
to dL/dt. Assuming a luminosity function of the form derived in the preceding 
section, it is easily shown that this would require a decay law approximately 
of the form 

L~I,|(1+ At), 


where A is a constant and ¢ the time after the outburst. However, the 
synchrotron mechanism does not easily lend itself to a decay law of this kind. 
For example, if the outburst produces a spectrum of relativistic electrons defined 
by the power law, nd oc(1/H*)dE (where Bw3), the subsequent decay of the 
- luminosity due to radiation losses at a fixed frequency and in a constant magnetic 
field is readily shown to be 

L=L,(1—Bt)*-?, 


where B is a constant related to the frequency chosen and the magnetic field. 

While it appears that a suitably varying magnetic field could produce the decay 

law suggested by observation, in fact no physically plausible variation has been 
H 
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found which will do so. Also, other features are expected in a decaying source, 
namely, a change in the radio spectrum with age (i.e. with absolute radio 
luminosity) and a high frequency cut-off in the radiation which is also a function 
of age. Present observational evidence does not appear to support the existence 
of these features and it seems that the luminosity function is probably not formed 
as the result of such a process. If a radio galaxy is a transient phenomenon, 
the luminosity function would appear to be determined primarily by the peak 
value of the emission rather than the decay law. 

Let us now consider the galaxies listed in Table 1, which are the most likely 
to be real identifications. These galaxies have no obvious common feature, as 
all types of system are represented ; they do, however, separate into two distinct 
classes, single galaxies and double galaxies. Members of these same classes have 
already been recognized in identifying the strongest radio sources; e.g. the 
radio sources Virgo A and Fornax A are identified with single galaxies and 
Cygnus A, Hydra A, Hercules A, and Perseus A with doubles; the source 
Centaurus A is a possible member of the latter class also. 

Identifications with single galaxies are in general uncertain, for there appears 
to be no reliable optical clue which indicates excess radio emission. The galaxy 
NGC4486, identified with Virgo A, has, it is true, a very blue “ jet’ extending 
from the nucleus. However, this is recognizable only because the galaxy is so 
close; similar features would not be visible usually in the single galaxies of 
Table 1 and, in fact, no abnormal features were noted. In Table 2 one galaxy, 
associated with the source 04-+-0/4, displays heavy absorption in the spiral arms, 
but it is not known if this is a significant form of abnormality. One cannot 
hope to derive much useful information from optical studies of the single galaxies 
listed. 

The double galaxies are rather more interesting. If these are real identi- 
fications, as appears likely in the majority of cases, they must be physical rather 
than optical doubles. The question arises whether they are colliding (Baade 
and Minkowski), separating (Ambartsumian), or orbiting. The first two 
possibilities can be separated from the third on a statistical basis by measurement 
of the radial velocities of the two components, and it is hoped that this information 
will be available before long. 

The collision hypothesis does not appear to stand up to a close examination ; 
the numbers of colliding galaxies in the Universe are expected to be quite small, 
in fact quite negligible outside the centres of dense clusters, as shown by Spitzer 
and Baade (1951), and, if the further limitation is imposed that only giant colliding 
galaxies produce markedly excess radio emission (to agree with the absolute 
luminosities of identified radio galaxies), the interpretation. is completely 
inadequate. The ideas of Ambartsumian (1956) are hard to understand, since 
he suggests that some of the double galaxies may be recently formed systems, 
separating and forming two nuclei; the application of these ideas to very old 
galaxies containing a high proportion of Population II stars, as some of the 
identified galaxies appear to be, does not appear legitimate. However, there is 
some observational support for the suggestion in an earlier identification made 
with the system associated with NGO6166. Dewhirst (1959) describes this 
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as containing within a common envelope two or possibly three nuclei in addition 
to the main system with a velocity difference between the brightest components 
’ of about 1500 km/s. 

In contrast to the general implausibility of these suggestions it appears 
that two closely orbiting galaxies are quite likely to build up strong magnetic 
fields as a result of interaction of the surrounding very tenuous gas The presence 
of a normal concentration of relativistic electrons could then lead to an abnormally 
high radio emission. There is no real evidence that this sort of process occurs, 
however, and NGC6166 could certainly not be explained in this way. Since 
none of these possibilities appears plausible in all cases, perhaps a radio source 
can be associated with any of the three. 

Of particular interest to the double galaxy situation are the results on the 
distribution of radio emission across a galaxy. Here again we appear to have 
evidence of duplicity. The Cygnus A source, with its two concentrations on 
either side of a double galaxy, is now well established, and recently Wade (1959) 
has shown that the extended corona of the Centaurus A source also appears to 
consist of two symmetrically situated concentrations with about the same 
apparent spatial separation as Cygnus A. The situation in Centaurus A is 
complicated by the existence of a small central source responsible for about 
25 per cent. of the total emission and apparently coinciding with the band of 
dust crossing NGC5128 (Mills 1953). A more detailed study of this small source 
by Twiss, Carter, and Little (in press) shows that it, too, is a double system. 
The possibility that these are all toroidal sources is not supported by the detailed 
distributions given, and the present evidence does appear to suggest that they are 
real doubles. 

In view of the scanty data it is rather profitless to speculate about the 
processes by which optical and radio double galaxies are formed. Rather, it is 
necessary to collect a great deal more information both at optical and radio 
wavelengths, particularly about the coincidences listed in Table 1. In this way 
it may be possible to separate the basic factors from a mass of probably irrelevant 
detail. 
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APPENDIX I 


Relations in an Einstein-de Sitter Universe 
This relativistic model universe is characterized by absence of curvature and 
cosmological constant. It is well known (e.g. Bondi 1952) that the appropriate 
metric is given by 
ds? =dt? —R?(dr? +-r2d6? +r? sin? 6d62), 
where 
R=cti. 


Denoting properties at the time of emission of the radiation by subscript 1, 
and at the time of reception by subscript 2, we may write down the following 
relations, using conventional notation. 


The present time ¢,—7%, where T=1/H. 
If 2=Anjny, 
1+2=R,/R, =(t,/t,)-#. 
The luminosity distance 
by 
1 =0r Fey=—clg| dr/R =3et,[1 —(t,/t,)*]. 
ty 


The number of sources of space density g, within q, is then given by 


4 
Na =37 P21} =36rp,078[1 —(t,/t,)4]8. ee (Al) 
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Whence the total number of sources in the observable universe is obtained by 
putting ¢,=0, thus 


AN =36rp,653 =3679,0°T?, 0... eee (A2) 
The flux density of a source of luminosity Z,, at a distance q,, is given by 
S=L,/47q7(1 +2)’. 


With a spectrum of the form Loci’, and in the absence of evolutionary changes 
in ZT we may write 


DL, =D y(v9/v1) 
or 


S=L,/4ngi(ltz)?ty. .......----0-- (A3) 
From (A1), (A2), and (A3), we have 


_Lagt (1 LW ott 
632% NF 2 


which, on using the relation p=6,/z.(S,/L)** (equation (9)), becomes 
S=(8)/M-tE—(V |M oP. ow ee ee eee (A4) 


Consider now the sources in unit solid angle at distance q,. These contribute 
a flux dB, given by 


and the brightness temperature 
tz _ Typo? (?= ° dq 
is =3,| B= Sxk |, (b+e)2+7 


2 Dipyd? ' 
~ Ank(5+2y) 


OTS iy erie cae (A5) 


FREE PATH FORMULAE FOR THE COEFFICIENT OF DIFFUSION D 
AND VELOCITY OF DRIFT W OF IONS AND ELECTRONS IN GASES 


By LG... Huxteys 
[Manuscript received May 24, 1960] 


Summary 


In derivations of formulae for D and W for electrons by the method of free paths 
(Huxley 1957a, 1957b) the assumption was made that along a number of successive 
free paths an electron travels at essentially the same speed c. As this assumption is 
not true of ions it would not be legitimate to apply the formulae to the motion of ions 
in gases without further discussion. However, the formulae are in fact valid for ionic 
motion and in what follows they are established in a more general form. 


I. COEFFICIENT OF DIFFUSION D 
Consider a large number p of free paths x, %,. . ., #, all traversed at the 
same speed ¢ by the ion or electron, but not necessarily consecutively. 


If 7, is the mean free path then 
n 
zt 


Let g be the mean speed of the ion (with speed e) relative to the molecules which 
move at random and let N be the number of molecules in unit volume. 

The sum of the times spent in traversing the paths x, is t=pl,/c, that is to 
Say, p=ct/l,. But p=gtNro*, consequently |,=(e/g)Nxo2=I,(c). For electrons 
g=c; l,=1/Nzro*. ois the limiting value of the impact parameter beyond which 
deflections of the ions are unimportant. 

Consider a diffusing group of n electrons or ions. In the absence of an electric 
or magnetic field the rate of increase of the mean of the squares of the distance 
of the particles from a fixed origin is given by the well-known formula 
dr?/dt=6D. 

In time ¢ let the kth particle traverse a succession of free paths x,,,, so that 
its vector position changes from ro, at time t=0 to r, at time t. That is, 

T.=To t+ LY Xgyme 
m=1 
Also, if rp is the position of the centroid of the group, 
nr, = i r, ; —— 
rere ten = eee 0, 
x =oUr, rT, = UD vyet+ear, 5 Ux x ; 
pat) eee as OE Pha met om Te a Xtom) et met Xen): 
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The total increment in the squares of the distances 7, is, therefore, since the 
middie term on the right-hand side is zero in the mean, 


Art= & (r¢—roe)=( ED Xm) (DE Xm) 


k=1 k=1 m=1 k=1 m=1 
2 
oc OS ele d p.m(Liy(m41) COS (94)k5m +xs¢m+2) COS (B2)i,m++ + -+)}; 


where (0;),,, is the angle between the mth free path and the (m-+j)th free path 
of the kth particle. 


Of the total number T of free paths x,,,, let the proportion of those traversed 
at speeds lying between ¢ and (c+de) be fde. 


The contribution of such free paths to 


2 
>> >>, Lim 


k=1 m=1 


is T2l5fde, where, as explained above, 


n=ho)=(5) (weet) and | * fde=1. 
0 


Their contribution to 


2D YD Lp (Fesme1 COS (94) km +L%inm+2 COS (CAs oe) 
k=1 m=1 


is (Tfde)21,8 ; where S is the mean value of the sum in the bracket, namely, 


[ x Lr,m+j COS (Os )esmler 
j=1 
associated with those free paths ,, traversed at speed c. In general the free 
paths 2;,,4; in S are not traversed at speed ¢. 

Thus the contribution to Ar? by those free paths traversed at speed ¢ is 
Arg =(Lfde)2lp(ly +8), and the average increment over a free path of that quantity 
is 21,(l) +8). Since the average time spent over a free path is J,/e, the mean rate 
of increase is 2c(J,+8) from each particle moving at speed c. It follows that the 
rate of increase of the mean of the squares of the distances r, taken over the 
distribution of velocities over the whole group of particles is 


d a i AS eee 
qr elo 8) =6), 
Thus 


where 


Cc 1 y — ioe. 
1=h,+8 ; =(5) (wot) ' and SL 3 maj 08 Genle= Se) 


The evaluation of S is simplified in the case of diffusing electrons because, if the 
free path @,,,, is traversed at speed ¢, the succeeding paths «,,,,,; are also traversed 
at almost the same speed. 
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Thus S=l,(cos 0,-+cos 0.+. . .+cos 0;+. . .), where cos; is the mean 
value of the cosine of the angles 0; between a free path traversed at speed ¢ and 
its jth succession. 

Consider cos 6,; if the angle between the second and third free 
paths is «, then 0,—(0,+«), and cos 6,=cos 0, cos «—sin 0, sin a, and 


cos §,=(cos 0,)(cos «)=(cos 0,)2, the term in the sines disappearing in the mean. 
Proceeding thus, it follows that (cos 0;)=(cos 9,)/. 
Consequently, S=I, cos 0,/(1—cos 9,) and 


D=H(le) with 1=I,/(1—cos 6,) and 4=1/Nxo®. .... (2) 
If the ions and electrons are scattered isotropically, then cos 0;—0 for all values 
of j. 
It may be remarked that S is also the mean value of the corresponding sum 


for the free paths preceding a free path 7,,,,,,; traversed at speed ¢, since the 
reversed sequence is dynamically possible. 


II. VELOCITY OF DRIFT IN A UNIFORM AND CONSTANT ELECTRIC FIELD 

In the absence of an electric field the number of free paths in a large group 
of n, free paths all traversed at speed c, whose lengths lie between # and x+-dz, 
is (with NA,(¢)=1/l,=(Nxo?)g/c) 

dn,=—n,.N A,(c) exp (—N A,(c)x)da = —N A,(c)n, da, 

where n,;=" exp (—NA,*). 
When, however, an electric field E is present the velocity of the ions or electrons 
changes over the course of a free path and at a distance w along a free path the 


velocity has become c-+Ac(#) where Ac(a#) is proportional to H. In what follows 
it will be assumed that Ac(x)/e<1 and terms dependent on EH? are rejected. 


Thus at the end of a free path of length a, the collisional cross section has 
become A,(c)+{dA,(c)/dce}Ac(~). Consequently, 


dn dA 
mrs =— [Ade HN) M1) 
so that the number of free paths that exceed a in length is 
x 
NN exp| —NAye | Ac(a)ae| : 
0 


and the number of free paths whose lengths lie between # and x+da becomes 


dA 
dn, = ber +N To hee)| jexp ( —NA ger Teel Aetna) dx. 
0 
Consequently, 


d. NdA « 
aie | MA +N ae *(Acta) — Nag Actoja) | exp (—N A,x)- da. 
0 } 
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In terms of the mean free path |,=1/N A, equation (3) becomes 


dn NAG AL db OUR abe? 
pee (—Hirte ale Actoyite —Ae(a)) a ONE oa fac (4) 


0 


Consider a group of » ions or electrons moving in a steady state of motion and 
let the number whose speeds le between ¢ and c+de be dn,=nf(c)de. In the 
course of each free path the ion or electron is displaced by the field E by a small 
amount in the direction of E so that the centroid of the group drifts with a 
velocity W. 

In travelling at speed c at an angle 0 to E along a free path of length a, 
an ion is deflected a distance 4(He/m)(ax/c)? sin 6 at right angles to the direction 0 
and is on this account advanced a distance 4(He/m)(x sin 0/c)? in the direction 
of the force Ee. 

In addition it travels the distance « cos 0 as in free flight without E. The 
increment in speed due to E is Ac(x)=(He/m)(a/c) cos 0 in the direction 0, the 
proportion of those free paths that begin in direction 6 and are terminated between 
x and #-+dz is obtained by substituting this value of Ac(«) in equation (4). 

The mean displacement along a single free path of an ion of the group dn, 
taken over all directions § and free paths x is, therefore, 


Rf Het = Vii 1 dly/ x Ee 
al a ghee 2 Ape Sof Ae He 
if I, {am 2) ee cos OT +R aol a 1) sie a 


Sey i peg ge oe 
3m \e 7m ede *m ec 


The average time taken to traverse a single free path at speed ¢ is J,/c, conse- 
quently the mean speed of drift of the centroid of the dn, electrons of the group 
n is 


It remains to consider the contribution to W, of the equivalent mean exten- 
sion of the free paths # that arises when the scattering of ions or electrons in 
encounters with molecules is not isotropic so that cos 8, is not zero. This 
extension was defined in Section I and was denoted by S=S(c). 

In effect, the displacement in the direction of E of an ion traversing a free 
path « between encounters at an angle § with E is increased from # cos 6 to 
(v-+8) cos 0, but this quantity vanishes in the mean taken over all directions 0. 
Nevertheless, the extension S is responsible for an addition to the mean displace- 
ment as given in expression (5). The field E produces a change in direction of 
the velocity c at the end of a free path , of an amount AQ =(He/m) (a/c) sin 0 
and an increase in the speed by an amount Ac(a) =(Ee/m)(x/c) cos 8. Thus S is 
increased to S-+(d8/de)Ac(#), and through the operation of the field E along « 
the ion or electron receives an additional displacement {S-+(dS /de)Ac(x)}A8 -sin 6 
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in the direction of E. The mean value of this displacement for all free paths of 
length # and direction 0 is 


af” [fas g Hew Z sin? 045% = © cos 6 sin 617 exp (— 7) sin 8a6ae 
Slo 


He (28 dS hel, 1 d 

9 ue ic) a es 7) qo(8%)- LPS Siw col Shae wl Wee eer (6) 
The total mean displacement of an electron or ion traversing a free path at 

speed ¢ is the sum of the final terms in expressions (5) and (6), namely, 


joe nee ere zo Held 
3¢ a Ce AP Ot, finn lida ae 


2mm ec de 

where /=/,+S as defined in equation (1). The mean time to traverse a free 
path is (J,/e), consequently the speed of drift of the centroid of those ions or 
electrons moving at speeds ¢ is 


Wiest He af 4 2), 


3m 2 de 


Consequently, the velocity of the centroid of the group as a whole is, 


Ee 1 
Wea a 


d Ee. 1 -devet 

aes (1g8\ ok rang eat feo 

oy (lc?) tw a w(a) ee ae (7) 
where the equivalent collisional cross section is A=1/NI and the average is 


taken with respect to the distribution of speeds ¢ (Huxley 1957a, where the 
derivation of the formula is restricted to the motions of electrons). 


From the expressions for the coefficient of diffusion D and velocity of drift W 
(equations (1) and (7)) it can be seen that the ratio W/D is given by 


vy ale © (2 *) 1) 
=F ot (1c?) /(ie) ( Stree ee (8) 
He 
~ lme? ] 


Consider the expression ¢~*d(lc?)/de in the special case in which the speeds of 
agitation of the ions or electrons are distributed according to the law of Maxwell, 
that is to say, dn/n=(4/a34/7) exp (—ce?/a2)e2de. 


The mean value required is, therefore, 


4 ed c*\ d 4 9 to e2 
——_ aS Al Ee 2 — 
aval aP ( =a) ae Chea a =e (r: sa) eta 
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The formula for W/D in this case reduces to 


W 3/ He 
pees BR TTIAL ATA (9) 


Suppose instead, that the dependence J(c) of J upon ¢ is of the form loce (inverse 
fifth power interaction), then the factor F also reduces to 3/2 and equation (9) 
is again valid. In general, however, equation (9) is not exact since /’ assumes 
values other than 3/2 when the law of distribution is not that of Maxwell or l(c) 
is not proportional to ¢. In general, for electronic motion, 43/2. 
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HETEROGENEOUS NUCLEATION BY AGGREGATES OF PARTICLES 
By R. B. Heap* and K. L. SUTHERLAND* fT 
[Manuscript received May 17, 1960] 


Summary 


The induction of a new phase by fine particles of a nucleating substance is compli- 
cated not only by variation in sizes but also by aggregation. The general effect of 
aggregation is to decrease the total number of particles but to reduce the supercooling 
(or supersaturation) necessary for growth of a new phase. Point contacts appear to 
be more fruitful than do line contacts, although both may be able to cause a substance 
to be more active than a large flat surface of the substrate. If the contact between 
particles is not a specially favoured position for nucleation, the aerosol as a whole appears 
to contain a higher proportion of active particles because of the persistence of initially 
active particles and the loss of less active ones. 


I. INTRODUCTION 

The nucleation of supersaturated solutions and vapours and of supercooled 
liquids by foreign particles depends upon the size and composition of the particles. 
Fletcher (1958) has related the extent of supersaturation or supercooling to the 
radii of spherical nuclei and the contact angle between the phases formed on 
their surfaces. In any practical work involving aerosols, particularly in the 
laboratory, it is necessary to use such concentrations that particles coagulate. 
Thus typical aerosols of silver iodide used in other studies in this laboratory 
contained 10° particles cm-* (approximate mean radius 10-§ em), the concentra- 
tion of which would fall by coagulation and by diffusive loss to the walls of the 
containing vessel to about one-third in an hour. Frequently it was necessary 
to study aerosols over periods longer than this when the effectiveness of the 


nuclei are being altered by the action of chemicals or of light. This paper 
examines theoretically the effect of coagulation. 


U. ForM oF AGGREGATES 

If the particles are spherical and of equal size, the aggregates form simple 
chains or structures with point contacts (Fig. 1 (a)). We shall confine this 
portion of our discussion to dimers. If the particles are crystals then the form 
of the aggregates is more complex. Only one situation will be examined : 
that of two crystals forming a step (Fig. 1 (b)). 

At edges or points of contact between surfaces the growth of a new phase 
will usually be enhanced by capillary condensation because of the favourable 
curvature of the surface of the phase being formed. If, however, the growing 
phase is crystalline then its fit on both of the surfaces forming the aggregate may 
become dominant. Thus in Figure 1 (b) it may be imagined that the substrate 
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for crystal nucleation may be suitable on the horizontal surface but unsuitable 
on the vertical surface. Under these conditions the aggregate may exhibit the 
characteristics only of the more active member of the pair. Such an aggregate 
shall be defined as ‘not in registry’. For the formation of a liquid phase 
the surfaces are always in registry. 


OSe- 2 Se 


(2) (b) 
Fig. 1—Forms of aggregates. (a) Aggregates from spheres, (b) stepped aggregates 
from plates. 


III. PROPERTIES OF AGGREGATES NOT IN REGISTRY 

Consider an aerosol containing p single particles in each of n classes. These 
classes may be determined experimentally by lowering the temperature or 
increasing the supersaturation by measured amounts, the less active particles 
acting only under the greater driving force. At a given time we shall want to 
know how many monomers, dimers, trimers, etc. are present and to classify 
these aggregates on their probable composition. If a polymer contains at 
least one particle from class 1 (the most active), then it shall continue to belong 
to class 1. If, however, the most active particle in the polymer belongs to class 7 
then the aggregate belongs to class 7. Aggregation therefore decreases the total 
number of particles and eliminates the least active classes. 

The kinetics of coagulation under Brownian movement was formulated 
by Smoluchowski (1916, 1917) on the basis that all particles are the same size 
and that all aggregates are regarded as having the same size and diffusion constant 
as the single particles. 

The rate of removal of particles is given by 


PL y/ Ot ==hV paar itt tte es se se oe (1) 
where k, the velocity constant, is defined in terms of the viscosity of the gas, the 


radius of particle, its sphere of influence, and the temperature (see Green and 


- Lane 1957). 
The total number of aggregates and single particles after a time ¢ is 


Steet let Ty) mee ain. Wate tg (2) 


where T is knpt and =v is the sum of single particles (v,), dimers (v.), trimers (vs), 
m-mers (v,,). The number of m-mers is 


VSSnp LPT (L TMT cease wines ee (3) 
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We now require the composition of each aggregate in terms of the probable 
most active monomer in it. Suppose we select m particles at random to form an 
m-mer. Since there are np particles, the number of different choices is 


2) 


The number of choices containing only particles from classes r+-1,7+2,.. ., ” 
is 
np a 

m 

and the number containing only particles from classes r, r+1,. . ., 1, is 
er? —rp eh 

m 

Thus the number containing only particles from groups r, r+1,. . .,”, and at 


least one particle from group r is the difference of these numbers, is 


(ae —rp | i (i oy) 


m m 


Since these and only these choices form m-mers which belong to class 7, 
the fraction of the m-mers which belongs to class r is 


Pyne?) (mem) (me) “ 


In any experiment the total number of particles within any one class is large 
and the number of particles within an aggregate small. Under these conditions 


is yi is 
aes 


so that (4) becomes 


pe —rp +p)™/m! —(np —rp)™/m! 
is (np)™/m! 
(n—r -+1)™—(n —r)™ 
o ph ot ae Ce ee (5) 


Equation (3) gives the number of each polymer at a given time and the total 
number of particles active in a given temperature class 7 is then 


“A= ¥,"y, Vala, te 0g) ge eee ee (6) 


A sample is considered containing five classes and their distribution between 
classes 1 to 5 (expressed as a percentage of the current number of particles) is 
shown in Figure 2. The percentage active at any temperature step for aggregates 
up to 16 particles is shown in Figure 3, which illustrates the steady march of the 
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sample into the most active class. When the total concentration of the sample 
has decreased to 20 per cent. of its original value (1=4, Fig. 2) the first 
temperature class has increased from 20 to 55 per cent. 
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80 
60 


40 


PERCENTAGE 


20 


Fig. 2.—The variation with time (7=knpt) of the total number of 
particles and of their distribution between classes 1-5 (expressed as a 
percentage of the current number of particles). 


100 


80 


60 


40 


CUMULATIVE PERCENTAGE 


20 


GLASS NO« 


Fig. 3.—The percentage active at given temperature ranges for 

aggregates containing m particles. Any particle active at the 

higher temperature (e.g. class 1) will also be active at a lower 
temperature (e.g. class 5). 


TV. PROPERTIES OF SPHERICAL DIMERS IN REGISTRY . 

Let the aerosol consist of particles each of two spheres of radius FR (Fig. 4). 

Let the surface of the condensed phase lying between the spheres be generated 
by a function y=f(x). We require that the surface free energy be a minimum 
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for a given volume V of condensed phase. Let the angle of contact, measured 
through the condensed phase, be 9. Then the difference in free energy between 
the system with and without a volume V of condensed phase is 

AG =VAG, +8 7¥12—SesV1s +SesVes, te (7) 
where the S and y are the surface areas and surface energies between the phases 
1, 2,3 as indicated. AG, is the free energy of transformation from one phase to 


Fig. 4.—Formation of a new phase (2) from the initial (1) on a nucleus (3) 
formed from two particles each of radius R. §=contact angle ; 
tan @,=tangent to the generating circle; tan @,—tangent to the 
generating curve y=f(z). 

a: a minimum at «=0, tan0,>0, tan 0,>0; 

8: a maximum at «=0, tan 0.>0, tan 0,<0; 

y: a maximum at x=0, tan 0.<0, tan 0,<0. 


another. The relation between contact angle and interfacial energies is assumed 
to be 
19 COS Oy ig——Yas: 2 les Ge oe ee ee (8) 


The volume and surface areas are 


% 
ig =2n| y’dx — §ra2(3R —2aX), 
0 


Sis =n” erie (9) 
0 


So3=47 Rip, 
where y, is dy/da. 
Hence we require that AG, given by 


AG=2r | [AG (yy? —20R +2") +27, 0{y(1 +42)! —R cos 0}]da 
0 


ate a (10) 
shall be a minimum. This condition is obtained from 
aK a aK _y 
a5 ce Pe i A ER Tot (11) 


where K is the integrand’of equation (10). 
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al pt a a fia a= =AGylY2=F, «++ (12) 


which is equivalent to 


d 
os ay 9 +y2)}=%y, 
that is, 
(CLS EU EI 10 20 an 1S) Fp onemces en 3 (13) 


where &, the integration constant, is determined from the condition that 


Y,=9 at w=0 or eet 
E=1/y.(1+Cyp). 


The solution of (13) must be written in several forms depending upon the value 
ore G 


Case (i) C>0 (subsaturation) 

The gradient of the generating curve is everywhere positive and the curve has 
a minimum at v=0. Equation (14) requires €>0 since f, y,>0. Integration 
of equation (13) gives 


1h Ee) 
2 peu |. dx, ONO Gaon a *oro OHO (15) 


where ¥,, the value of y for s=0, is found by solving equation (13) with y, equal 
to zero, that is, 


Wome PEON ALES ME (16) 
This equation transforms to 
Rez|(ab)* (1 —22)de | 
(abso | Ve aye—pPy’ Rate eee (17) 
where 
v=—a/R, 2=(ab)*y, 
P=(b/2a-+1)+(B/4a®rbja)t, Po (18) 
f?=(b/2a +1) — (7/42 +0/a)}, 
a=CR, b=ER. 


This may be evaluated in terms of elliptic integrals 
4 
(F »—LUF(k, Jn) —F(k,9)} (Hl, 47) —B(ky@)}) «++ (19) 
where F(k,47), E(k,47) are the complete elliptic integrals of the first and second 


kinds, F(k,¢), H(k,p) are incomplete elliptic integrals, and 


sin 9 =(1/k)(1 —abu?/e?), 
BES (AS fel) ay aoa aed ee (20) 
u=y/Rk. 
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The parameter b is determined from the condition that the curve shall intersect 
the circles (see Fig. 4, curve «) at an angle 0 (the contact angle). If tan 0, is 
the gradient to the elliptic curve and tan 0, that to the circle, then §,=0,—0, 
and 


tan 0,=y,={b?u2/(1 —abu*)? —1}*=tan 0,—0 
_ (1—9,)/u,—tan 0 
TL —o,) tan O/u, 


where (u;, v;) is the point of intersection of the generating curve. Therefore 


b=(e2-+1)8/{u, -au%(e241)}. 2.2.2... (22) 


Case (ii) C=0 (saturation) 
Reverting to equation (12) then 


Pally?) —y A 0; as Fae ee es (23) 
the solution to which is 

a=(1/b) cosh bt, ... VES Pee were eee (24) 
where 

b=(170,) arsink (1 /e); ee. Sh ee (25) 


and « is defined.in equation (21). 


Case (iii) C<0 (supersaturation) 
Put a=—c where c>0. It may be shown that 6 must be greater than 0. 
Three cases arise, shown in Figure 4, and equation (17) becomes 


Daca (1+-2?)dz 
(ac)#v eee A ve (26) 
(x) A minimum at «=0, tan 0,>0 (Fig. 4, curve «). 
The solution is 


be ale iN Fs. 
(2) = (a) PG, 4=) —Flk,e)} Bld) (hg), -. (27) 
where 
e? =(b/20—1) + (b?/4e2—be)}, 
f2=(b|2e—1) —(0"/4e2 —b/o)}, 


. a be ,\3 
es ae ee 
sin @ x(t aa Be ES, 
=| bu; aa] ot (1—v,)/u; tan 0 
(1 +-beu?)? 1+(1—0,) tan 6/u,’ 
joe 
beux(e2 +1)t—y 

When ¢ is close to 2m the differences of the tabulated functions cease to be 
accurate. The differences are then expanded by Taylor’s theorem about 4m 
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(6) A maximum at w«=0 and tan 0,>0 (Fig. 4, curve 8). 
The equation to be solved is 


al (1+2?)dz 
(aio — | rane y SPivMe “gifs hel (aiers, a (29) 
and the solution is 
bent 
(72) V=SP (yo) tHE g)) oe eeceveeeeees (30) 


where the various symbols are given by equation (28). 
(y) A maximum at «=0 and tan 0,<0 (Fig. 4, curve y). 
The solution is given by equations (30) and (28). 


Case (iv) 
A special case arises when 0,=0. The curve becomes a straight line given by 
MENDON Same neice Sasi seai sk Sets, ete Tools « (31) 
Case (v) 


A special case arises when b becomes infinite, i.e. if the intersection of curve 
and circle is given by 1—cu(c?+1)?=0. The curve becomes a circle given by 


DN aga | ch a en aCe (32) 


(a) Calculation of the Form of the Liquid Surface 
The calculations are confined to the substances silver iodide and water. 
The contact angle for AgI/Liquid H,O/H,O (vap.) is 10°. It is assumed, following 
Fletcher, that this angle applies whether the phases are ice/vapour or ice/water, 
although this is very doubtful. It does, however, enable a comparison to be 
made with the curve for spherical particles obtained by Fletcher’s method (1958). 


We are interested in supersaturated conditions so that no calculations were 
made in the under-saturation region. The procedure in calculation was to take 
a point (u,,v;) on the generating circle representing the particle and then to 
determine ¢ by successive approximations from either of equations ( 27) or (30). 
Equations (24), (31), and (32) enable three curves to be calculated explicitly 
and give a guide to choice of ¢ values for the approximative methods. The 
results of calculation are shown in Figure 5. It will be noticed that the c values 
rise. to a@ maximum and then fall: there are two surfaces for each ¢ value 


(—AGy/Y12)- 


(b) The Free Energy of Embryos 
The theory of nucleation of a new phase on a substrate requires that transient 
groups of molecules (embryos) form on the surface. The relatively large free 
surface energies of the smaller embryos cause these to be unstable : the number 
present of each size at a given temperature and supersaturation may then be 
determined from Boltzmann’s law. With increasing size the free energy barrier 
rises to a maximum and once an embryo exceeds a critical size, it grows without 
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Fig. 5.—Shapes of the menisci from saturation (ec=0) through 
increasing supersaturation (c increasing): O0=10°, generating 
eircle of unit radius. 


limit. Equation (10) expresses the contributions of the free energy changes 
due to surface and to phase transformation and may be rewritten as 


v 


AG =2nAd Re | * wav —2nAG,R02(3 —2v,) 


0 
Y% 
Arya | u(1-+u2)idv —4rry,,R0, cos 0 
= — 3nAG,R5v;(3 —v,;) —4ry1.R0, cos 0 
2 
42nAG,R8 } * Pau 
Uy v 
CU ‘i 2 4 
+4ny,gF? | ean ditt eee (33) 
Uo v 


Since w, is known as a function of w only (equation (13)), the solution may 
be obtained. 


The integrals may be written into a combination of the forms 


2dz 24dz 


dz 
| V/ eee) { V{(e—a)(A—f%)} J o/{(P—2)(22 f° 
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The first of these integrals gives an elliptic integral of the first kind ; the second 
integral gives an elliptic integral of the second kind. The third integral is 
reduced by using the relation 


d 
fie ie 2(1-Lk2) S ae+aet| ae, ce (34) 


where Z is the denominator in the integrals above and k is the modulus of the 
elliptic integral (defined in equations (20) and (28) ). 
The final expressions for the free energies are: * 


Case (ii): a=0 
AG = — 3rAG,Rv?(3 —v) —4rry1,R?0 cos 8 
+ {y1.h?/b? —rAG,R?/2b5}(2vb +-sinh 2vb). oe oD) 
Case (iiix): a=—c, minimum at «#(v)=—0 


AG = — 2rAG,Rv?(3 —v) —4r7y,,R?v cos 0 


e?k? 
+2rAG,R? oe <6 


m,n) —F ke} +h1+" AEP ete 7) 20,9) 


sin @ cos @ (1 —k? sin? ¢)? 


+ eyelash, An) Ee || ee RR cr (36) 


Case (iii8): a=—c, maximum at #(v)=0 


AG = — 2rAG,Rv?2(3 —v) —4r1.R*0 cos 0 


ie2 
+2nAGy Rar ee sin 9 cos ¢(1—k? sin? ¢)! 


Faye) +{14 “I tatte) 


b 
+4ryik spe). ose Bh > ono WP Seat aN (37) 
Case (iv): 9,=0 

AG =27r AG, R3v{u2 —}0(3 —v)} +4rry1,.R?{u—cos 0}. «+ --- (38) 

Case (v): b= : iti , 

AG =§nAG,R\ -(? — 3) (; —2 +) —v(3 ats 

il 1 

+4nyiak\ — AC —v)—V COS at. er ace (39) 


These equations may all be written in the form 
AG =AG,R3H+ERAB. eee cece cece (40) 


The values of H and F were tabulated for the values of a and ¢ shown in Figure 5. 


* The subscript ¢ in v; is dropped beyond this point. 
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(c) Calculation of Relation between R and Supersaturation : 
It is usual to calculate the supersaturation at which the rate of nucleation 
becomes appreciable. For a spherical particle the rate, per particle, is 


J =4710°5R? exp (—AG*/ET), Ssans0---. 0s. (41) 


where AQG* is the free energy of the critical embryo. The pre-exponential term 
is only known to within a factor of 100, but this is relatively unimportant as the 
rate increases very rapidly when the free energy of the embryo is in the neighbour- 
hood of AG*. It will be clear from the curves which follow that the critical 
supersaturation may be readily obtained graphically without direct recourse to 
an equation of the type of (41) which requires AG*~60kT to satisfy it (Fletcher 
1960). 


Bak 19%? 


Fig. 6.—The variation in free energy of a liquid embryo with 

size. The ordinate gives the point of intersection of the 

generating curves (see Fig. 4). R=10-*cm; 6=10°; 

temperature 0 °C; the curves are labelled with the values of 
logy P/P,, the value 0-0313 being critical. 


(i) Condensation.—If P is the partial pressure of water vapour in the gas 
and P, that of the saturated vapour at the same temperature, then 


AG es =-( RT) VO) in (Pee oe (42) 


where V, is the molar volume of water, R the gas constant, and 7 the absolute 
temperature. We shall adopt the same values as did Fletcher for Y12=75 -6 erg cm 
at 0°C. Using 0=10°, the variation of AG, at a given radius for a monomer 
particle, is calculated as a function of v,, i.e. as the embryo increases in volume 
the meniscus will move around the sphere. Figure 6 shows the curves obtained 
at four different supersaturations. At log, P/P,—0-04 the embryo grows 
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Spontaneously even from the smallest size. At 0:03 a hump has appeared in 
the curve and at this supersaturation the difference in free energy between the 
minimum (spontaneous capillary condensation) and the maximum is so large 
that the nucleation rate is infinitesimally small. A free energy barrier (the hump 
in the curve) for the required nucleation rate is 2 x10-" ergs per embryo, so that 
for an appreciable rate of nucleation only a plateau will be apparent, on the 
. seale shown. This is found at log,, P/P,—0-0313 for dimers of particle radius 
103° cm: 

Figure 7 compares the variation of P/P, with the radii of the particles 
forming a dimer to that for monomer particles. It is only for particle radii 
between 10-* and 10-’ cm that appreciable differences appear. 

2:4 


10.7 {O70 1055 1074 
R (CM) 
Fig. 7 Fig. 8 


Fig. 7.—The supersaturation, P/P,, at which water condenses on a dimer consisting of particles 
of radius R. Temperature 0°C; contact angle 10°. The broken curve is for a monomer 
of radius R. 

Fig. 8.—The variation in free energy with size of embryos formed by sublimation at different 
supersaturations; R=10~° cm; §=10°; curves are labelled with the values of AT (the 
supercooling), the value 23-4 °C being critical. 


(ii) Sublimation.—The motive free energy is given by 
AG, =—(RT/V,) in (PyfPy), we eee eee (43) 


where V. is the molar volume of ice, P, is the vapour pressure of a flat water 
surface, and P, that of a flat ice surface. Again, as treated by Fletcher (1958), 
we shall neglect effects of disregistry and lattice strain except insofar as these | 
~ are included in the contact angle relating the surface free energies. The vapour 

pressure P, may be related to the supercooling (AZ =273 —T°) necessary for 


nucleation b 
10249 (P,/P) =aAT, eee, .onfej onle; wir ehs. eu heiye) (44) 
where «—0-00972/degC. A value of 100 erg em-2 (Fletcher 1958) is assumed 
for y1.. Figure 8 shows the variation in free energy for the growth of an embryo 
on a dimer consisting of particles of 10~* cm radius. Figure 9 compares the 
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supercooling required for different radii of monomer and dimer. It 1s apparent 
that the disability conferred on a particle by a finite contact angle is removed 
by the formation of dimers. There is, in principle, a considerable gain in the use 
of dimers formed by spheres in the useful range (10~° to 10-5 em) of sizes where 
a number of particles can be produced per unit weight of seeding material. 

(iii) Freezing.—A similar calculation to that for sublimation is made for 
the process of nucleation by a dimer suspended in supercooled water. Here the 
Y12 value for the ice-water interface is taken as 20 erg em~* (Fletcher 1958) 
and the results are shown in Figure 10. Again there appears to be considerable 
advantage in the use of dimers or polymers. 


AMT(DEG C) 
AT (DEG C) 


10-7 imc 10-> 10-4 eu 


165° 1075 1074 
SM R (CM) 
Fig. 9 Fig. 10 


Fig. 9.—The supercooling AJ’ at which ice is formed by sublimation on: (a) a dimer 
consisting of spherical particles of radius R; (b) a monomer (spherical particles) of 
radius R; (c) a ledge of width and height J, and length J; (d) a ledge of width and 
height /, and length 5/7; (curves (a) and (d) are superimposed). Contact angle 10° for 
(a), (6), (c), (d). (e) A monomer (spherical particles) of radius R but contact angle 0°. 
Fig. 10.—The supercooling AJ’ at which ice is formed by freezing on a dimer consisting 
of particles of radius R. Contact angle 10°. The broken curve is for a monomer of 
radius R. 


V. PROPERTIES OF A LEDGE FORMED BY PLATE DIMERS 

If the aerosol of silver iodide is formed by evaporation, the particles are 
usually crystalline hexagonal plates. Suppose that a dimer forms by one plate 
lying partially over another. For simplicity consider a ledge, such as that shown 
in Figure 1 (b), in which the height and width are equal to J and its length is nl. 
For the calculation n will be chosen as 1 and 5. 

It will be assumed that the embryo is contained by a cylindrical surface and 
having two ends cut square. This form of surface will not be that giving the 
minimum free surface energy and consequently the dimers will be more effective 
than we calculate. At undersaturations the embryo lies within the step and if 
the contact angle is greater than 1x, Wylie (1952) has shown that it rests in the 
ledge as a segment of a sphere. At or near saturation the embryo meets the 
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ledge edges A, B and during further growth the contact angle 0 is indeterminate 
until the embryo forms the angle 0 with the surfaces outside the ledges, when 
it may again move over the surface. 


As long as the centres of the circles forming the embryo surface lie below 
AB, the curvature is convex, and the embryo is constrained between AB, the 
free surface energy is given by 


l 


AG =AGyeo\¢ arsin +3? — see} 


(V/2)o 
+2y12) 0 (azsin Ti —I cos i) +o? arsin mee 
—1eet ey tae), iit Hes bay pl. ee Sa (45) 


where «, the length of the embryo, is less than or equal to nl and 0 is the radius 
of the surface. For the cylindrical surface taken, it may be shown that, provided 
the contact angle is less than 34°, the embryo fills the whole length of the ledge, 
ie. ow —nl, even at the maximum curvature of the surface. For surfaces of small 
curvature, like the flat surface joining AB, the contact angle may be as high as 
in before the embryo segments. 


(a) Calculation of the Relation between Size of Ledge and Nucleation 

Calculations were confined to sublimation nuclei using the constants of 
Section IV (ce) (ii) and a value of n=1 and 5. The most rapid process was to 
plot the variations of AG with the volume of the embryo (the term attached to 
AG, in equation (45)). The curves are similar in shape to Figure 8; Figure 9 
shows the variation of supercooling with step height and length. The dimer 
forming a ledge, of height and width 1, is not as effective as a dimer formed from 
spheres of radius R but most of this difference will arise from the choice of a free 
surface energy which is not minimal (see preceding section). The dimer formed 
with step lengths of 5/ are apparently more effective but to form a step of length 
5l the particles must have dimensions of at least 5/ by l by U or more probably 
Bl by 51 by I, i.e. to the equivalents of particles of 1-71R and 2-92R respectively. 
Plotting the curve for the ledge dimer 51 using either of these values of & 
practically superimposes it on the ledge curve for a step l, ie. the activity at a 
given ledge height per unit volume of dimer is constant. It is likely that ledge 
dimers of all sizes are as effective as dimers formed from spheres. It is certain 
that ledges are advantageous when the particles are large, for then the effect 
of a finite contact angle is almost removed. 


VI. CONCLUSIONS 
An aerosol allowed to aggregate will show a decrease in the total number of 
active particles but with an increasing proportion of particles active at higher 
temperatures. Where point contacts or ledges form (and, by inference, corners) 
these are likely to be very active in two senses. They generally raise the 
temperature at which nucleation is possible on a given substrate. This is of 
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considerable importance since the number of natural clouds at temperatures 
close to 0°C exceeds those at the more usual temperature (—10 °C) where 
nucleation with silver iodide is effective. For a given mass of silver iodide it 
might well be more effective to dispense slightly coagulated smokes. If the 
particles in the smokes are tabular there may not be any advantage in coagulation, 
except for very coarse smokes, which are probably uneconomic in application 
anyway. Methods of producing monodisperse aerosols are being studied so 
that the various models may be differentiated. 
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SHORT COMMUNICATIONS 


THERMALLY ACTIVATED FERROMAGNETIC DOMAIN WALL MOTION* 
By F. D. STAcEYt 


The properties of small ferromagnetic domains are influenced strongly by 
thermal agitation if the energies involved in transitions between different states 
of magnetization (rotations and domain wall movements) are only a few times kT 
(k being Boltzmann’s constant and 7 the absolute temperature). The study 
of such small-scale magnetic phenomena, now known as “ micromagnetics ”’ 
(see, for example, Brown 1959a), has developed rapidly in recent years in response 
to commercial interests in fine powder permanent magnets and thin film memory 
cores. The possibility of observing directly thermally induced domain wall 
movements is therefore a matter of some interest. It appears that Olmen and 
Mitchell (1959) have made such an observation, although they did not specifically 
claim to have done so. The purpose of this note is to present an elementary 
theory of thermally activated domain wall motion and to show that the results 
reported by Olmen and Mitchell are in agreement with it. 

The delayed response of magnetic domains to changes in external field is 
termed ‘‘ magnetic viscosity ’’ but some ambiguity has arisen from the two 
different physical processes involved. Domain wall velocities have been observed 
in metals by Sixtus and Tonks (1931, 1932) and in ferrites by Galt (1954). In 
both cases the velocity v in a field H is 


v=G(H —H,), 


H, being a critical field (static coercive force) and @ is a constant of the order 
104 em sec—! oersted-!. In metals the velocity is limited principally by eddy 
currents. Galt explained the velocities observed in ferrites in terms of the 
dissipation of anisotropy energy during ionic reordering in moving walls. Galt 
(1954) and George (1959) have referred to this latter process as magnetic viscosity. 

The expression is also used for much slower changes in magnetic domains, 
such as were observed by Street and Woolley (1949) in ‘“‘ Alnico ”? and subsequently 
recognized as important phenomena in rock magnetism (e.g. by Néel 1955). 
This longer period magnetic viscosity is controlled by thermal activation, which 
must be invoked to explain the slow domain wall motion observed by Olmen 
and Mitchell. 

Consider a domain wall crossing a sequence of potential barriers, which will 
be assumed for simplicity to have equal heights E,. In the presence of a field H 
in the direction of the expanding domain, the energy which must be supplied to 
the wall in order to impel it across a barrier is 

yd eet 4 | (1) 


* Manuscript received March 22, 1960. 
+ Geophysics Department, Australian National University, Canberra. 
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where energy («H) is supplied by the field during penetration of or ‘ climbing 
up’? the barrier. The wall can move freely only if H<0, but if H>0, movement 
can be induced thermally. The probability P of the wall acquiring in unit time 
sufficient thermal energy to cross a barrier is 


P=CO exp (—B/kT). «2. cess eee cveces (2) 


Alternatively, P is the number of occasions per second on which the wall acquires 
thermal energy HZ which could impel it across the barrier. There are several 
calculations of the frequency factor C. Néel (1955) suggested 10° to 10% ies 
Stacey (1959a) obtained C =6 x 101 sec! at room temperature, but this must be 
regarded as an upper limit, as it expresses the probability of sufficient energy 
becoming instantaneously available rather than the probability that the wall 
will move sufficiently rapidly to make use of it (Stacey 1959b). A completely 
different formula for specific application to single domains has been given by 
Brown (1959b). In the following calculation two alternatives are assumed : 
(i) Cis a large but arbitrary constant with negligible dependence upon parameters 
such as temperature. This leads to equations (5A) and (6A). (ii) C is pro- 
portional to 7, as suggested by Stacey (1959a). This leads to the alternative 
equations (5B) and (6B). The difference is not very significant in the present 
case. 

If the average separation of potential barriers in the direction of motion is d, 
then, assuming that the delay of the wall at each barrier is much greater than the 
time taken to move between barriers, the wall velocity v is given by 


o=Pd=Cd exp £=(L,—ab kT}. yaa aanme (3) 
Putting Cd=A, a/k=B, and E,/k=BH,, we obtain a simpler form 
=A oxp (B/PVH—Ao” ween eee (4) 


Equation (4) has the exponential dependence of wall velocity upon applied field 
which is required by the observations of Olmen and Mitchell. It must be noted 
that H —H, is necessarily negative in the range of interest as A is a large constant. 

The temperature range of the observations was too limited to allow detailed 
comparison with (4) but the variation of v with 7 can be used to estimate the 
constant Hy, which is the field required to move the wall at 0 °K, or to move it. 
very rapidly at any temperature. By differentiating equation (4) we obtain 

diin'v)/d?=—B(A Ay? ee ee ee (5A) 
if the frequency factor C is independent of 7, or 

d(ln v)/dT =—B(H—H,)/f?+1/T ........ (5B) 
if CO is proportional to 7. 

Examination of the data of Olmen and Mitchell shows that for their film 
d(In v)/d7=0-09 degC-1 at T=320°K and H=4-Ooersteds and that 
B/T =8-8 oersted-1, so that from (5A) H,=7:-3 oersteds, and from (5B) 
H,=7-1, oersteds. It is not possible to distinguish experimentally between these 
alternatives but, as the coercive force of the film when measured at 60 c/8 was 
7 oersteds and would therefore be slightly higher still at even more rapid switching, 


ae must be regarded as an excellent agreement with the theoretical estimates 
Otel. 
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Having determined Hy, the absolute values of measured velocity may be 
used to estimate A. From Olmen and Mitchell’s data, at 4-0 oersteds switching 
field and 320 °K., spike velocity is 3 x10-2 em/sec. Taking true wall velocity as 
one-seventh of this, as noted by Olmen and Mitchell, we find 


A =Cd=1-7 x10" cm/sec 
if C is independent of 7, or 


ae Od= bo KAO CmiseG® lak es ene ote (6B) 
if C is proportional to 7. 


At this point an inadequacy of the simple theory must be noted. It has 
been assumed that the wall is inflexible and that it only meets one barrier at a 
time, neither of which can be true. Rather we should expect the wall to encounter 
a more or less straight line of elementary barriers and tend to fold round each of 
the barriers, making individual jumps past them. Im this case the wall would 
effectively be divided into a large but unknown number of semi-independent 
sections, each opposed by a single elementary barrier. This allows a physically 
reasonable interpretation of the constant « in equation (3), since we may put 


«=2ail,, 


where J, is the saturation magnetization of the film, a is the average wall area 
of each effectively independent section of wall, and ¢ is the effective half-thickness 
of the elementary potential barriers. However, this does not alter the foregoing 
equations, since we can consider the wall velocity of each section of wall separately 
in terms of the same equations. 

If a direct microscopic observation could be made of the individual jumps 
of a slow-moving domain wall, so that d would be measured, the fundamental 
constant C could be estimated from equations (6). It may be noted as a pre- 
liminary result that a value of C as low as 101° sec does not appear to be con- 
sistent with the above figures, as it gives much too large a value of d. Taking 
the theoretical upper limit, C=6 x10” sec-! (at room temperature) we would 
expect to find d~10 cm. It therefore appears possible that such an observation 
could be made. 


The helpfulness of suggestions by Dr. P. Gaunt of Sheffield University and 
Dr. W. A. Rachinger of the Aeronautical Research Laboratories, Melbourne, is 
eratefully acknowledged. 
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ON THE RELATION BETWEEN LUMINOSITY DISTANCE AND 
DOPPLER SHIFT IN RELATIVISTIC COSMOLOGY* 


By R. VAN DER BORGHTT 


Introduction 

An approximate relation between luminosity distance D and Doppler 
shift 5 in cosmology is usually obtained by a succession of complicated expansions 
in series. 

In the present paper it is shown that an expression for D, even to a higher 
approximation in 5, can be obtained in a much simpler way. 

The relation between luminosity distance D and red-shift $ commonly used 
in relativistic cosmology is (MecVittie 1956)t 


hi +h. 
=P (1+ a °8), eee eo (1) 
where 
h,=R(ty)|R(t)) =Ry/Rp, i.e. the Hubble constant, 
pe Stenaiveh (2) 
and 


R being a function of ¢. 

The exact expression for D is 
dis r 
D=sZ 1T+ar?/4’ O) 6, ©) 10) w ise. we fe fe. te 0) ee) 6) 6 (4) 


where r is a function of & given by the null-geodesic equation 


C . Zi! ihe Ls dr 
iP roma heeace res ths en (3) 


The constant «—-+1 for an elliptic space, 
=0 for a flat space, 
=-—1 for a hyperbolic space. 


The left-hand side of (5) could be written 


of at, | ae 
, Rit) es 


and therefore . 


* Manuscript received March 30, 1960. 
} Visiting Fellow, A.N.U., Mount Stromlo Observatory, Canberra. 


{ MeVirirm, G. C. (1956).—‘‘ General Relativity and Cosmology.” (Chapman & Hall: 
London.) 
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where 
‘ tore 
=C6 Sed De Coca oh ence ear ae OeNe 7 
2 BR (7) 
If «=-+1 it follows from (6) that 
r=2 tan ty, 
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and substituting this in (4) we have the following expression for the luminosity 


distance 
D=(R6/R) sin y. 


It is also easily seen that: 
if «=0, D=(R5/R)U, 
if «=—1, Veh iral le) ain bee nce. 5 hs a eeeke (10) 


Elliptic and Hyperbolic Spaces 
For these spaces 


where 
‘ R, 
F(R) = "2 (°| =i Pe (12) 
sinh pr RE 


the sin corresponding to the elliptic space and sinh to the hyperbolic one. 
Using a Taylor expansion for Ff, we have 


F(R) =F) + (FR) (R—R,) +(im), (RR)? (3) (RR), 


dR), dR? 2 \are é 
oulaenre Oe Re (13) 
where, firstly, F(R,)=0, 
secondly, 
ae cos /, fm dk\ 1. 
dk  cosh\ |, RR)RR 
that is, 
(i) = Pea ee ee (14) 
Ry 0*'0 
thirdly, 
i Ry = (Ro d 1 
See a C di) =1 8 Cc nae a) e gree  s) 
dk? sinh pn RR) RR? cosh pn RR/AR\RR 
But : - 
dave (le) 2 eae eA 6) 
(25) <3 (a) 


Substituting this in (15), 


de pe eye ae ad (o ee eee 
—.=+¢ ( ‘ RR R2R cosh hy RR R2R3 b) 
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and 
es . 
Fel oo et) 2 ee (18) 
dk"), Koko 
Fourthly, 
dt? = 3 COS /, Ro dh —1 
dk? cosh iL RR) R8R? 
: Ry ee 
ea iG | le _2R-*R-1 Rk) > 
sinh Rr RR R 
, sin / (® dR\R?+RR 
= sinh (’| n RR) RR 
to 08 (,[% ak L2os 28 eee 
Al munk RE RR PBR RRA = 
and 
ar 8 ¢ Do etek 3Re 
2) a tesgl-— tle 8) as 
aU EARR SR BRS Reboe de eee ae 


Substituting (14), (18), and (19) in (13) we have 


F(R) =— (Ry) +c Lalty (RR) 
(R) an Dee male: 


+f Rein Al On Pere ler. se], 
0 


“Rime RL BR} RoR RRA” HERI) 6 
aie ete ee (20) 
Substituting this in (11) and remembering that 
1—R,/R=—64, 
and 
R=R,/(1+5)=R,(1—3-++82—38 . . .), 
we have 
cR, Rokk (—8)? 
D=— 0 8 0 0 2 3 
R | )+e BE R(l1—3 +8 idee (255 
2 2K RyRy 3K] )(—3) 
“ — o(1—3 +8 
ee rab Ry Re hoe R Ro i rai aaa 
or, after a few reductions, 
Cc heaen 
D=$+c++ 28 
hy as 2h? 


Cc hi +h, e [h, is 3 
-{5 ene 28 a teat Pm 08 (21) 
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This is the luminosity distance-red-shift relation, correct to the third order 


in 6. The plus sign in the last term is used for elliptic spaces and the minus sign 
for hyperbolic ones. 


Flat Spaces 
In this case the luminosity distance is given by 


Peep is oR 
Dawe = i.e. B-0| S ares (22) 
Then, firstly, 
dF C 
(iz), ee ee (23) 
secondly, 
2 
aa =A ee eee: (24) 
thirdly, 
bal =-( 2 2h 3K aa we t5'2 5) 
CE) De a ieee a cpa ot he ae 


It is evident, by comparison of (14), (18), (19) and (23), (24), (25) that 
the luminosity distance—red-shift relation is 


¢ hi+h ape hs Ce 
Dory ees ke ae :) ta | 
inl 2h 3 iid ene 
sey Ais co (26) 
where «=—-+1 for an elliptic space, 
a=0 for a flat space, 
a—=-—1l1 for a hyperbolic space, 
and 
hy=Ro/Ro- 
The first two terms in (26) give 
c hoch, 3 
See Ont) Pak ii omhada lero wiaiane 27 


which is the formula commonly used in relativistic cosmology and which does 
not distinguish between the three types of spaces. 


STUDY OF THE FLARE-SURGE EVENT OF SEPTEMBER 7, 1958* 
By Karen M. Lowmant and D. E. Brinest 


Giovanelli and McCabe (1958) indicate that, while most flare-surges appear 
dark against the solar disk throughout the greater part of their existence, it is not 
uncommon for a surge in its earliest stages, just as it leaves the flare, to be bright. 
Rapidly moving material that remains bright against the solar disk for a longer 
time is much less frequently observed. Of the cases of longer lasting bright 
surges that have been reported as exceptional events by Dodson and Hedeman 
(1949, 1952), Ellison (1950), Dodson and McMath (1952), Bray, Loughhead, 
Burgess, and McCabe (1957), and Reid (1959), those reported by Dodson and 
McMath, by Bray et al., and by Reid appear to be detached and rapidly moving 
portions of flares, whereas the remainder did not attain flare brightness. 


The surge event of September 7, 1958 remained bright against the solar 
disk for several minutes, but was clearly less bright than the flare from which it 
appeared to originate. A somewhat unusual programme of filming rate and 
exposure times used by the Climax observers contributed significantly to the 
detail with which the progress of the surge could be followed. During the early 
stages of the surge, four photographs were made per minute with an exposure 
time adapted to bring out the features on the solar disk, while a fifth exposure 
each minute was sufficiently long to show prominence features beyond the limb. 
About 4 min after flare maximum the exposure rate was increased to two frames 
every 2 sec, one timed to bring out disk features, the other to bring out limb 
features. There were several frames that were exposed for an extra-long time, 
probably as a result of faulty shutter operation. On these very long exposures 
faint features in the surge could be detected that otherwise would have gone 
unnoticed. 

Immediately after flare maximum, bright surge material was observed 
tracing a path across the Sun’s disk, from the edge of the flare to the limb, then 
beyond the limb. Some material seemed to part from the main mass and drift 
more nearly radially toward the limb, while the main body of the surge intersected 
a tangent to the limb at an angle of approximately 50°. 


In order to study the motion of the surge, we measured the radial distance 
of the leading edge of the surge from the limb, both before and after it crossed the 
limb. These measurements of radial distance, plotted against time, are shown 
in Figure 1, the velocity in the early stages being 470 km/sec. 


* The research reported in this paper was supported by the National Academy of Sciences as 


part of the International Geophysical Year Programme, with the assistance of a grant from the 
Ford Foundation. 


} Manuscript received December 23, 1959. 
{ The High Altitude Observatory of the University of Colorado, Boulder, Colorado, U.S.A. 
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The surge crossed the limb at 1452 U.T. through a low prominence which 
had appeared on the limb about 1 hr earlier. 


29 sketches of the leading edge of the surge, chosen at 30 sec intervals, gives a 
representation of the trajectory of the surge. 
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Fig. 1.—Variation of surge position with time ; surge of September 7, 1958, at P.A. 72°. 
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surge was a more diffuse and less well defined, slightly slower-moving mass of 
material which rose even higher, to over 100,000 km beyond the limb. 

At about 1454 : 30-1455 U.T., 24-3 min after the first surge crossed the 
limb, a second, somewhat larger and more spectacular surge began building up 


Fig. 2—Composite sketch showing the flare at 1449 U.T. (black area) and the 
trajectory of the first surge from 1449 U.T. until 1502 U.T., with positions of 
the leading edge drawn at 30 second intervals. 
exposures ; the dotted outline from longer exposures. 


The solid outline is from normal 


gradually beyond the limb, slightly north of the first, and rising essentially 
perpendicular to the limb. 

The initial motion of this second surge was slow and the surge itself difficult 
to distinguish from a pre-existing prominence on the limb. In its earlier phase 
this surge was arch-shaped, with a strikingly bright knot of emission located 
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within the centre of the top of the arch. The arch gradually took on a more 
columnar shape with an enlarged top containing the bright knot. This ‘‘ column ”’ 
rose between the above-mentioned prominence and the first surge. As it rose 
the arch-like characteristics of the lower levels disappeared and the prominence 
assumed its former outline on the limb. 

Six minutes after the first suggestion of this slow-starting second surge, 
some very fast material seemed to break through it and speed upward with a 
velocity (between 1500 : 30 and 1502 : 19 U.T.) of nearly 520 km/sec. Thereafter, 


fay 


KG 


Fig. 3.—Composite sketch of the leading edge of the second surge at 60 second 
intervals, between 1452 and 1506 U.T. Darkened outline shows the initial 
stable prominence. 


at about 130,000 km from the limb, it slowed down but continued to move, 
ultimately reaching a distance of more than 140,000km. A few extra-long 
exposed frames of the film show material in the region between 170,000 and 
200,000 km beyond the limb. Figure 3 is a composite of 15 sketches of the 
leading edge of the second surge, at 1-min intervals. 

We next examined the flare around 1455 : 30 U.T., which should have been 
the time that material would have left the flare, at 520 km/sec, in order to appear 
at the limb at 1458:15 U.T. However, no noticeable special event within the 


TABLE 1 
DEDUCED ELECTRON DENSITIES 


Heightvabove: the Electron Density 
Chromosphere (em=*) B 
(km) 
40,000 + 10,000 4-2x 108 
70,000 2b OL0s 
80,000 129 108 


flare could be found on the photographs. The only direct evidence found to 
connect this second surge with the flare were slight traces of dark and bright 
material apparently extending from the flare into the base of the surge. 

Of the various solar radio noise emission features that were reported as 
beginning about the time of the maximum phase of the flare by Boulder, Cornell, 
Ottawa, and Fort Davis, an unclassified burst of intensity 1, having some features 
of a type II burst, and lasting from 1449 -7-1455 U.T.  soneaeee to be the most 
likely to be associable with either surge. This burst was recorded at Fort Davis 
in the 100-580 Me/s frequency range, along with a small group of moderate 
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intensity type III bursts at 1449-50 U.T. of the type described by Giovanelli 
(1958), and an unclassified burst of slight intensity at 1457 U.T. 


By comparing Figure 1 with a reproduction of the Fort Davis sweep- 
frequency record kindly furnished us by Dr. Alan Maxwell, we have been able to 
deduce a rough correspondence between the optically determined position of the 
front of the first surge and the frequency of maximum emission in the simul- 
taneous unclassified burst. Then, by equating the second harmonic of plasma 
frequency to the observed burst frequency, in accordance with the March 1, 
1957 observation of Giovanelli and Roberts (1958), and estimating the height 
of the surge above the chromosphere by the distance, measured normal to the 
limb, from the centre of the flare to the front of the surge, we deduce the 
electron densities given in Table 1. 

These densities, and the corresponding scale height of approximately 
56,000 km, seem quite reasonable for a region in the corona overlying a large 
amount of solar activity. 
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DISTRIBUTION OF AURORAS IN THE SOUTHERN HEMISPHERE* 
By F. R. Bonpf and F. JACKAT 


Several attempts to delineate the southern zone of maximum frequency of 
visibility of auroras have led to widely discrepant results. That of White and 
Geddes (1939) suffered through lack of observations ; that of Vestine and Snyder 
(1945) was based on observed relations between morphology of geomagnetic 
disturbance and auroras in the northern hemisphere, and application to the 
southern hemisphere using very limited geomagnetic data. It may be, too, 
that the distribution of auroras changes slowly with solar activity. 

The data obtained in Antarctica during the International Geophysical Year 
1957-58 offer the possibility of a new study of the geographic distribution of 
auroras around the epoch of solar activity maximum. To be entirely satisfactory 
this study must take into account a number of factors not yet considered even 
in relation to the now extensive body of northern hemisphere data. The 
geographic distribution is different for different auroral forms ; there is a marked 
diurnal variation in latitude and a strong dependence of latitude on level of 
world-wide magnetic disturbance (Jacka 1953). It is also desirable that the 
distribution of auroras be specified with a resolution in latitude much less than 
the width of the curve of frequency of occurrence versus latitude. 


As a first step in such a study I.G.Y. data from the Australian National 
Antarctic Research Expeditions’ (ANARE) stations Mawson and Macquarie 
Island have been used to make a tentative estimate of the location of the isoaurore 
of maximum frequency of occurrence of auroras with discernible lower borders. 


From measured angles of elevation and azimuth to points on the auroral 
lower borders, assumed at a constant height of 105 km, the positions of all 
auroras were plotted on synoptic maps at 15 minute intervals. From these 
maps the frequency of occurrence versus geomagnetic latitude was determined 
on a number of geomagnetic meridians. The frequency distributions all have a 
mode within a few degrees of the observer’s latitude. Because some faint 
auroras observable overhead may not be observable at lower elevations it is 
possible that the observed modal latitude is shifted a little towards the observer’s 
latitude. We have therefore taken as our estimate of the latitude of the mode, 
the mid point of the range in which the frequency exceeds 2/3 the maximum 
observed frequency. These estimated positions are marked by crosses in 
Figure 1. 

These data are very meagre, so we tentatively introduce the hypothesis 
that the isoaurores are symmetrical about the geomagnetic meridian through 
the magnetic dip pole. Guided by the Mawson and Macquarie Island data and 


* Manuscript received April 4, 1960. 
j Antarctic Division, Department of External Affairs, Melbourne. . 
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this hypothesis the full line curve on Figure 1 was drawn as a tentative repre- 
sentation of the isoaurore of maximum frequency of occurrence of auroras with 
discernible lower borders. 
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Fig. 1.—Full line curve : estimate of the isoaurore of maximum frequency of occurrence of auroras 

with discernible lower borders based on observed positions denoted by crosses in the regions of 

Macquarie Island and Mawson and assuming symmetry about the geomagnetic meridian through 

the dip pole. Dashed curve: Hultqvist’s projection from 5-6 Earth radii in the geomagnetic 
equatorial plane. 


The cross near geographic longitude 345° H. indicates the position of most 
frequent occurrence of auroral arcs estimated from the data summarized by 
Bvans and Thomas (1959) for Halley Bay. This point is in reasonable agreement 


with our curve. 
If Figure 1 is not grossly in error, there is one other I.G.Y. station (Byrd) 


from which the maximum frequency isoaurore should be observed at high 
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elevation and several (marked by squares on Figure 1) from which it might be 
observed at low elevation and consequently low reliability. 

The synoptic maps also permit the measurement of orientation of auroras 
relative to the meridians. The orientation angles show a small diurnal variation 
but the frequency distribution of angles is fairly narrow and roughly symmetrical. 
The mean orientations of ares and simply shaped bands on a number of meridians 
near Mawson and Macquarie Island are indicated by short double-ended arrows 
on Figure 1. These are approximately parallel with the estimated isoaurore. 
We expect then that the mean value of azimuth (taking values 180° apart as 
equivalent) to the nearest (greatest elevation) point of arcs is perpendicular to 
the isoaurores. Two such azimuths (from data of Hatherton and Salmon (1959)) 
are indicated by arrows from Scott Base and Hallett in Figure 1. Agreement 
with our estimated isoaurore is reasonable. 


In the region well inside the maximum frequency isoaurore, even the general 
form of the geographical distribution of auroras is uncertain. There is, however, 
mounting evidence for the existence of an inner zone of maximum frequency. 


At this stage it is perhaps desirable to seek some correspondence with 
northern hemisphere data. Hultqvist (1958, 1959) has shown that projections 
from circles in the geomagnetic equatorial plane, along the actual magnetic field 
lines, onto the Earth’s surface, produce oval curves which very closely parallel 
the auroral isochasms in the northern hemisphere. Using Hultqvist’s data, we 
find that a circle of 5-6 Earth radii projects onto the dashed curve in Figure 1. 
The discrepancy between the curve and our maximum frequency isoaurore is 
so great as to cast doubts on the soundness of this approach to the problem of 
revealing a correspondence between the two hemispheres. On the other hand, 
Hultqvist’s projections are based on geomagnetic field data of epoch 1945 which 
were very scanty in high latitudes, especially in the southern hemisphere. 


The authors are now seeking data from other antarctic stations to enable 
more comprehensive study of the distribution and morphology of auroras in the 
southern hemisphere. 
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THE STRESS FIELDS AROUND SOME DISLOCATION ARRAYS* 
By A. K. Heapt 


The stress fields around some edge dislocation arrays have been calculated 
using the approximation of Leibfried (1951) that the discrete dislocations are 
replaced by a continuous dislocation density. Haasen and Leibfried (1954) have 
calculated one such stress distribution (corresponding to (ii) below) via the stress 
function which they evaluated as a superposition integral over the dislocation 
distribution of the stress function of an infinitesimal dislocation. The results 
given here were evaluated directly as the superposition of the stress fields of the 
infinitesimal dislocations making up the distribution. Like most of the integrals 
which occur in this approximation, they can be evaluated most readily by contour 
integration (Haasen and Leibfried 1954). A multiple coordinate system is used 
and is such that the results can be written compactly and approximations and 
limiting cases can be easily obtained. 

(i) Blocks at z=-+La with n positive dislocations between «=b and x=a 
and n negative dislocations between «=—b and «=—a which are held apart by 
an applied shear stress o (Head and Louat 1955 (vi)). 


The components of stress at a point P(x,y) are 
27 Ry 
gamo( 208) rs Op cos in a (9, +0,+8,+0, )} 
uf 


3 
aa as cos {yn +3(91+ 2) — 5 (91 +92)} 


UL] 


+ 2 sin Mert e—0—03)|, YO a (1) 


e2\#{ Ry 

ay =0(Pe] fr 212 cos {y—4F(91 + 2 +8; +9,)} 
3 

a a cos {4 +3(91 +2) — 5101-00} + (2) 

3 Ry 

Oxy =0( S62) | — =% sin {n—3(p1 +92 +61 +69)} 

a P1P2 


3 
+ i sin (n-+Hes+92)— 9(01+83)} 
1 
+ cos $(9,;+9.—0, —0 »} hte g historian (3) 


where (R,y), (71,91), (12,92), (P1sP1)) (P22) are polar coordinates, as shown in 
Figure 1. The expression (3) for the shear stress is the total shear stress including 


the applied shear stress o. 


* Manuscript received February 12, 1960. 
+ Division of Tribophysics, C.S.I.R.0., Melbourne. 


614 SHORT COMMUNICATIONS 


(ii) Haasen and Leibfried (1954) have considered the special case of (i) 
for b—0. This can be derived from (i) by putting R=e;=p2, 7=?i1=P2 and 
then (1)-(3) are equivalent to their result. 


-a 


Fig. 2.—Coordinate system for equations (4)—(6). 


Fig. 3.—Coordinate system for Section (iv). 


(ili) If the negative dislocations in (i) are removed to infinity by putting 
T3=02—2R=co, 9, =y=—0,=0, we are left with Figure 2, a group of positive 
dislocations forced against a block by a shear stress o (Head and Louat 1955 (iv)), 
and (1)-(3) become 


c..=balelr)|f cos (0+) — " cos $(@ —30) +4 sin He-O, .. (4) 


Cwy=Lo(oin)| — Ecos H0+9)+ Leos He —80)], eee (5) 


oy =holelr)#| 2 sin 4(0+¢) pe sin $(@ —30) +2 cos He-0). se 6) 


The stresses near the head of the pile-up will be given approximately by 
putting 
oz DL, 90, ral, y= ino 
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(where L is the length of the pile-up). Equations (4)-(6) then become the 
approximate expressions derived by Stroh (1954) from consideration of a pile-up 
of discrete dislocations. 

(iv) Figure 3 is example (vii) from Head and Louat (1955) with » positive 
dislocations beyond «=—a and n negative dislocations beyond vw=a, driven 
together by a shear stress oc but prevented from coalescing by blocks at #=-La. 
The stress field in this case is also given by (1)-(3) if the symbols are reinterpreted 
as in Figure 3. 
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THE OPTIMUM LINE WIDTH FOR A REFLECTION CAVITY MASER* 
By G. J. TROUPT 


Introduction 

It has been pointed out previously (Troup 1959, hereinafter referred to as I) 
that there is an optimum line width for the transition used in a reflection cavity 
maser and that this optimum line width gives maximum bandwidth for a given 
gain. However, it was not explicitly deduced in I that this optimum line width 
actually gives the maximum value of (Gain)? x Bandwidth product for a given 
cavity maser system. This note makes the appropriate deduction. 


The Line Width for Maximum (Gain)* x Bandwidth Product 
Commencing with equation (3) of I, we have that the gain @ of a reflection 
cavity maser is given to a good approximation by 


GA (| Om |+Q)z| Om |—Qr)7) we eee (1) 
where @, is the loaded Q of the cavity and | Q,, | is the modulus of the negative 
molecular Q. 

The bandwidth B of a reflection cavity maser, assuming a Lorentz line shape 
for the amplifying transition, is given by (equation (4) of I) 

B=d[(| Qn (/Or)—TIIA +8) On LT ee ee (2) 
where 3 is the line width of the amplifying transition and fis the centre frequency 
of both cavity and molecular responses. 

Using equations (1) and (2), we get 
GB~B[(| Qn Ox) +148 | Om lift eee (3) 


* Manuscript received April 29, 1960. 
+ Weapons Research Establishment, Salisbury, S. Aust. 
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Now for reasonably large gains (@>10 say) | Q,,|~@Q,, so that 


GIB~28(118 | Q,, |[f)mt. eee cee ee ee ec ee eee (4) 
Now the molecular Q, Q,,, is of the form (equation (2) of I) 
| Q,5.( OLN, oa os PR ee (5) 


where K is a constant for the particular active material and cavity system, and 
N* is the number density of excess molecules in the emissive state. 
Using equation (5), we rewrite equation (4) as 
GiB 238[ 14-6 fh) a ee ee ve (6) 

Differentiation of this expression with respect to 5 shows that G*B has a maximum 
value 

(GIB) pas. = (Nf) ee ee ee (7) 
when S8=(K/N*f)-!=Sopt.. The quantity Sop. is identical with the value of 6 
deduced in I to give maximum bandwidth at a given gain. This means that 
suitable adjustment of the transition line width 6, keeping all other parameters 
constant, can maximize the (Gain)? x Bandwidth product. 


Conclusion 

In the paramagnetic maser, the simplest way of adjusting 6 is clearly to use 
a magnetic field deliberately made inhomogeneous. This of course requires that 
5 <Sopt. itially. The discovery of the narrowness of the lines of Fe** in Ti,O, 
(Butcher and Gill 1959; Low 1960) makes the technique described above 
particularly applicable to this material. Moreover, Townes e¢ al. at Columbia 
have obtained a (Gain)? x Bandwidth product of 250 Me/s at 4500 Me/s using 
ruby (Cr?+ in Al,O,) as the active medium, by the use of magnetic field inhomo- 
geneity (Butcher and Gill 1959). This latter result is taken as experimental 
confirmation of the existence of an optimum line width for a given maser system. 
It also confirms that the use of an inhomogeneous magnetic field is a practicable 
way of adjusting the line width. 


Note added in Proof—Maiman (1960) has come to the same conclusions 
regarding optimum line width, and has pointed out some of the consequences for 
paramagnetic masers. 


The author is indebted to Dr. J. E. Houldin, Physics Department, Chelsea 
Polytechnic, London, for encouraging the submission of this note for publication. 
Permission to publish has been granted by the Chief Scientist, Research and 
Development Branch, Department of Supply, Australia. 
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PHOTOPROTONS FROM TANTALUM* 
By J. H. Carver,t RK. B. Taytor,t and W. TURCHINETzt 


Photoproton emission from tantalum has been studied by measuring the 
yield of the 5-5 hr isomer of 1®Hf produced in the 1®!Ta(y,p)!®™Hf reaction. 
The experimental arrangements were similar to those used in a previous investi- 
gation (Carver and Turchinetz 1958) of the photodisintegration of tantalum, 
except that the shielding of the scintillation counter, a 14 in. diameter by 2 in. 
long NalI(T1) crystal, has now been improved along the lines described by Pringle, 


3 


N 


; COUNTS/CHANNEL (ARBITRARY UNITS) 


200 (KEV) 


12} 10 20 30 40 5c 
CHANNEL 


Fig. 1.—y-Ray spectra of 5-5 hr Hf and 2-1 hr *°Ta, Spectrum B was 
measured 20 hr after spectrum A and corresponds mainly to '°Hf. 


Turchinetz, and Funt (1955). Larger quantities of tantalum (65 g) were used 
-in the present experiments and, with these improvements, it was possible to 
separate the 5-5 hr 1®Hf activity from the similar, and more abundant, 2:1 hr 
1784, activity produced in the 18!Ta(y,3n)!78Ta reaction. 

In the energy region from 200 to 500 kV 1*°Hf has y-ray lines at 216, 332, 
443, and 501 keV, and 178Ta has lines at 214, 326, 332, and 427 keV (Strominger, 


* Manuscript received May 23, 1960. 
+ Research School of Physical Sciences, Australian National University, Canberra. 
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Hollander, and Seaborg 1958). At lower energies the y-ray spectrum is dominated 
by 8:15 hr 18Ta produced in the 1*!Ta(y,n)'*Ta reaction. The difference 
between the y-ray spectra of 4Hf and 176Ta is illustrated in Figure 1, which 
shows two of the observed spectra. Spectrum B of Figure 1 was measured 20 hr 
after spectrum A and corresponds mainly to 1%°Hf. An analysis of the decay of 
these y-ray lines, an example of which is shown in Figure 2, indicates that pee 
are two components with half-lives of 5-5 and 2-1 hr, corresponding to 1%Hf 
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Fig. 2.—Decay curve for the 5-5 hr 18Hf and 2-1 hr 178Ta y-rays 
of Figure 1. 


and *’*Ta respectively. It is interesting to note that the half-life of !78Ta obtained 
in this way is the same as that measured by Wilkinson (1950), whereas if no 
correction is made for the 5-5 hr component a value of 2-5 hr is obtained 
(Carver and Turchinetz 1958). In this analysis a small correction has been 
made for the bremsstrahlung produced by the B-particles of 8-15 hr 18°Ta, which 
is produced in the 18!Ta(v,n) reaction. This correction was made by irradiating 
a tantalum sample at 15 MeV, where there is no yield of 17°Ta or 18Hf, and 
normalizing the bremsstrahlung tail to the intensity of the 55 keV X-ray peak. 
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The measured yield function and the derived cross-section curve for the 
*8!Ta(y,p)'8™Hf reaction are shown in Figure 3. The points on the yield function 
were obtained from a least squares analysis of the 5-5 hr Hf plus 2-1 hr 
™8Ta y-rays. The absolute cross-section scale was derived by comparison of 
these yields with the 10 min ®Cu activity induced in copper foils using the 
absolute °°Cu(y,n) cross-section data of Berman and Brown (1954). This method 
is similar to that previously used to determine the (,n), (y,2n), and (y,3n) cross 
sections in 1§!Ta (Carver and Turchinetz 1958). 
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Fig. 3.—(a) The yield curve and (b) the derived cross section 
for the 1°Ta(y,p)!®°™Hf reaction. 


The threshold for the 18!Ta(y,p) reaction is 6-7 MeV (Wapstra 1958), but 
owing to the effect of the Coulomb barrier there is no appreciable yield below a 
y-ray energy of ~16MeV. The 1*!Ta(y,p) reaction leading to the 5-5 hr 
isomer of Hf has a peak cross section of 1-0 mbn at a y-ray energy of 28 MeV 
and an integrated cross section to 32 MeV of (10 +2) MeV mbn* which corresponds 
to 0-3 per cent. of the total integrated y-ray absorption cross section. Owing 
- to the difference between the shapes of the (y,p) and the total y-ray cross section 
the yield for the (,p) reaction leading to the isomeric state is 0-1 per cent. of the 
total yield for all y-ray processes. Most of the previous data about (y,p) cross 


* Approximately equal contributions to the error come from the analysis of the decay curves 
and from the comparison of the 1®Hf plus 17*Ta yields with the Cu yield. There is a smaller 
error, 6 per cent., in Berman and Brown’s (1954) estimate of the integrated ®*Cu(y,n) cross section 


on which the present determination is based. 
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sections in the heavy nuclei have been obtained with 23 MeV bremsstrahlung 
(e.g. Toms and Stephens 1955). The present results, together with those obtained 
for caesium and iodine by Taylor (1960) in this laboratory, indicate that higher 
energy y-rays make a considerable contribution to the photoproton emission. 


Toms and Stephens (1955) used photographic plates to measure the total 
photoproton yield from tantalum for 23 MeV bremsstrahlung. Combining their 
result with the present measurements suggests that at 23 MeV the 5-5 hr isomeric 
state is populated in ~10 per cent. of the (y,p) disintegrations. The spin of the 
5-5 hr level of 18°Hf is probably 9 and the ground states of ®°Hf and +*’Ta have 
spins 0 and 7/2 respectively (Strominger, Hollander, and Seaborg 1958), so that 
the present results are consistent with the general rule that in the production 
of an isomeric pair the member with spin closest to that of the target nucleus is 
favoured. If the bulk of the (y,p) reactions proceeds through direct transitions 
leading to states in 1®Hf of only moderate excitation it is not surprising that the 
magnitude of the preference for the low spin state is larger than in the 
181Ta(v,3n)'78Ta, 178mTa reactions, where the high spin state is populated in 
30 per cent. of the disintegrations (Carver and Turchinetz 1958). 
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CORRIGENDA 
VotumE 13, Numper 2A, Paces 451-9 


‘* Susceptibility of Chromium-based Alloys of Transition Elements ” 
Pages 451, 453, and 458: For Neel read Néel. 
Page 453, line 20: The expression within the braces should read: 
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Page 453, line 22: Insert as proportional before to. 
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